Viral Quasispecies Research Articles

Quantifying perinatal transmission of Hepatitis B viral quasispecies by tag linkage deep sequencing

Despite full immunoprophylaxis, mother-to-child transmission (MTCT) of Hepatitis B Virus still occurs in approximately 2–5% of HBsAg positive mothers. Little is known about the bottleneck of HBV transmission and the evolution of viral quasispecies in the context of MTCT. Here we adopted a newly developed tag linkage deep sequencing method and analyzed the quasispecies of four MTCT pairs that broke through immunoprophylaxis. By assigning unique tags to individual viral sequences, we accurately reconstructed HBV haplotypes in a region of 836 bp, which contains the major immune epitopes and drug resistance mutations. The detection limit of minor viral haplotypes reached 0.1% for individual patient sample. Dominance of “a determinant” polymorphisms were observed in two children, which pre-existed as minor quasispecies in maternal samples. In all four pairs of MTCT samples, we consistently observed a significant overlap of viral haplotypes shared between mother and child. We also demonstrate that the data can be potentially useful to estimate the bottleneck effect during HBV MTCT, which provides information to optimize treatment for reducing the frequency of MTCT.

Mother-to-Child HIV Transmission Bottleneck Selects for Consensus Virus with Lower Gag-Protease-Driven Replication Capacity.

In the large majority of cases, HIV infection is established by a single variant, and understanding the characteristics of successfully transmitted variants is relevant to prevention strategies. Few studies have investigated the viral determinants of mother-to-child transmission. To determine the impact of Gag-protease-driven viral replication capacity on mother-to-child transmission, the replication capacities of 148 recombinant viruses encoding plasma-derived Gag-protease from 53 nontransmitter mothers, 48 transmitter mothers, and 47 infected infants were assayed in an HIV-1-inducible green fluorescent protein reporter cell line. All study participants were infected with HIV-1 subtype C. There was no significant difference in replication capacities between the nontransmitter (n = 53) and transmitter (n = 44) mothers (P = 0.48). Infant-derived Gag-protease NL4-3 recombinant viruses (n = 41) were found to have a significantly lower Gag-protease-driven replication capacity than that of viruses derived from the mothers (P < 0.0001 by a paired t test). High percent similarities to consensus subtype C Gag, p17, p24, and protease sequences were also found in the infants (n = 28) in comparison to their mothers (P = 0.07, P = 0.002, P = 0.03, and P = 0.02, respectively, as determined by a paired t test). These data suggest that of the viral quasispecies found in mothers, the HIV mother-to-child transmission bottleneck favors the transmission of consensus-like viruses with lower viral replication capacities.IMPORTANCE Understanding the characteristics of successfully transmitted HIV variants has important implications for preventative interventions. Little is known about the viral determinants of HIV mother-to-child transmission (MTCT). We addressed the role of viral replication capacity driven by Gag, a major structural protein that is a significant determinant of overall viral replicative ability and an important target of the host immune response, in the MTCT bottleneck. This study advances our understanding of the genetic bottleneck in MTCT by revealing that viruses transmitted to infants have a lower replicative ability as well as a higher similarity to the population consensus (in this case HIV subtype C) than those of their mothers. Furthermore, the observation that "consensus-like" virus sequences correspond to lower in vitro replication abilities yet appear to be preferentially transmitted suggests that viral characteristics favoring transmission are decoupled from those that enhance replicative capacity.

Identification of HCV Resistant Variants against Direct Acting Antivirals in Plasma and Liver of Treatment Na\xefve Patients

Current standard-of-care treatment of chronically infected hepatitis C virus (HCV) patients involves direct-acting antivirals (DAA). However, concerns exist regarding the emergence of drug -resistant variants and subsequent treatment failure. In this study, we investigate potential natural drug-resistance mutations in the NS5B gene of HCV genotype 1b from treatment-naïve patients. Population-based sequencing and 454 deep sequencing of NS5B gene were performed on plasma and liver samples obtained from 18 treatment- naïve patients. The quasispecies distribution in plasma and liver samples showed a remarkable overlap in each patient. Although unique sequences in plasma or liver were observed, in the majority of cases the most dominant sequences were shown to be identical in both compartments. Neither in plasma nor in the liver codon changes were detected at position 282 that cause resistance to nucleos(t)ide analogues. However, in 10 patients the V321I change conferring resistance to nucleos(t)ide NS5B polymerase inhibitors and in 16 patients the C316N/Y/H non-nucleoside inhibitors were found mainly in liver samples. In conclusion, 454-deep sequencing of liver and plasma compartments in treatment naïve patients provides insight into viral quasispecies and the pre-existence of some drug-resistant variants in the liver, which are not necessarily present in plasma.

Catch Shiny Droplets in Suspension—Finding the Needle in a Haystack

In a recent issue of Cell Chemical Biology, Chaipan etal. (2017) described a high-throughput screening methodology to identify epitopes on HIV-1 particles recognized by broadly neutralizing antibodies. The approach utilizes a droplet-based microfluidics platform combining robust phenotypic single-virus sorting with next-generation sequencing of viral quasispecies.

QUENTIN: reconstruction of disease transmissions from viral quasispecies genomic data.

Genomic analysis has become one of the major tools for disease outbreak investigations. However, existing computational frameworks for inference of transmission history from viral genomic data often do not consider intra-host diversity of pathogens and heavily rely on additional epidemiological data, such as sampling times and exposure intervals. This impedes genomic analysis of outbreaks of highly mutable viruses associated with chronic infections, such as human immunodeficiency virus and hepatitis C virus, whose transmissions are often carried out through minor intra-host variants, while the additional epidemiological information often is either unavailable or has a limited use. The proposed framework QUasispecies Evolution, Network-based Transmission INference (QUENTIN) addresses the above challenges by evolutionary analysis of intra-host viral populations sampled by deep sequencing and Bayesian inference using general properties of social networks relevant to infection dissemination. This method allows inference of transmission direction even without the supporting case-specific epidemiological information, identify transmission clusters and reconstruct transmission history. QUENTIN was validated on experimental and simulated data, and applied to investigate HCV transmission within a community of hosts with high-risk behavior. It is available at https://github.com/skumsp/QUENTIN. pskums@gsu.edu or alexz@cs.gsu.edu or rahul@sfsu.edu or yek0@cdc.gov. Supplementary data are available at Bioinformatics online.

The pepper virome: natural co-infection of diverse viruses and their quasispecies

BackgroundThe co-infection of diverse viruses in a host plant is common; however, little is known about viral populations and their quasispecies in the host.ResultsHere, we report the first pepper viromes that were co-infected by different types of viral genomes. The pepper viromes are dominated by geminivirus DNA-A followed by a novel carlavirus referred to as Pepper virus A. The two pepper cultivars share similar viral populations and replications. However, the quasispecies for double-stranded RNA virus and two satellite DNAs were heterogeneous and homogenous in susceptible and resistant cultivars, respectively, indicating the quasispecies of an individual virus depends on the host.ConclusionsTaken together, we provide the first evidence that the host plant resistant to viruses has an unrevealed antiviral system, affecting viral quasispecies, not replication.

Quasispecies composition and evolution of a typical Zika virus clinical isolate from Suriname

The arthropod-borne Zika virus (ZIKV) is currently causing a major international public health threat in the Americas. This study describes the isolation of ZIKV from the plasma of a 29-year-old female traveler that developed typical symptoms, like rash, fever and headache upon return from Suriname. The complete genome sequence including the 5′ and 3′ untranslated regions was determined and phylogenetic analysis showed the isolate clustering within the Asian lineage, close to other viruses that have recently been isolated in the Americas. In addition, the viral quasispecies composition was analyzed by single molecule real time sequencing, which suggested a mutation frequency of 1.4 × 10−4 for this ZIKV isolate. Continued passaging of the virus in cell culture led to the selection of variants with mutations in NS1 and the E protein. The latter might influence virus binding to cell surface heparan sulfate.

Lessons learned from humoral responses of HIV patients.

Since 2009 many broadly neutralizing antibodies against HIV have been identified, yet there is still no vaccine capable of inducing such antibodies in humans. This review considers the early observations of HIV sera neutralization in light of more recent studies and highlights areas for future research. Large clinical cohort studies using standardized neutralization assays and pseudoviruses derived from primary isolates have shown that 10-30% of HIV infections result in some level of serum neutralization breadth. However, less than 10% of individuals develop a greater breadth of neutralization and are termed elite neutralizers. During HIV infection, many individuals develop strain-specific neutralization against their viral quasispecies, and similar immunogen-matched activity can now be induced in animal models. However, only in a minority of infections do broadly neutralizing antibodies develop. Therefore, understanding how the viral diversity, host immune environment, and antibody repertoires intersect to support the generation of neutralization breadth in elite neutralizers could provide guidelines as to how to improve immunization responses.

De novo assembly of highly polymorphic metagenomic data using in situ generated reference sequences and a novel BLAST-based assembly pipeline

BackgroundThe accuracy of metagenomic assembly is usually compromised by high levels of polymorphism due to divergent reads from the same genomic region recognized as different loci when sequenced and assembled together. A viral quasispecies is a group of abundant and diversified genetically related viruses found in a single carrier. Current mainstream assembly methods, such as Velvet and SOAPdenovo, were not originally intended for the assembly of such metagenomics data, and therefore demands for new methods to provide accurate and informative assembly results for metagenomic data.ResultsIn this study, we present a hybrid method for assembling highly polymorphic data combining the partial de novo-reference assembly (PDR) strategy and the BLAST-based assembly pipeline (BBAP). The PDR strategy generates in situ reference sequences through de novo assembly of a randomly extracted partial data set which is subsequently used for the reference assembly for the full data set. BBAP employs a greedy algorithm to assemble polymorphic reads. We used 12 hepatitis B virus quasispecies NGS data sets from a previous study to assess and compare the performance of both PDR and BBAP. Analyses suggest the high polymorphism of a full metagenomic data set leads to fragmentized de novo assembly results, whereas the biased or limited representation of external reference sequences included fewer reads into the assembly with lower assembly accuracy and variation sensitivity. In comparison, the PDR generated in situ reference sequence incorporated more reads into the final PDR assembly of the full metagenomics data set along with greater accuracy and higher variation sensitivity. BBAP assembly results also suggest higher assembly efficiency and accuracy compared to other assembly methods. Additionally, BBAP assembly recovered HBV structural variants that were not observed amongst assembly results of other methods. Together, PDR/BBAP assembly results were significantly better than other compared methods.ConclusionsBoth PDR and BBAP independently increased the assembly efficiency and accuracy of highly polymorphic data, and assembly performances were further improved when used together. BBAP also provides nucleotide frequency information. Together, PDR and BBAP provide powerful tools for metagenomic data studies.

De novo assembly of viral quasispecies using overlap graphs.

A viral quasispecies, the ensemble of viral strains populating an infected person, can be highly diverse. For optimal assessment of virulence, pathogenesis, and therapy selection, determining the haplotypes of the individual strains can play a key role. As many viruses are subject to high mutation and recombination rates, high-quality reference genomes are often not available at the time of a new disease outbreak. We present SAVAGE, a computational tool for reconstructing individual haplotypes of intra-host virus strains without the need for a high-quality reference genome. SAVAGE makes use of either FM-index–based data structures or ad hoc consensus reference sequence for constructing overlap graphs from patient sample data. In this overlap graph, nodes represent reads and/or contigs, while edges reflect that two reads/contigs, based on sound statistical considerations, represent identical haplotypic sequence. Following an iterative scheme, a new overlap assembly algorithm that is based on the enumeration of statistically well-calibrated groups of reads/contigs then efficiently reconstructs the individual haplotypes from this overlap graph. In benchmark experiments on simulated and on real deep-coverage data, SAVAGE drastically outperforms generic de novo assemblers as well as the only specialized de novo viral quasispecies assembler available so far. When run on ad hoc consensus reference sequence, SAVAGE performs very favorably in comparison with state-of-the-art reference genome-guided tools. We also apply SAVAGE on two deep-coverage samples of patients infected by the Zika and the hepatitis C virus, respectively, which sheds light on the genetic structures of the respective viral quasispecies.

Oldest Varroa tolerant honey bee population provides insight into the origins of the global decline of honey bees

The ecto-parasitic mite Varroa destructor has transformed the previously inconsequential Deformed Wing Virus (DWV) into the most important honey bee viral pathogen responsible for the death of millions of colonies worldwide. Naturally, DWV persists as a low level covert infection transmitted between nest-mates. It has long been speculated that Varroa via immunosuppression of the bees, activate a covert infection into an overt one. Here we show that despite Varroa feeding on a population of 20–40 colonies for over 30 years on the remote island of Fernando de Noronha, Brazil no such activation has occurred and DWV loads have remained at borderline levels of detection. This supports the alternative theory that for a new vector borne viral transmission cycle to start, an outbreak of an overt infection must first occur within the host. Therefore, we predict that this honey bee population is a ticking time-bomb, protected by its isolated position and small population size. This unique association between mite and bee persists due to the evolution of low Varroa reproduction rates. So the population is not adapted to tolerate Varroa and DWV, rather the viral quasispecies has simply not yet evolved the necessary mutations to produce a virulent variant.

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

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

Quasispecies evolution of the prototypical genotype 1 porcine reproductive and respiratory syndrome virus early during in vivo infection is rapid and tissue specific

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major infectious threat to the pig industry worldwide. Increasing evidence suggests that microevolution within a quasispecies population can give rise to high sequence heterogeneity in PRRSV; potentially impacting the pathogenicity of the virus. Here, we report on micro-evolutionary events taking place within the viral quasispecies population in lung and lymph node 3 days post infection (dpi) following experimental in vivo infection with the prototypical Lelystad PRRSV (LV). Sequence analysis revealed 16 high frequency single nucleotide variants (SNV) or differences from the reference LV genome which are assumed to be representative of the consensus inoculum genome. Additionally, 49 other low frequency SNVs were also found in the inoculum population. At 3 dpi, a total of 9 and 10 SNVs of varying frequencies could already be detected in the LV population infecting the lung and lymph nodes, respectively. Interestingly, of these, three and four novel SNVs emerged independently in the two respective tissues when compared to the inoculum. The remaining variants, though already present at lower frequencies in the inoculum, were positively selected and their frequency increased within the quasispecies population. Hence, we were able to determine directly from tissues infected with PRRSV the repertoire of genetic variants within the viral quasispecies population. Our data also suggest that microevolution of these variants is rapid and some may be tissue-specific.

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

The association between the evolution of hepatitis B virus (HBV) quasispecies and the development of hepatocellular carcinoma (HCC) is unknown. 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 9years 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. 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. 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.

Quasispecies of Hepatitis C Virus Participate in Cell-Specific Infectivity

It is well documented that a variety of viral quasispecies are found in the patients with chronic infection of hepatitis C virus (HCV). However, the significance of quasispecies in the specific infectivity to individual cell types remains unknown. In the present study, we analyzed the role of quasispecies of the genotype 2a clone, JFH1 (HCVcc), in specific infectivity to the hepatic cell lines, Huh7.5.1 and Hep3B. HCV RNA was electroporated into Huh7.5.1 cells and Hep3B/miR-122 cells expressing miR-122 at a high level. Then, we adapted the viruses to Huh7 and Hep3B/miR-122 cells by serial passages and termed the resulting viruses HCVcc/Huh7 and HCVcc/Hep3B, respectively. Interestingly, a higher viral load was obtained in the homologous combination of HCVcc/Huh7 in Huh7.5.1 cells or HCVcc/Hep3B in Hep3B/miR-122 cells compared with the heterologous combination. By using a reverse genetics system and deep sequence analysis, we identified several adaptive mutations involved in the high affinity for each cell line, suggesting that quasispecies of HCV participate in cell-specific infectivity.

Impaired human immunodeficiency virus type 1 replicative fitness in atypical viremic non-progressor individuals

BackgroundProgression rates from initial HIV-1 infection to advanced AIDS vary significantly among infected individuals. A distinct subgroup of HIV-1-infected individuals—termed viremic non-progressors (VNP) or controllers—do not seem to progress to AIDS, maintaining high CD4+ T cell counts despite high levels of viremia for many years. Several studies have evaluated multiple host factors, including immune activation, trying to elucidate the atypical HIV-1 disease progression in these patients; however, limited work has been done to characterize viral factors in viremic controllers.MethodsWe analyzed HIV-1 isolates from three VNP individuals and compared the replicative fitness, near full-length HIV-1 genomes and intra-patient HIV-1 genetic diversity with viruses from three typical (TP) and one rapid (RP) progressor individuals.ResultsViremic non-progressors and typical patients were infected for >10 years (range 10–17 years), with a mean CD4+ T-cell count of 472 cells/mm3 (442–529) and 400 cells/mm3 (126–789), respectively. VNP individuals had a less marked decline in CD4+ cells (mean −0.56, range −0.4 to −0.7 CD4+/month) than TP patients (mean −10.3, −8.2 to −13.1 CD4+/month). Interestingly, VNP individuals carried viruses with impaired replicative fitness, compared to HIV-1 isolates from the TP and RP patients (p < 0.05, 95% CI). Although analyses of the near full-length HIV-1 genomes showed no clear patterns of single-nucleotide polymorphisms (SNP) that could explain the decrease in replicative fitness, both the number of SNPs and HIV-1 population diversity correlated inversely with the replication capacity of the viruses (r = −0.956 and r = −0.878, p < 0.01, respectively).ConclusionIt is likely that complex multifactorial parameters govern HIV-1 disease progression in each individual, starting with the infecting virus (phenotype, load, and quasispecies diversity) and the intrinsic ability of the host to respond to the infection. Here we analyzed a subset of viremic controller patients and demonstrated that similar to the phenomenon observed in patients with a discordant response to antiretroviral therapy (i.e., high CD4+ cell counts with detectable plasma HIV-1 RNA load), reduced viral replicative fitness seems to be linked to slow disease progression in these antiretroviral-naïve individuals.

Deep sequencing and variant analysis of an Italian pathogenic field strain of equine infectious anaemia virus.

Equine infectious anaemia virus (EIAV) is a lentivirus with an almost worldwide distribution that causes persistent infections in equids. Technical limitations have restricted genetic analysis of EIAV field isolates predominantly to gag sequences resulting in very little published information concerning the extent of inter-strain variation in pol, env and the three ancillary open reading frames (ORFs). Here, we describe the use of long-range PCR in conjunction with next-generation sequencing (NGS) for rapid molecular characterization of all viral ORFs and known transcription factor binding motifs within the long terminal repeat of two EIAV isolates from the 2006 Italian outbreak. These isolates were from foals believed to have been exposed to the same source material but with different clinical histories: one died 53days post-infection (SA) while the other (DE) survived 5months despite experiencing multiple febrile episodes. Nucleotide sequence identity between the isolates was 99.358% confirming infection with the same EIAV strain with most differences comprising single nucleotide polymorphisms in env and the second exon of rev. Although the synonymous:non-synonymous nucleotide substitution ratio was approximately 2:1 in gag and pol, the situation is reversed in env and ORF3 suggesting these sequences are subjected to host-mediated selective pressure. EIAV proviral quasispecies complexity invivo has not been extensively investigated; however, analysis suggests it was relatively low in SA at the time of death. These results highlight advantages of NGS for molecular characterization of EIAV namely it avoids potential artefacts generated by traditional composite sequencing strategies and can provide information about viral quasispecies complexity.

Interferon-I: The Pièce de Résistance of HIV-1 Transmission?

Despite the extensive viral quasispecies that develops in an individual during the course of HIV-1 infection, transmission is typically established by a single donor viral variant. Recent studies now provide insight into the phenotypic properties influencing this selection process at transmission, including the contribution of resistance to type I interferons.

A comparative study on the characterization of hepatitis B virus quasispecies by clone-based sequencing and third-generation sequencing

Hepatitis B virus (HBV) has a high mutation rate due to the extremely high replication rate and the proofreading deficiency during reverse transcription. The generated variants with genetic heterogeneity are described as viral quasispecies (QS). Clone-based sequencing (CBS) is thought to be the ‘gold standard’ for assessing QS complexity and diversity of HBV, but an important issue about CBS is cost-effectiveness and laborious. In this study, we investigated the utility of the third-generation sequencing (TGS) DNA sequencing to characterize genetic heterogeneity of HBV QS and assessed the possible contribution of TGS technology in HBV QS studies. Parallel experiments including 3 control samples, which consisted of HBV full gene genotype B and genotype C plasmids, and 10 patients samples were performed by using CBS and TGS to analyze HBV whole-genome QS. Characterization of QS heterogeneity was conducted by using comprehensive statistical analysis. The results showed that TGS had a high consistency with CBS when measuring the complexity and diversity of QS. In addition, to detect rare variants, there were strong advantages conferred by TGS. In summary, TGS was considered to be practicable in HBV QS studies and it might have a relevant role in the clinical management of HBV infection in the future.Emerging Microbes & Infections (2017) 6, e100; doi:10.1038/emi.2017.88; published online 8 November 2017

Software Dedicated to Virus Sequence Analysis "Bioinformatics Goes Viral".

Computer-assisted technologies of the genomic structure, biological function, and evolution of viruses remain a largely neglected area of research. The attention of bioinformaticians to this challenging field is currently unsatisfying in respect to its medical and biological importance. The power of new genome sequencing technologies, associated with new tools to handle “big data”, provides unprecedented opportunities to address fundamental questions in virology. Here, we present an overview of the current technologies, challenges, and advantages of Next-Generation Sequencing (NGS) in relation to the field of virology. We present how viral sequences can be detected de novo out of current short-read NGS data. Furthermore, we discuss the challenges and applications of viral quasispecies and how secondary structures, commonly shaped by RNA viruses, can be computationally predicted. The phylogenetic analysis of viruses, as another ubiquitous field in virology, forms an essential element of describing viral epidemics and challenges current algorithms. Recently, the first specialized virus-bioinformatic organizations have been established. We need to bring together virologists and bioinformaticians and provide a platform for the implementation of interdisciplinary collaborative projects at local and international scales. Above all, there is an urgent need for dedicated software tools to tackle various challenges in virology.