Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS

HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors - synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.

Single nucleotide polymorphism (SNP) frequencies and distribution reveal complex genetic composition of seven novel natural isolates of Cydia pomonella granulovirus

The co-evolution between baculoviruses and their insect hosts results in selection of virus populations. To explore this phenomenon at the molecular level, seven natural isolates of Cydia pomonella granulovirus (CpGV) collected from orchards in northwest China were studied using Illumina next generation sequencing (NGS). A total of 540 genome positions with single nucleotide polymorphisms (SNPs) were detected in comparison with known CpGV isolates. New members of previously defined phylogenetic genome groups A, D and E of CpGV, as well as two novel phylogenetic lines, termed genome group F and G, were identified. Combining SNP frequency distribution with the prevalence of genome group-specific SNPs, revealed that six isolates of CpGV were mixtures of different ratios of at least two genotypes, whereas only one isolate, CpGV-WW, was genetically highly homogeneous. This study significantly extends our current understanding of the genetic diversity of CpGV and opens new lines of application of this virus.

Method for serial passage of infectious hematopoietic necrosis virus (IHNV) in rainbow trout

Transmission is a fundamental component of pathogen fitness. A better understanding of pathogen transmission can greatly improve disease management. In particular, controlled studies of multiple rounds of natural transmission (i.e. serial passage) can provide powerful epidemiological and evolutionary inferences. However, such studies are possible in only a few systems because of the challenges in successfully initiating and maintaining transmission in the laboratory. Here we developed an efficient and reproducible cohabitation method for conducting controlled experiments investigating the effects of serial passage on infectious hematopoietic necrosis virus (IHNV) in rainbow trout. This method was used to investigate the transmission efficiency and kinetics of viral shedding of IHNV over 3 serial passages. Transmission efficiency decreased from 100 to 62.5% over the passage steps and was associated with a decrease in virus shedding into water. A shift in the peak of viral shedding was also observed, from Day 2 post immersion for passage 0 to at least 24 h later for all subsequent passages. Finally, the characterization of viruses after 1 round of transmission and propagation on cells showed no change in glycoprotein (G gene) sequences or viral virulence compared to the ancestral virus stock. The methods developed provide valuable tools for reproducible population-level studies of IHNV epidemiology and evolution.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

A novel framework for inferring parameters of transmission from viral sequence data

Transmission between hosts is a critical part of the viral lifecycle. Recent studies of viral transmission have used genome sequence data to evaluate the number of particles transmitted between hosts, and the role of selection as it operates during the transmission process. However, the interpretation of sequence data describing transmission events is a challenging task. We here present a novel and comprehensive framework for using short-read sequence data to understand viral transmission events, designed for influenza virus, but adaptable to other viral species. Our approach solves multiple shortcomings of previous methods for this purpose; for example, we consider transmission as an event involving whole viruses, rather than sets of independent alleles. We demonstrate how selection during transmission and noisy sequence data may each affect naive inferences of the population bottleneck, accounting for these in our framework so as to achieve a correct inference. We identify circumstances in which selection for increased viral transmission may or may not be identified from data. Applying our method to experimental data in which transmission occurs in the presence of strong selection, we show that our framework grants a more quantitative insight into transmission events than previous approaches, inferring the bottleneck in a manner that accounts for selection, both for within-host virulence, and for inherent viral transmissibility. Our work provides new opportunities for studying transmission processes in influenza, and by extension, in other infectious diseases.

Pyrosequencing method for sensitive detection of HBV drug resistance mutations

Hepatitis B (HBV) drug resistance assay is important for guiding therapy after the development of virologic breakthrough for patients receiving nucleoside/-tide analog therapy. However, the existing genotyping tools are either costly or lack sensitivity to detect mixed genotypes, and an improved method of resistant mutation detection is needed. An assay protocol for clinical application using pyrosequencing method was developed, capable of detecting all known validated HBV polymerase gene mutations that impart resistance to lamivudine, adefovir, tenofovir, and entecavir. Sixty-eight serum samples with known HBV resistance genotypes, previously tested with either Sanger sequencing assay or commercial line probe assay, were used for validation. Where there were discrepancies between the two methods, clonal sequencing by Sanger's method was used for confirmation. The modified pyrosequencing method accurately identified all the cloned polymerase genotypes and was able to distinguish as little as 5% of the mutant populations. This assay can be performed on serum sample with HBV DNA as low as 13.5 IU/mL. The cost per test was less than existing commercial assay. HBV drug resistance pyrosequencing assay was accurate, more sensitive and cheaper compared with the existing methods. It can detect minor populations of drug-resistant clones earlier, before the drug resistant clones become dominant, allowing the opportunity for an earlier change of therapy.

Comparison of lesional skin c-KIT mutations with clinical phenotype in patients with mastocytosis

BACKGROUND

Activating c-KIT mutations cause abnormal mast cell growth and appear to play a role in mastocytosis. However, the correlation of c-KIT mutations with disease phenotypes is poorly characterized.

AIM

To evaluate the correlation of c-KIT mutations with clinical presentations and laboratory findings.

METHODS

Total cellular RNA was isolated from the skin lesions of 43 adults and 7 children with mastocytosis, and PCR amplicons of cDNA were sequenced for c-KIT mutations.

RESULTS

The most common activating mutation, KIT-D816V, was identified in 72% of adults and 57% of children. Additional activating mutations, namely, V560G and the internal tandem duplications (ITDs) 502-503dupAY, were detected in 12% of adults and 8% of children. V560G occurred more commonly in our patients than previously reported, and it appeared to be associated with more advanced disease. Otherwise, the presence or absence of activating mutations did not correlate with skin lesion morphology, disease extent or total serum tryptase levels. Four adults had expression only of wild-type KIT, while two others had expression of a truncated KIT lacking tyrosine kinase activity; yet these patients were clinically indistinguishable from those patients with activating c-KIT mutations.

CONCLUSIONS

Activating c-KIT mutations exist in a significant portion of patients with mastocytosis, but not all patients showed expression of these mutations. Except for V560G, the presence or absence of activating c-KIT mutations did not predict the extent of disease. These observations suggest that although activating c-KIT mutations are associated with mast cell growth, other genes probably play a role in the cause of mastocytosis.

Sieve analysis of breakthrough HIV-1 sequences in HVTN 505 identifies vaccine pressure targeting the CD4 binding site of Env-gp120

Although the HVTN 505 DNA/recombinant adenovirus type 5 vector HIV-1 vaccine trial showed no overall efficacy, analysis of breakthrough HIV-1 sequences in participants can help determine whether vaccine-induced immune responses impacted viruses that caused infection. We analyzed 480 HIV-1 genomes sampled from 27 vaccine and 20 placebo recipients and found that intra-host HIV-1 diversity was significantly lower in vaccine recipients (P ≤ 0.04, Q-values ≤ 0.09) in Gag, Pol, Vif and envelope glycoprotein gp120 (Env-gp120). Furthermore, Env-gp120 sequences from vaccine recipients were significantly more distant from the subtype B vaccine insert than sequences from placebo recipients (P = 0.01, Q-value = 0.12). These vaccine effects were associated with signatures mapping to CD4 binding site and CD4-induced monoclonal antibody footprints. These results suggest either (i) no vaccine efficacy to block acquisition of any viral genotype but vaccine-accelerated Env evolution post-acquisition; or (ii) vaccine efficacy against HIV-1s with Env sequences closest to the vaccine insert combined with increased acquisition due to other factors, potentially including the vaccine vector.

On the effective depth of viral sequence data

Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

Single Genome Sequencing of Expressed and Proviral HIV-1 Envelope Glycoprotein 120 (gp120) and nef Genes

The current study provides detailed protocols utilized to amplify the complete HIV-1 gp120 and nef genes from single copies of expressed or integrated HIV present in fresh-frozen autopsy tissues of patients who died while on combined antiretroviral therapy (cART) with no detectable plasma viral load (pVL) at death ( Lamers et al., 2016a and 2016b; Rose et al., 2016 ). This method optimizes protocols from previous publications ( Palmer et al., 2005 ; Norström et al., 2012 ; Lamers et al., 2015 ; 2016a and 2016b; Rife et al., 2016 ) to produce single distinct PCR products that can be directly sequenced and includes several cost-saving and time-efficient modifications.

Error rates, PCR recombination, and sampling depth in HIV-1 whole genome deep sequencing

Deep sequencing is a powerful and cost-effective tool to characterize the genetic diversity and evolution of virus populations. While modern sequencing instruments readily cover viral genomes many thousand fold and very rare variants can in principle be detected, sequencing errors, amplification biases, and other artifacts can limit sensitivity and complicate data interpretation. For this reason, the number of studies using whole genome deep sequencing to characterize viral quasi-species in clinical samples is still limited. We have previously undertaken a large scale whole genome deep sequencing study of HIV-1 populations. Here we discuss the challenges, error profiles, control experiments, and computational test we developed to quantify the accuracy of variant frequency estimation.

Population genomics of intrapatient HIV-1 evolution

Many microbial populations rapidly adapt to changing environments with multiple variants competing for survival. To quantify such complex evolutionary dynamics in vivo, time resolved and genome wide data including rare variants are essential. We performed whole-genome deep sequencing of HIV-1 populations in 9 untreated patients, with 6-12 longitudinal samples per patient spanning 5-8 years of infection. The data can be accessed and explored via an interactive web application. We show that patterns of minor diversity are reproducible between patients and mirror global HIV-1 diversity, suggesting a universal landscape of fitness costs that control diversity. Reversions towards the ancestral HIV-1 sequence are observed throughout infection and account for almost one third of all sequence changes. Reversion rates depend strongly on conservation. Frequent recombination limits linkage disequilibrium to about 100 bp in most of the genome, but strong hitch-hiking due to short range linkage limits diversity.