Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.

Significance: With changes in land use and increased urbanization, the frequency with which pathogens jump species barriers to emerge in new hosts is expected to rise. Knowing which viruses may be more likely to become transmissible among humans, as opposed to only generating dead-end spillover infections, would be of considerable benefit to pandemic planning. Using multivariate modeling and multimodel inference, we sought to both identify and quantify those biological features of viruses that best determine interhuman transmissibility. This analysis revealed that chronic, nonsegmented, non–vector-borne, nonenveloped viruses with low host mortality had the highest likelihood of being transmissible among humans whereas genomic features had little predictive power. Our analysis therefore reveals that multiple virological features determine the likelihood of successful emergence.

Chikungunya virus (CHIKV) is an arthritogenic alphavirus causing epidemics of acute and chronic arthritic disease. Herein we describe a comprehensive RNA-Seq analysis of feet and lymph nodes at peak viraemia (day 2 post infection), acute arthritis (day 7) and chronic disease (day 30) in the CHIKV adult wild-type mouse model. Genes previously shown to be up-regulated in CHIKV patients were also up-regulated in the mouse model. CHIKV sequence information was also obtained with up to ≈8% of the reads mapping to the viral genome; however, no adaptive viral genome changes were apparent. Although day 2, 7 and 30 represent distinct stages of infection and disease, there was a pronounced overlap in up-regulated host genes and pathways. Type I interferon response genes (IRGs) represented up to ≈50% of up-regulated genes, even after loss of type I interferon induction on days 7 and 30. Bioinformatic analyses suggested a number of interferon response factors were primarily responsible for maintaining type I IRG induction. A group of genes prominent in the RNA-Seq analysis and hitherto unexplored in viral arthropathies were granzymes A, B and K. Granzyme A^-/- and to a lesser extent granzyme K^-/-, but not granzyme B^-/-, mice showed a pronounced reduction in foot swelling and arthritis, with analysis of granzyme A^-/- mice showing no reductions in viral loads but reduced NK and T cell infiltrates post CHIKV infection. Treatment with Serpinb6b, a granzyme A inhibitor, also reduced arthritic inflammation in wild-type mice. In non-human primates circulating granzyme A levels were elevated after CHIKV infection, with the increase correlating with viral load. Elevated granzyme A levels were also seen in a small cohort of human CHIKV patients. Taken together these results suggest granzyme A is an important driver of arthritic inflammation and a potential target for therapy. Trial Registration: ClinicalTrials.gov NCT00281294 The largest chikungunya virus (CHIKV) epidemic ever recorded began in 2004 in Africa and spread across Asia reaching Europe and recently the Americas, with millions of cases reported. We undertook a detailed analysis of the mRNA expression profile during acute and chronic arthritis in an adult wild-type mouse model of CHIKV infection and disease. Gene induction profiles showed a high concordance with published human data, providing some validation of the mouse model. The host response was overwhelmingly dominated by type I interferon response genes, even after type I interferon induction was lost. The analysis also provided information on CHIKV RNA, with no adaptive viral genome changes identified. An important goal of the analysis was to identify new players in arthritic inflammation. Granzyme A was prominent in the RNA-Seq data and granzyme A deficient mice showed reduced arthritis, with no effects on viral loads. Arthritic disease could also be ameliorated in wild-type mice with a granzyme A inhibitor. Elevated circulating granzyme A levels were seen in non-human primates infected with CHIKV and in human CHIKV patients. Granzyme A thus emerges to be a major driver of CHIKV-mediated arthritic inflammation and a potential target for anti-inflammatory interventions.

It is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART. This Europe-wide case–control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%–25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression. Two hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35–5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76–6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12–5.18, P = 0.024). A dose–effect relationship between virological failure and mutational load was found. Pre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

Viruses use the cellular machinery of their hosts for replication. It has therefore been proposed that the nucleotide and dinucleotide compositions of viruses should match those of their host species. If this is upheld, it may then be possible to use dinucleotide composition to predict the true host species of viruses sampled in metagenomic surveys. However, it is also clear that different taxonomic groups of viruses tend to have distinctive patterns of dinucleotide composition that may be independent of host species. To determine the relative strength of the effect of host versus virus family in shaping dinucleotide composition, we performed a comparative analysis of 20 RNA virus families from 15 host groupings, spanning two animal phyla and more than 900 virus species. In particular, we determined the odds ratios for the 16 possible dinucleotides and performed a discriminant analysis to evaluate the capability of virus dinucleotide composition to predict the correct virus family or host taxon from which it was isolated. Notably, while 81% of the data analyzed here were predicted to the correct virus family, only 62% of these data were predicted to their correct subphylum/class host and a mere 32% to their correct mammalian order. Similarly, dinucleotide composition has a weak predictive power for different hosts within individual virus families. We therefore conclude that dinucleotide composition is generally uniform within a virus family but less well reflects that of its host species. This has obvious implications for attempts to accurately predict host species from virus genome sequences alone. IMPORTANCE Determining the processes that shape virus genomes is central to understanding virus evolution and emergence. One question of particular importance is why nucleotide and dinucleotide frequencies differ so markedly between viruses. In particular, it is currently unclear whether host species or virus family has the biggest impact on dinucleotide frequencies and whether dinucleotide composition can be used to accurately predict host species. Using a comparative analysis, we show that dinucleotide composition has a strong phylogenetic association across different RNA virus families, such that dinucleotide composition can predict the family from which a virus sequence has been isolated. Conversely, dinucleotide composition has a poorer predictive power for the different host species within a virus family and across different virus families, indicating that the host has a relatively small impact on the dinucleotide composition of a virus genome.

Hepadnaviruses (hepatitis B viruses [HBVs]) are the only animal viruses that replicate their DNA by reverse transcription of an RNA intermediate. Until recently, the known host range of hepadnaviruses was limited to mammals and birds. We obtained and analyzed the first amphibian HBV genome, as well as several prototype fish HBVs, which allow the first comprehensive comparative genomic analysis of hepadnaviruses from four classes of vertebrates. Bluegill hepadnavirus (BGHBV) was characterized from in-house viral metagenomic sequencing. The African cichlid hepadnavirus (ACHBV) and the Tibetan frog hepadnavirus (TFHBV) were discovered using in silico analyses of the whole-genome shotgun and transcriptome shotgun assembly databases. Residues in the hydrophobic base of the capsid (core) proteins, designated motifs I, II, and III, are highly conserved, suggesting that structural constraints for proper capsid folding are key to capsid protein evolution. Surface proteins in all vertebrate HBVs contain similar predicted membrane topologies, characterized by three transmembrane domains. Most striking was the fact that BGHBV, ACHBV, and the previously described white sucker hepadnavirus did not form a fish-specific monophyletic group in the phylogenetic analysis of all three hepadnaviral genes. Notably, BGHBV was more closely related to the mammalian hepadnaviruses, indicating that cross-species transmission events have played a major role in viral evolution. Evidence of cross-species transmission was also observed with TFHBV. Hence, these data indicate that the evolutionary history of the hepadnaviruses is more complex than previously realized and combines both virus-host codivergence over millions of years and host species jumping. IMPORTANCE Hepadnaviruses are responsible for significant disease in humans (hepatitis B virus) and have been reported from a diverse range of vertebrates as both exogenous and endogenous viruses. We report the full-length genome of a novel hepadnavirus from a fish and the first hepadnavirus genome from an amphibian. The novel fish hepadnavirus, sampled from bluegills, was more closely related to mammalian hepadnaviruses than to other fish viruses. This phylogenetic pattern reveals that, although hepadnaviruses have likely been associated with vertebrates for hundreds of millions of years, they have also been characterized by species jumping across wide phylogenetic distances.