The quasispecies nature and biological implications of the hepatitis C virus

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.

Impact of Intrahost NS5 Nucleotide Variations on Dengue Virus Replication

Due to the nature of RNA viruses, their high mutation rates produce a population of closely related but genetically diverse viruses, termed quasispecies. To determine the role of quasispecies in DENV disease severity, 22 isolates (10 from mild cases, 12 from fatal cases) were obtained, amplified, and sequenced with Next Generation Sequencing using the Illumina MiSeq platform. Using variation calling, unique wildtype nucleotide positions were selected and analyzed for variant nucleotides between mild and fatal cases. The analysis of variant nucleotides between mild and fatal cases showed 6 positions with a significant difference of p < 0.05 with 1 position in the structural region, and 5 positions in the non-structural (NS) regions. All variations were found to have a higher percentage in fatal cases. To further investigate the genetic changes that affect the virus’s properties, reverse genetics (rg) viruses containing substitutions with the variations were generated and viral growth properties were examined. We found that the virus variant rgNS5-T7812G (G81G) had higher replication rates in both Baby hamster kidney cells (BHK-21) and Vero cells while rgNS5-C9420A (A617A) had a higher replication rate only in BHK-21 cells compared to wildtype virus. Both variants were considered temperature sensitive whereby the viral titers of the variants were relatively lower at 39°C, but was higher at 35 and 37°C. Additionally, the variants were thermally stable compared to wildtype at temperatures of 29, 37, and 39°C. In conclusion, viral quasispecies found in isolates from the 2015 DENV epidemic, resulted in variations with significant difference between mild and fatal cases. These variations, NS5-T7812G (G81G) and NS5-C9420A (A617A), affect viral properties which may play a role in the virulence of DENV.

Towards a Systems Immunology Approach to Understanding Correlates of Protective Immunity against HCV

Over the past decade, tremendous progress has been made in systems biology-based approaches to studying immunity to viral infections and responses to vaccines. These approaches that integrate multiple facets of the immune response, including transcriptomics, serology and immune functions, are now being applied to understand correlates of protective immunity against hepatitis C virus (HCV) infection and to inform vaccine development. This review focuses on recent progress in understanding immunity to HCV using systems biology, specifically transcriptomic and epigenetic studies. It also examines proposed strategies moving forward towards an integrated systems immunology approach for predicting and evaluating the efficacy of the next generation of HCV vaccines.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

A SARS-CoV-2 host infection model network based on genomic human Transcription Factors (TFs) depletion

In December 2019 a new beta-coronavirus was isolated and characterized by sequencing samples from pneumonia patients in Wuhan, Hubei Province, China. Coronaviruses are positive-sense RNA viruses widely distributed among different animal species and humans in which they cause respiratory, enteric, liver and neurological symptomatology. Six species of coronavirus have been described (HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1) that cause cold-like symptoms in immunocompetent or immunocompromised subjects and two strains of sometimes fatal zoonotic origin that cause severe acute respiratory syndrome (SARS-CoV and MERS-CoV). The SARS-CoV-2 strain is the emerging seventh member of the coronavirus family, which is actually determining a global emergency.

In silico analysis is a promising approach for understanding biological events in complex diseases and due to serious worldwide emergency and serious threat to global health, it is extremely important to use bioinformatics methods able to study an emerging pathogen like SARS-CoV-2.

Herein, we report on in silico comparative analysis between complete genome of SARS-CoV, MERS-CoV, HCoV-OC43 and SARS-CoV-2 strains, to identify the occurrence of specific conserved motifs on viral genomic sequences which should be able to bind and therefore induce a subtraction of host's Transcription Factors (TFs) which lead to a depletion, an effect comparable to haploinsufficiency (a genetic dominant condition in which a single copy of wild-type allele at a locus, in heterozygous combination with a variant allele, is insufficient to produce the correct quantity of transcript and, therefore, of protein, for a correct standard phenotypic expression).

In this competitive scenario, virus versus host, the proposed in silico protocol identified the TFs same as the distribution of TFBSs (Transcription Factor Binding Sites) on analyzed viral strains, potentially able to influence genes and pathways with biological functions confirming that this approach could brings useful insights regarding SARS-CoV-2. According to our results obtained by this in silico approach it is possible to hypothesize that TF-binding motifs could be of help in the explanation of the complex and heterogeneous clinical presentation in SARS-CoV-2 and subsequently predict possible interactions regarding metabolic pathways, and drug or target relationships.

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.

Using Next Generation Sequencing to Study the Genetic Diversity of Candidate Live Attenuated Zika Vaccines

Zika virus (ZIKV) is a mosquito-transmitted positive-sense RNA virus in the family Flaviviridae. Candidate live-attenuated vaccine (LAV) viruses with engineered deletions in the 3’ untranslated region (UTR) provide immunity and protection in animal models of ZIKV infection, and phenotypic studies show that LAVs retain protective abilities following in vitro passage. The present study investigated the genetic diversity of wild-type (WT) parent ZIKV and its candidate LAVs using next generation sequencing analysis of five sequential in vitro passages. The results show that genomic entropy of WT ZIKV steadily increases during in vitro passage, whereas that of LAVs also increased by passage number five but was variable throughout passaging. Additionally, clusters of single nucleotide variants (SNVs) were found to be present in the pre-membrane/membrane (prM), envelope (E), nonstructural protein NS1 (NS1), and other nonstructural protein genes, depending on the specific deletion, whereas in the parent WT ZIKV, they are more abundant in prM and NS1. Ultimately, both the parental WT and LAV derivatives increase in genetic diversity, with evidence of adaptation following passage.

HCV genotype profile in Brazil of mono-infected and HIV co-infected individuals: A survey representative of an entire country

INTRODUCTION

To be eligible for government-provided treatment in Brazil, all HCV-infected individuals are required to be genotyped shortly after diagnosis. We describe the HCV genotype (G) profiles by geographic region, gender, age and HIV co-infection.

METHODS

We assessed 29,071 genotypes collected from HCV-infected individuals from March 2016 to March 2018 (Abbott Real-Time HCV Genotype). We randomly selected 12,336 samples for HIV co-infection testing using an EIA rapid test kit (TR DPP HIV 1/2 Bio-Manguinhos). Descriptive statistical analyses were performed using R.

RESULTS

Overall, HCV genotype distribution was 40.9% G1A, 30.2% G1B, 23.8% G3, 3.8% G2, 0.7% G4, 0.1% G5 and 0.6% with multiples genotypes. G1A prevalence was 44.4% among males and 35.8% among females. G1B and G2 were more prevalent in older individuals than G1A and G3. G3 was more prevalent in the South region. Of samples tested for HIV co-infection, 15% were HIV+. Median age among HCV/HIV co-infected individuals was 50 years old compared to 57 years old among mono-infected individuals. Distinct HCV genotype prevalence between HCV/HIV co-infected and HCV mono-infected individuals were respectively: G1A 60.6% versus 37.8%, G1B 15.2% versus 32.9%, and G3 18.9% versus 24.7%. G4 was detected among co-infected young men (3.5% versus 0.2% among mono-infected).

CONCLUSION

The increasing prevalence of G3, as inferred by the younger ages of the HCV-infected individuals, poses an extra challenge with regards to disease progression. Distinct genotypical profiles between HCV mono-infection and HCV/HIV co-infection warrant future research in order to better understand and help mitigate HCV chains of transmission.

Inter- and intra-lineage genetic diversity of wild-type Zika viruses reveals both common and distinctive nucleotide variants and clusters of genomic diversity

Zika virus (ZIKV) strains belong to the East African, West African, and Asian/American phylogenetic lineages. RNA viruses, like ZIKV, exist as populations of genetically-related sequences whose heterogeneity may impact viral fitness, evolution, and virulence. Genetic diversity of representative ZIKVs from each lineage was examined using next generation sequencing (NGS) paired with downstream entropy and single nucleotide variant (SNV) analysis. Comparisons showed that inter-lineage diversity was statistically supported, while intra-lineage diversity. Intra-lineage diversity was significant for East but not West Africa strains. Furthermore, intra-lineage diversity for the Asian/American lineage was not supported for human serum isolates; however, a placenta isolate differed significantly. Relative in the pre-membrane/membrane (prM/M) gene of several ZIKV strains. Additionally, the East African lineage contained a greater number of synonymous SNVs, while a greater number of non-synonymous SNVs were identified for American strains. Further, inter-lineage SNVs were dispersed throughout the genome, whereas intra-lineage non-synonymous SNVs for Asian/American strains clustered within prM/M and NS1 gene. This comprehensive analysis of ZIKV genetic diversity provides a repository of SNV positions across lineages. We posit that increased non-synonymous SNV populations and increased relative genetic diversity of the prM/M and NS1 proteins provides more evidence for their role in ZIKV virulence and fitness.

Kinetics of hepatitis B surface antigen quasispecies during REP 2139-Ca therapy in HBeAg-positive chronic HBV infection

Chronic HBV infection results in various clinical manifestations due to different levels of immune response. In recent years, hepatitis B treatment has improved by long-term administration of nucleos(t)ide analogues (NUCs) and peg-interferon. Nucleic acid polymers (NAPs; REP 2139-Ca and REP 2139-Mg) are new antiviral drugs that block the assembly of subviral particles, thus preventing the release of HBsAg and allowing its clearance and restoration of functional control of infection when combined with various immunotherapies. In the REP 102 study (NCT02646189), 9 of 12 patients showed substantial reduction of HBsAg and seroconversion to anti-HBs in response to REP 2139-Ca, whereas 3 of 12 patients did not show responses (>1 log reduction of HBsAg and HBV DNA from baseline). We characterized the dynamic changes of HBV quasispecies (QS) within the major hydrophilic region (MHR) of the 'pre-S/S' open reading frame including the 'a' determinant in responders and nonresponders of the REP 102 study and four untreated matched controls. HBV QS complexity at baseline varied slightly between responders and nonresponders (P = .28). However, these responders showed significant decline in viral complexity (P = .001) as REP 2139-Ca therapy progressed but no significant change in complexity was observed among the nonresponders (P = .99). The MHR mutations were more frequently observed in responders than in nonresponders and matched controls. No mutations were observed in 'a' determinant of major QS population which may interfere with the detection of HBsAg by diagnostic assays. No specific mutations were found within the MHR which could explain patients' poor HBsAg response during REP 2139-Ca therapy.

Associations between serum HBX quasispecies and their integration in hepatocellular carcinoma

HBV quasispecies are closely related to the course and outcome of liver disease. However, whether the complexity and diversity of HBX quasispecies affects its integration in the liver cell and thereby enhances the resultant carcinogenesis is still not clear. 15 HCC patients were recruited; genomic DNA and HBV DNA were extracted from liver cancer tissue and serum respectively. The integrated HBX fragment in liver cancer tissue was amplified by Alu repeat sequence-polymerase chain reaction (Alu-PCR) and sequenced. The serum HBX gene was amplified by nested PCR and sequenced. Quasispecies complexity and diversity, phylogenetic characteristics, lymphocyte count and survival time between HBX-integrated and HBX-unintegrated patients were evaluated. Results showed that the integrated HBX fragment was detected in the tumor tissue of nine patients, and the integration rate was 60.00% (9/15). Compared with the HBX-unintegrated patients, the HBX-integrated patients had a higher quasispecies complexity (P=0.028 and 0.004, at the nucleotide and amino acid levels, respectively). The HBX-integrated patients had a tendency of higher quasispecies diversity, lower lymphocyte count and the survival time. A total of 12 mutation sites were revealed in the HBX-integrated fragment after alignment with the reference sequence. In these, the HBX-integrated groups had significantly higher mutation frequencies at C1497T, A1630G, G1721A, A1762T/G1764A and A1774G. This study revealed influence factors of HBX integration both in virus and the host. The increased complexity and diversity of HBX quasispecies might destroy the host immune balance, and lead to HBX integration ultimately.

Assessing Species Diversity Using Metavirome Data: Methods and Challenges

Assessing biodiversity is an important step in the study of microbial ecology associated with a given environment. Multiple indices have been used to quantify species diversity, which is a key biodiversity measure. Measuring species diversity of viruses in different environments remains a challenge relative to measuring the diversity of other microbial communities. Metagenomics has played an important role in elucidating viral diversity by conducting metavirome studies; however, metavirome data are of high complexity requiring robust data preprocessing and analysis methods. In this review, existing bioinformatics methods for measuring species diversity using metavirome data are categorised broadly as either sequence similarity-dependent methods or sequence similarity-independent methods. The former includes a comparison of DNA fragments or assemblies generated in the experiment against reference databases for quantifying species diversity, whereas estimates from the latter are independent of the knowledge of existing sequence data. Current methods and tools are discussed in detail, including their applications and limitations. Drawbacks of the state-of-the-art method are demonstrated through results from a simulation. In addition, alternative approaches are proposed to overcome the challenges in estimating species diversity measures using metavirome data.

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.

Serotype-specific interactions among functional domains of dengue virus 2 nonstructural proteins (NS) 5 and NS3 are crucial for viral RNA replication

Four serotypes of mosquito-borne dengue virus (DENV), evolved from a common ancestor, are human pathogens of global significance for which there is no vaccine or antiviral drug available. The N-terminal domain of DENV NS5 has guanylyltransferase and methyltransferase (MTase), and the C-terminal region has the polymerase (POL), all of which are important for 5'-capping and RNA replication. The crystal structure of NS5 shows it as a dimer, but the functional evidence for NS5 dimer is lacking. Our studies showed that the substitution of DENV2 NS5 MTase or POL for DENV4 NS5 within DENV2 RNA resulted in a severe attenuation of replication in the transfected BHK-21 cells. A replication-competent species was evolved with the acquired mutations in the DENV2 and DENV4 NS5 MTase or POL domain or in the DENV2 NS3 helicase domain in the DENV2 chimera RNAs by repeated passaging of infected BHK-21 or mosquito cells. The linker region of seven residues in NS5, rich in serotype-specific residues, is important for the recovery of replication fitness in the chimera RNA. Our results, taken together, provide genetic evidence for a serotype-specific interaction between NS3 and NS5 as well as specific interdomain interaction within NS5 required for RNA replication. Genome-wide RNAseq analysis revealed the distribution of adaptive mutations in RNA quasispecies. Those within NS3 and NS5 are located at the surface and/or within the NS5 dimer interface, providing a functional significance to the crystal structure NS5 dimer.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

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.

Heads or Tails: Genotyping of Hepatitis C Virus Concerning the 2k/1b Circulating Recombinant Form

As different hepatitis C virus (HCV) genotypes respond differently to initiated therapy, correct HCV genotyping is essential. A potential risk for misclassification of the intergenotypic HCV circulating recombinant form (CRF) 2k/1b strains exists, depending on the genotyping method used. The aim was to investigate the differences in HCV genotyping methods with regard to CRF 2k/1b and to gain insight in the prevalence of the CRF 2k/1b. Genotyping results by Versant HCV Genotype Assay were compared with nonstructural protein 5B (NS5B) sequencing. In total, from November 2001 until March 2015, 3296 serum samples were analyzed by Versant HCV Genotype Assay. As misclassified CRF is harbored among HCV genotype 2, we further focused our search on 142 (4.3%) samples positive for HCV genotype 2. On 116 (81.7%) retrieved samples, the NS5B sequencing was performed. Twelve out of the 116 retrieved samples (10.3%) were classified as CRF 2k/1b by sequencing of the NS5B region. Ten of these 12 samples were originally misclassified as genotype 2a or 2c, while 2 of them were misclassified as genotype 2. Our results show that the current prevalence of CRF 2k/1b is underestimated. The importance of correct HCV genotyping is emphasized, considering the tailored choice of treatment regimen and overall prognosis.

Copy number variation of genes involved in the hepatitis C virus-human interactome

Copy number variation (CNV) is a newly discovered form of intra-species genetic polymorphism that is defined as deletions or duplications of genome segments ranging from 1 kbp to several Mbp. CNV accounts for the majority of the genetic variation observed in humans (CNV regions cover more than 10% of the human genome); therefore, it may significantly influence both the phenotype and susceptibility to various diseases. Unfortunately, the impact of CNV on a number of diseases, including hepatitis C virus (HCV) infection, remains largely unexplored. Here, we analyzed 421 human genes encoding proteins that have been shown to interact with HCV proteins or genomic RNA (proteins from the HCV-human interactome). We found that 19 of the 421 candidate genes are located in putative CNV regions. For all of these genes, copy numbers were determined for European, Asiatic and African populations using the multiplex ligation-dependent amplification (MLPA) method. As a result, we identified 4 genes, IGLL1, MLLT4, PDPK1, PPP1R13L, for which the CN-genotype ranged from 1 to 6. All of these genes are involved in host-virus interaction; thus, their polymorphism has a potential impact on the development of HCV infection and/or therapy outcome.

Viral quasispecies complexity measures

Mutant spectrum dynamics (changes in the related mutants that compose viral populations) has a decisive impact on virus behavior. The several platforms of next generation sequencing (NGS) to study viral quasispecies offer a magnifying glass to study viral quasispecies complexity. Several parameters are available to quantify the complexity of mutant spectra, but they have limitations. Here we critically evaluate the information provided by several population diversity indices, and we propose the introduction of some new ones used in ecology. In particular we make a distinction between incidence, abundance and function measures of viral quasispecies composition. We suggest a multidimensional approach (complementary information contributed by adequately chosen indices), propose some guidelines, and illustrate the use of indices with a simple example. We apply the indices to three clinical samples of hepatitis C virus that display different population heterogeneity. Areas of virus biology in which population complexity plays a role are discussed.

Sendai virus intra-host population dynamics and host immunocompetence influence viral virulence during in vivo passage

In vivo serial passage of non-pathogenic viruses has been shown to lead to increased viral virulence, and although the precise mechanism(s) are not clear, it is known that both host and viral factors are associated with increased pathogenicity. Under- or overnutrition leads to a decreased or dysregulated immune response and can increase viral mutant spectrum diversity and virulence. The objective of this study was to identify the role of viral mutant spectra dynamics and host immunocompetence in the development of pathogenicity during in vivo passage. Because the nutritional status of the host has been shown to affect the development of viral virulence, the diet of animal model reflected two extremes of diets which exist in the global population, malnutrition and obesity. Sendai virus was serially passaged in groups of mice with differing nutritional status followed by transmission of the passaged virus to a second host species, guinea pigs. Viral population dynamics were characterized using deep sequence analysis and computational modeling. Histopathology, viral titer and cytokine assays were used to characterize viral virulence. Viral virulence increased with passage and the virulent phenotype persisted upon passage to a second host species. Additionally, nutritional status of mice during passage influenced the phenotype. Sequencing revealed the presence of several non-synonymous changes in the consensus sequence associated with passage, a majority of which occurred in the hemagglutinin-neuraminidase and polymerase genes, as well as the presence of persistent high frequency variants in the viral population. In particular, an N1124D change in the consensus sequences of the polymerase gene was detected by passage 10 in a majority of the animals. In vivo comparison of an 1124D plaque isolate to a clone with 1124N genotype indicated that 1124D was associated with increased virulence.

Mechanisms of Virologic Failure with Direct-Acting Antivirals in Hepatitis C and Strategies for Retreatment

The current standard of care for hepatitis C therapy is the combination of direct-acting antiviral (DAA) agents. These orally administered medications target the viral proteins and halt the hepatitis C virus lifecycle. Despite high cure rates with these novel drugs, virologic failure with DAAs are of mounting concern as real-world sustained virologic response 12 rates seem lower than expected. The mechanisms of virologic failure to DAAs are likely multifactorial, including baseline resistance variants, the efficacy of the agents used, and host factors. Salvage therapy for DAA virologic failures is an area of emerging research.

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.

Impact of IL28B and PNPLA3 polymorphisms on treatment outcomes in patients infected with genotype 6 hepatitis C virus

BACKGROUND/AIMS

Interleukin-28B (IL28B) and patatin-like phospholipase domain containing 3 (PNPLA3) gene polymorphisms are associated with hepatitis C virus (HCV) clearance and fatty liver, respectively. We aimed to test if their polymorphisms are associated with virologic responses in Chinese chronic hepatitis C (CHC) patients.

METHODS

This was a retrospective-prospective cohort study. Consecutive patients infected by genotype 1 and 6 HCV received antiviral therapy were included. Host IL-28B rs12979860/rs8099917 and PNPLA3 rs738409 genotype were tested. The primary outcome was sustained virologic response (sustained virologic response [SVR]: undetectable HCV RNA 24 weeks post-treatment).

RESULTS

From 305 patients had positive antibody to HCV, 52 and 31 patients infected by genotype 1 and 6 HCV, respectively were recruited. Mean age was 58 ± 11 years; 70% were male. Mean baseline HCV RNA was 6.8 ± 2.7 log IU/ml. The SVR for patients infected by genotype 1 and 6 HCV was 67.3% and 90.3%, respectively. The proportions of IL28B genotypes were 78%, 21%, and 1% for TT/TG/GG at rs8099917, and 81%, 18%, and 1% for CC/TC/TT at rs12979860, respectively. The proportions of PNPLA3 rs738409 genotypes were 16%, 36%, and 48% for GG/GC/CC. IL28B genotype was significantly associated with SVR in patients infected by genotype 1 but not genotype 6 HCV, with 80% versus 38% of patients infected by genotype 1 achieved SVR carried TT versus TG/GG at rs8099917, respectively (P=0.003). PNPLA3 genotype was not associated with SVR.

CONCLUSIONS

IL28B gene with rs8099917 T allele as an independent predictor of SVR in Chinese CHC patients infected by genotype 1 but not genotype 6 HCV.

Intensive temporal mapping of hepatitis C hypervariable region 1 quasispecies provides novel insights into hepatitis C virus evolution in chronic infection

Hepatitis C virus (HCV) is an RNA virus which exists as swarms of closely related viruses known as quasispecies (QS). A number of studies have demonstrated associations between QS hypervariable region 1 (HVR1) characteristics (diversity and complexity) and treatment success. We investigated HCV QS change in chronic infection over intervals of 2-4 weeks in 23 chronically infected individuals to describe the natural history of virus evolution and establish whether HCV QS characteristics could be used to individualize treatment regimens at a molecular level. HVR1 QS diversity, complexity and divergence continue to change in an unpredictable fashion in chronic infection even where there is little phylogenetic change, which is likely to preclude the use of these features in treatment individualization. Our phylogenetic analysis identified no change in the HVR1 QS in 12 subjects, minor change in four subjects and we describe a time-ordered phylogeny for the first time over a period as short as 16 weeks in seven subjects. We identified the existence of multiple subpopulation infections using partitioned analysis of QS and illustrated how subpopulations were sequentially replaced in a number of subjects. We illustrated marked variation in the nucleotide substitution per codon position between patients with sequence change and those without change in the phylogenetic tree. Analysis of codon-specific selection pressures identified a number of codons under purifying selection, suggesting that these code for structurally conserved amino acids. We also identified sections of the HVR1 under positive selection with marked sequence heterogeneity, suggesting that these may be potential epitope-binding sites.

Influence of host resistance on viral adaptation: hepatitis C virus as a case study

Genetic and cellular studies have shown that the host’s innate and adaptive immune responses are an important correlate of viral infection outcome. The features of the host’s immune response (host resistance) reflect the coevolution between hosts and pathogens that has occurred over millennia, and that has also resulted in a number of strategies developed by viruses to improve fitness and survival within the host (viral adaptation). In this review, we discuss viral adaptation to host immune pressure via protein–protein interactions and sequence-specific mutations. Specifically, we will present the “state of play” on viral escape mutations to host T-cell responses in the context of the hepatitis C virus, and their influence on infection outcome.

IFNL3 (IL28B) and IFNL4 polymorphisms are associated with treatment response in Thai patients infected with HCV genotype 1, but not with genotypes 3 and 6

Recent studies have shown an association between single nucleotide polymorphisms (SNPs) in the interferon lambda-3 (IFNL3 or IL-28B) and IFNL4 genes and treatment response to hepatitis C virus genotype 1 (HCV-1) infection. The importance of these SNPs for HCV genotype 3 (HCV-3), and particularly HCV genotype 6 (HCV-6), remains to be elucidated. We analyzed a cohort of 225 Thai individuals with chronic HCV infection treated with pegylated-interferon and ribavirin, of whom 69 (30.7%), 114 (50.7%) and 42 (18.6%) patients were infected with HCV-1, HCV-3, and HCV-6, respectively. DNA extracted from blood samples was analyzed for the SNPs rs12979860 and ss469415590. The distribution of CC, CT, and TT genotypes of rs12979860 was 189 (84%), 28 (12.4%) and 8 (3.6%), respectively, while the distribution of TT/TT, ΔG/TT, and ΔG/ΔG genotypes of ss469415590 was 192(85.3%), 28(12.5%), and 5(2.2%), respectively. Significantly lower frequencies of the favorable genotypes CC (for rs12979860) and TT/TT (for ss469415590) were found in the HCV-1 group in comparison with the other groups. The favorable genotypes were associated significantly with rapid and sustained virological response in the HCV-1 group. However, they were only associated with rapid virological response in the HCV-3 and HCV-6 groups. Furthermore, both SNPs were associated equally with the treatment outcome in the HCV-1 group. In contrast, the role of these SNPs in predicting treatment response was attenuated in the HCV-3 and HCV-6 groups. Thus, identification of these SNPs may be useful only in patients with refractory HCV-1 infection.

Comparison of quasispecies diversity of HCV between chronic hepatitis c and hepatocellular carcinoma by Ultradeep pyrosequencing

Backgrounds. Hepatitis C virus (HCV) exists as population of closely related genetic variants known as quasispecies. HCV quasispecies diversity is strongly influenced by host immune pressure on virus. Quasispecies diversity is expected to decline as host immune response to HCV decreases over natural course of progressing from chronic hepatitis C (CHC) to hepatocellular carcinoma (HCC). Methods. Ultradeep pyrosequencing (UDPS) was used to evaluate degree of quasispecies diversity in 49 patients infected with HCV including 26 with CHC and 23 with HCC. Whole structural protein of HCV genome was subjected to UDPS. Results. Shannon's indices for quasispecies diversity in HCV E1 were significantly lower in patients with HCC than in those with CHC. 14 amino acid positions differed significantly between two groups. Area under curve of ROC analysis for differentiating HCC from CHC was >0.8 for all of 14 amino acid positions. Conclusion. HCV quasispecies diversity as indicator of declining host immune functions was easily assessed by UDPS technology. Shannon's indices in 14 amino acid positions were found to differentiate between patients with CHC and those with HCC. Our data propose that degree of HCV quasispecies measured by UDPS might be useful to predict progression of HCC in chronic HCV patients.

In vitro inhibition of hepatitis C virus by antisense oligonucleotides in PBMC compared to hepatoma cells

AIM

To assess the efficiency of phosphorothioate antisense oligodeoxynucleotide 1 (S-ODN1) on HCV translation inhibition in PBMC compared to hepatoma cells in vitro for the first time.

MATERIALS AND METHODS

The study included 34 treatment naive HCV patients. IRES domain III and IV sequence variations were tested in 45 clones from 9 HCV patients. PBMC of HCV positive patients were subjected to S-ODN in vitro. Concomitantly HepG2 cells infected by the same patient's serum were also treated with S-ODN1 for 24 and 48 hours. Cellular RNA was tested for HCV plus and minus strands by reverse transcription polymerase chain reaction (RT-PCR).

RESULTS

Sequence variations were seen in HCV IRES domain III only while domain IV was conserved among all the tested patient's clones. S-ODN1 successfully inhibited HCV translation in HepG2 cells, while in PBMC inhibition was partial.

CONCLUSION

HCV IRES domain IV is more conserved than domain IIId in genotype 4 HCV patients. S-ODN against HCV IRES domain IV was not efficient to inhibit HCV translation in PBMC under the study conditions. Further studies testing other S-ODN targeting other HCV IRES domains in PBMC should be done.

Inference with viral quasispecies diversity indices: clonal and NGS approaches

Given the inherent dynamics of a viral quasispecies, we are often interested in the comparison of diversity indices of sequential samples of a patient, or in the comparison of diversity indices of virus in groups of patients in a treated versus control design. It is then important to make sure that the diversity measures from each sample may be compared with no bias and within a consistent statistical framework. In the present report, we review some indices often used as measures for viral quasispecies complexity and provide means for statistical inference, applying procedures taken from the ecology field. In particular, we examine the Shannon entropy and the mutation frequency, and we discuss the appropriateness of different normalization methods of the Shannon entropy found in the literature. By taking amplicons ultra-deep pyrosequencing (UDPS) raw data as a surrogate of a real hepatitis C virus viral population, we study through in-silico sampling the statistical properties of these indices under two methods of viral quasispecies sampling, classical cloning followed by Sanger sequencing (CCSS) and next-generation sequencing (NGS) such as UDPS. We propose solutions specific to each of the two sampling methods-CCSS and NGS-to guarantee statistically conforming conclusions as free of bias as possible.

CONTACT

josep.gregori@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Interaction of Hepatitis C Viral Proteins with Cellular Oncoproteins in the Induction of Liver Cancer

Hepatitis C virus infection is a major health problem all over the world. A large proportion of patients infected by HCV develop liver cirrhosis or cancer. However, the mechanism(s) remain to be elucidated. Since HCV does not carry any known oncogene, it is thought that interaction between virally encoded proteins and host proteins is responsible for carcinogenesis. Many crucial interactions between HCV-encoded proteins and host proteins have been reported. In this review we focus on the interaction of viral proteins with important regulators of cell cycle—oncoproteins YB-1, p53, and cyclin D1—which play a major role in cell proliferation, apoptosis, DNA repair, and genomic stability. Genetic variants of HCV accumulate in patients and alter these interactions of host cell proteins. It is a battle between the virus and host and the final outcome depends on the winner; if the host succeeds in clearing the virus the patient may not develop serious liver diseases. On the other hand, if the virus dominates by evolving quasispecies which code for altered proteins that interact differently with host proteins, or induce mutations in host protooncogenes, then the patient may develop liver cirrhosis and/or liver cancer.

[Molecular aspects of the antiviral response against hepatitis C virus implicated in vaccines development]

Hepatitis C is a contagious liver disease caused by hepacivirus of the Flaviviridae family. It has a RNA genome, a unique highly variable molecule. It encodes ten proteins which are necessary to infect cells and multiply. Replication occurs only in hepatocytes. Because of its wide genomic variability and the absence of symptoms, it is difficult to make an early diagnosis and successful treatment. In this review we analyze the molecular mechanism by which the virus infects the hepatocytes and causes the disease. We focused the analysis on different therapies, with the possibility of improving treatment with the use of new specific vaccines. We highlight the use of new therapies based on nucleic acids, mainly DNA vectors. In the near future, once this treatment is adequately evaluated in clinical trials, and the costs are calculated, it could be a very beneficial alternative to conventional methods.

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.

Genetic analysis of West Nile virus isolates from an outbreak in Idaho, United States, 2006-2007

West Nile virus (WNV) appeared in the U.S. in 1999 and has since become endemic, with yearly summer epidemics causing tens of thousands of cases of serious disease over the past 14 years. Analysis of WNV strains isolated during the 2006–2007 epidemic seasons demonstrates that a new genetic variant had emerged coincidentally with an intense outbreak in Idaho during 2006. The isolates belonging to the new variant carry a 13 nt deletion, termed ID-Δ13, located at the variable region of the 3′UTR, and are genetically related. The analysis of deletions and insertions in the 3′UTR of two major lineages of WNV revealed the presence of conserved repeats and two indel motifs in the variable region of the 3′UTR. One human and two bird isolates from the Idaho 2006–2007 outbreaks were sequenced using Illumina technology and within-host variability was analyzed. Continued monitoring of new genetic variants is important for public health as WNV continues to evolve.

Dynamics of defective hepatitis C virus clones in reinfected liver grafts in liver transplant recipients: ultradeep sequencing analysis

Hepatitis C virus (HCV) reinfects liver allografts in transplant recipients by replicating immediately after transplantation, causing a rapid increase in blood serum HCV RNA levels. We evaluated dynamic changes in the viral genetic complexity after HCV reinfection of the graft liver; we also identified the characteristics of replicating HCV clones using a massively parallel ultradeep sequencing technique to determine the full-genome HCV sequences in the liver and serum specimens of five transplant recipients with genotype 1b HCV infection before and after liver transplantation. The recipients showed extremely high genetic heterogeneity before transplantation, and the HCV population makeup was not significantly different between the liver and blood serum specimens of the individuals. Viral quasispecies complexity in serum was significantly lower after liver transplantation than before it, suggesting that certain HCV clones selectively proliferated after transplantation. Defective HCV clones lacking the structural region of the HCV genome did not increase in number, and full-genome HCV clones selectively increased in number immediately after liver transplantation. A re-increase in the same defective clone existing before transplantation was detected 22 months after transplantation in one patient. Ultradeep sequencing technology revealed that the genetic heterogeneity of HCV was reduced after liver transplantation. Dynamic changes in defective HCV clones after liver transplantation indicate that these clones have important roles in the HCV life cycle.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Metagenomic abundance estimation and diagnostic testing on species level

One goal of sequencing-based metagenomic community analysis is the quantitative taxonomic assessment of microbial community compositions. In particular, relative quantification of taxons is of high relevance for metagenomic diagnostics or microbial community comparison. However, the majority of existing approaches quantify at low resolution (e.g. at phylum level), rely on the existence of special genes (e.g. 16S), or have severe problems discerning species with highly similar genome sequences. Yet, problems as metagenomic diagnostics require accurate quantification on species level. We developed Genome Abundance Similarity Correction (GASiC), a method to estimate true genome abundances via read alignment by considering reference genome similarities in a non-negative LASSO approach. We demonstrate GASiC’s superior performance over existing methods on simulated benchmark data as well as on real data. In addition, we present applications to datasets of both bacterial DNA and viral RNA source. We further discuss our approach as an alternative to PCR-based DNA quantification.

Response-guided therapy for patients with hepatitis C virus genotype 6 infection: a pilot study

The optimal duration of treatment with pegylated interferon (PEG-IFN) plus ribavirin (RBV) in patients with hepatitis C virus (HCV) genotype 6 is unknown. This study was aimed at determining treatment response on the basis of rapid virological response (RVR) of HCV genotype 6 in comparison with genotypes 1 and 3. Sixty-six treatment naïve patients were treated with PEG-IFN-α2a (180 μg/week) plus weight-based RBV (1000-1200 mg/day). Patients with genotype 1 n = 16) and genotype 3 (n = 16) were treated for a fixed duration of 48 and 24 weeks, respectively. Patients with genotype 6 (n = 34) who achieved RVR were treated for 24 weeks (response-guided therapy) and the remaining patients were treated for 48 weeks (standard therapy). The mean baseline HCV RNA levels were not statistically different between groups (6.4 ± 0.8, 6.0 ± 1.0 and 6.5 ± 0.8 Log(10) IU/mL for genotypes 1, 3 and 6, respectively). Patients with genotypes 1, 3 and 6 achieved RVR in 43.8%, 87.5% and 73.5% of cases, respectively. One patient with genotype 1 and 3 with genotype 6 were considered nonresponders and discontinued therapy. Sustained virological response (SVR) was achieved in 62.5%, 81.3% and 76.5% of patients with genotypes 1, 3 and 6, respectively. The SVR rate in patients with genotype 6 who underwent response-guided therapy was 88%. This pilot study suggested that the SVR rate of HCV genotype 6 was at an intermediate level between those of genotypes 3 and 1. Treatment with PEG-IFN plus RBV for 24 weeks may be sufficient for patients with genotype 6 who achieve RVR. Prospective randomized trials are required to evaluate this response-guided strategy in a larger number of patients with genotype 6.

Dynamics of hepatitis B virus quasispecies in association with nucleos(t)ide analogue treatment determined by ultra-deep sequencing

BACKGROUND AND AIMS

Although the advent of ultra-deep sequencing technology allows for the analysis of heretofore-undetectable minor viral mutants, a limited amount of information is currently available regarding the clinical implications of hepatitis B virus (HBV) genomic heterogeneity.

METHODS

To characterize the HBV genetic heterogeneity in association with anti-viral therapy, we performed ultra-deep sequencing of full-genome HBV in the liver and serum of 19 patients with chronic viral infection, including 14 therapy-naïve and 5 nucleos(t)ide analogue(NA)-treated cases.

RESULTS

Most genomic changes observed in viral variants were single base substitutions and were widely distributed throughout the HBV genome. Four of eight (50%) chronic therapy-naïve HBeAg-negative patients showed a relatively low prevalence of the G1896A pre-core (pre-C) mutant in the liver tissues, suggesting that other mutations were involved in their HBeAg seroconversion. Interestingly, liver tissues in 4 of 5 (80%) of the chronic NA-treated anti-HBe-positive cases had extremely low levels of the G1896A pre-C mutant (0.0%, 0.0%, 0.1%, and 1.1%), suggesting the high sensitivity of the G1896A pre-C mutant to NA. Moreover, various abundances of clones resistant to NA were common in both the liver and serum of treatment-naïve patients, and the proportion of M204VI mutants resistant to lamivudine and entecavir expanded in response to entecavir treatment in the serum of 35.7% (5/14) of patients, suggesting the putative risk of developing drug resistance to NA.

CONCLUSION

Our findings illustrate the strong advantage of deep sequencing on viral genome as a tool for dissecting the pathophysiology of HBV infection.

Innate immunity and hepatitis C virus infection: a microarray's view

Hepatitis C virus (HCV) induces a chronic infection in more than two-thirds of HCV infected subjects. The inefficient innate and adaptive immune responses have been shown to play a major pathogenetic role in the development and persistence of HCV chronic infection. Several aspects of the interactions between the virus and the host immune system have been clarified and, in particular, mechanisms have been identified which underlie the ability of HCV to seize and subvert innate as well as adaptive immune responses. The present review summarizes recent findings on the interaction between HCV infection and innate immune response whose final effect is the downstream inefficient development of antigen-specific adaptive immunity, thereby contributing to virus persistence.

Biased mutational pattern and quasispecies hypothesis in H5N1 virus

Like other RNA viruses, influenza viruses are subject to high mutation rates. Carrying segmented RNA genomes, their genetic variability is even higher. We aimed at analyzing the mutational events occurring during the infection of chickens by the Highly Pathogenic Avian Influenza (HPAI) H5N1 virus. We therefore studied the different sequences of two surface proteins, hemagglutinin (HA) and neuraminidase (NA), as well as two internal proteins, PB2 and NS. Three organs (lung, spleen, brain) were obtained from a chicken, experimentally infected with a lethal dose of HPAI H5N1 virus. Cloning these PCR fragments enabled us to investigate the mutations undergone by the virus after several replicative cycles. The first outcome is the presence of a strong mutational bias, resembling host-driven ADAR1 adenosine deamination, which is responsible for 81% of all mutations. Whereas the frequency of RNA dependent RNA polymerase-related mutations is compatible with the survival of the virus, the ADAR1-like activity usually strongly increases the mutation frequency into a level of “error catastrophe” in theory incompatible with virus survival. Nevertheless, the virus was successfully infective. HPAI H5N1 virus displayed traits in agreement with the quasispecies theory. The role of this quasispecies structure in successful infection and the superposition with the ADAR1-like response is discussed.

Genetic heterogeneity of hepatitis C virus in association with antiviral therapy determined by ultra-deep sequencing

Background and Aims

The hepatitis C virus (HCV) invariably shows wide heterogeneity in infected patients, referred to as a quasispecies population. Massive amounts of genetic information due to the abundance of HCV variants could be an obstacle to evaluate the viral genetic heterogeneity in detail.

Methods

Using a newly developed massive-parallel ultra-deep sequencing technique, we investigated the viral genetic heterogeneity in 27 chronic hepatitis C patients receiving peg-interferon (IFN) α2b plus ribavirin therapy.

Results

Ultra-deep sequencing determined a total of more than 10 million nucleotides of the HCV genome, corresponding to a mean of more than 1000 clones in each specimen, and unveiled extremely high genetic heterogeneity in the genotype 1b HCV population. There was no significant difference in the level of viral complexity between immediate virologic responders and non-responders at baseline (p = 0.39). Immediate virologic responders (n = 8) showed a significant reduction in the genetic complexity spanning all the viral genetic regions at the early phase of IFN administration (p = 0.037). In contrast, non-virologic responders (n = 8) showed no significant changes in the level of viral quasispecies (p = 0.12), indicating that very few viral clones are sensitive to IFN treatment. We also demonstrated that clones resistant to direct-acting antivirals for HCV, such as viral protease and polymerase inhibitors, preexist with various abundances in all 27 treatment-naïve patients, suggesting the risk of the development of drug resistance against these agents.

Conclusion

Use of the ultra-deep sequencing technology revealed massive genetic heterogeneity of HCV, which has important implications regarding the treatment response and outcome of antiviral therapy.

Quasispecies of dengue virus

Pathogenic viruses have RNA genomes that cause acute and chronic infections. These viruses replicate with high mutation rates and exhibit significant genetic diversity, so-called viral quasispecies. Viral quasispecies play an important role in chronic infectious diseases, but little is known about their involvement in acute infectious diseases such as dengue virus (DENV) infection. DENV, the most important human arbovirus, is a causative agent of dengue fever (DF) and dengue hemorrhagic fever (DHF). Accumulating observations suggest that DENV exists as an extremely diverse virus population, but its biological significance is unclear. In other virus diseases, quasispecies affect the therapeutic strategies using drugs and vaccines. Here, I describe the quasispecies of DENV and discuss the possible role of quasispecies in the pathogenesis of and therapeutic strategy against DENV infection in comparison with other viruses such as Hepatitis C virus, human immunodeficiency virus type 1, and poliovirus.

Inferring viral quasispecies spectra from 454 pyrosequencing reads

BACKGROUND

RNA viruses infecting a host usually exist as a set of closely related sequences, referred to as quasispecies. The genomic diversity of viral quasispecies is a subject of great interest, particularly for chronic infections, since it can lead to resistance to existing therapies. High-throughput sequencing is a promising approach to characterizing viral diversity, but unfortunately standard assembly software was originally designed for single genome assembly and cannot be used to simultaneously assemble and estimate the abundance of multiple closely related quasispecies sequences.

RESULTS

In this paper, we introduce a new Viral Spectrum Assembler (ViSpA) method for quasispecies spectrum reconstruction and compare it with the state-of-the-art ShoRAH tool on both simulated and real 454 pyrosequencing shotgun reads from HCV and HIV quasispecies. Experimental results show that ViSpA outperforms ShoRAH on simulated error-free reads, correctly assembling 10 out of 10 quasispecies and 29 sequences out of 40 quasispecies. While ShoRAH has a significant advantage over ViSpA on reads simulated with sequencing errors due to its advanced error correction algorithm, ViSpA is better at assembling the simulated reads after they have been corrected by ShoRAH. ViSpA also outperforms ShoRAH on real 454 reads. Indeed, 7 most frequent sequences reconstructed by ViSpA from a real HCV dataset are viable (do not contain internal stop codons), and the most frequent sequence was within 1% of the actual open reading frame obtained by cloning and Sanger sequencing. In contrast, only one of the sequences reconstructed by ShoRAH is viable. On a real HIV dataset, ShoRAH correctly inferred only 2 quasispecies sequences with at most 4 mismatches whereas ViSpA correctly reconstructed 5 quasispecies with at most 2 mismatches, and 2 out of 5 sequences were inferred without any mismatches. ViSpA source code is available at http://alla.cs.gsu.edu/~software/VISPA/vispa.html.

CONCLUSIONS

ViSpA enables accurate viral quasispecies spectrum reconstruction from 454 pyrosequencing reads. We are currently exploring extensions applicable to the analysis of high-throughput sequencing data from bacterial metagenomic samples and ecological samples of eukaryote populations.

Chronic hepatitis C

The goal of antiviral therapy for patients with chronic hepatitis C virus (HCV) infection is to attain a sustained virologic response (SVR), which is defined as undetectable serum HCV-RNA levels at 6 months after the cessation of treatment. Major improvements in antiviral therapy for chronic hepatitis C have occurred in the past decade. The addition of ribavirin to interferon-alfa therapy and the introduction of pegylated interferon (PEG-IFN) have substantially improved SVR rates in patients with chronic hepatitis C. The optimization of HCV therapy with PEG-IFN and ribavirin continues to evolve. Studies are ongoing that use viral kinetics to tailor therapy to an individual's antiviral response and determine the ideal length of treatment to maximize the chance of SVR. Improved SVR can be achieved with new specific inhibitors that target the HCV NS3/4A protease and the NS5B polymerase. Several long-term follow-up studies have shown that SVR, when achieved, is associated with a very low risk of virologic relapse. Furthermore, antiviral therapy can reduce the morbidity and mortality rates associated with chronic hepatitis C by reducing fibrosis progression, the incidence of cirrhosis, and hepatocellular carcinoma.

Beyond the consensus: dissecting within-host viral population diversity of foot-and-mouth disease virus by using next-generation genome sequencing

The diverse sequences of viral populations within individual hosts are the starting material for selection and subsequent evolution of RNA viruses such as foot-and-mouth disease virus (FMDV). Using next-generation sequencing (NGS) performed on a Genome Analyzer platform (Illumina), this study compared the viral populations within two bovine epithelial samples (foot lesions) from a single animal with the inoculum used to initiate experimental infection. Genomic sequences were determined in duplicate sequencing runs, and the consensus sequence of the inoculum determined by NGS was identical to that previously determined using the Sanger method. However, NGS revealed the fine polymorphic substructure of the viral population, from nucleotide variants present at just below 50% frequency to those present at fractions of 1%. Some of the higher-frequency polymorphisms identified encoded changes within codons associated with heparan sulfate binding and were present in both foot lesions, revealing intermediate stages in the evolution of a tissue culture-adapted virus replicating within a mammalian host. We identified 2,622, 1,434, and 1,703 polymorphisms in the inoculum and in the two foot lesions, respectively: most of the substitutions occurred in only a small fraction of the population and represented the progeny from recent cellular replication prior to onset of any selective pressures. We estimated the upper limit for the genome-wide mutation rate of the virus within a cell to be 7.8 × 10(-4) per nucleotide. The greater depth of detection achieved by NGS demonstrates that this method is a powerful and valuable tool for the dissection of FMDV populations within hosts.

Mutant spectra in virus behavior

RNA viruses replicate as complex mutant spectra, also termed ‘mutant clouds’, known as viral quasispecies. While this is a widely observed viral population structure, it is less known that a number of biologically relevant features of this important group of viral pathogens depend on (or are strongly influenced by) the complexity and composition of mutant spectra. Among them, fitness increase or decrease depending on intrapopulation complementation or interference, selection triggered by memory genomes, pathogenic potential of viruses, disease evolution and the response to antiviral treatments. Quasispecies represent the recognition of complex behavior in viruses, and it is an oversimplification to equate such a population structure with the classic polymorphism of population biology. Darwinian principles acting on genome collectivities that replicate with high error rates provide a unique population structure prone to flexible and largely unpredictable behavior.

Cross-genotypic polyclonal anti-HCV antibodies from human ascitic fluid

Many anti-HCV antibodies are available, but more are needed for research and clinical applications. This study examines whether ascitic fluid from cirrhotic patients could be a source of reagent-grade antibodies. Ascitic fluid from 29 HCV patients was screened by ELISA for anti-HCV antibodies against three viral proteins: core, NS4B, and NS5A. Significant patient-to-patient variability in anti-HCV antibody titers was observed. Total ascitic fluid IgG purified by Protein-A chromatography reacted with HCV proteins in immunoblots, cell extracts, and replicon-expressing cells. Affinity-purification using synthetic peptides as bait allowed the preparation of cross-genotypic antibodies directed against pre-selected regions of HCV core, NS4B, and NS5A proteins. The performance of the polyclonal antibodies was comparable to that of monoclonal antibodies. Anti-NS4B antibody preparations reacted with genotype 1a, 1b, and 2a NS4B proteins in immunoblots and allowed NS4B to be localized in replicon-expressing cells. Ascitic fluid is an abundant source of human polyclonal cross-genotypic antibodies that can be used as an alternative to blood. This study shows the utility of selectively purifying human polyclonal antibodies from ascitic fluid. Affinity purification allows antibodies to be selected that are comparable to monoclonal antibodies in their ability to react with targeted regions of viral proteins.

Thermostability of subpopulations of H2N3 influenza virus isolates from mallard ducks

Maintenance of avian influenza virus in waterfowl populations requires that virions remain infectious while in the environment. Temperature has been shown to negatively correlate with persistence time, which is the duration for which virions are infectious. However, thermostability can vary between isolates regardless of subtype, and it is not known whether this variation occurs when host and geographic location of isolation are controlled. In this study, we analyzed the thermostabilities of 7 H2N3 viruses isolated from mallard ducks in Alberta, Canada. Virus samples were incubated at 37 degrees C and 55 degrees C, and infectivity titers were calculated at different time points. Based on the rate of infectivity inactivation at 37 degrees C, isolates could be grouped into either a thermosensitive or thermostable fraction for both egg- and MDCK-grown virus populations. Titers decreased more rapidly for isolates incubated at 55 degrees C, and this loss of infectivity occurred in a nonlinear, 2-step process, which is in contrast with the consensus on thermostability. This suggests that stock samples contain a mixture of subpopulations with different thermostabilities. The rate of decrease for the sensitive fraction was approximately 14 times higher than that for the stable fraction. The presence of subpopulations is further supported by selection experiments and plaque purification, both of which result in homogenous populations that exhibit linear decreases of infectivity titer. Therefore, variation of thermostability of influenza virus isolates begins at the level of the population. The presence of subpopulations with high thermostability suggests that avian viruses can persist in water longer than previously estimated, thus increasing the probability of transmission to susceptible hosts.