Analysis of high-depth sequence data for studying viral diversity: a comparison of next generation sequencing platforms using Segminator II

Background: Next generation sequencing provides detailed insight into the variation present within viral populations, introducing the possibility of treatment strategies that are both reactive and predictive. Current software tools, however, need to be scaled up to accommodate for high-depth viral data sets, which are often temporally or spatially linked. In addition, due to the development of novel sequencing platforms and chemistries, each with implicit strengths and weaknesses, it will be helpful for researchers to be able to routinely compare and combine data sets from different platforms/chemistries. In particular, error associated with a specific sequencing process must be quantified so that true biological variation may be identified. Results: Segminator II was developed to allow for the efficient comparison of data sets derived from different sources. We demonstrate its usage by comparing large data sets from 12 influenza H1N1 samples sequenced on both the 454 Life Sciences and Illumina platforms, permitting quantification of platform error. For mismatches median error rates at 0.10 and 0.12%, respectively, suggested that both platforms performed similarly. For insertions and deletions median error rates within the 454 data (at 0.3 and 0.2%, respectively) were significantly higher than those within the Illumina data (0.004 and 0.006%, respectively). In agreement with previous observations these higher rates were strongly associated with homopolymeric stretches on the 454 platform. Outside of such regions both platforms had similar indel error profiles. Additionally, we apply our software to the identification of low frequency variants. Conclusion: We have demonstrated, using Segminator II, that it is possible to distinguish platform specific error from biological variation using data derived from two different platforms. We have used this approach to quantify the amount of error present within the 454 and Illumina platforms in relation to genomic location as well as location on the read. Given that next generation data is increasingly important in the analysis of drug-resistance and vaccine trials, this software will be useful to the pathogen research community. A zip file containing the source code and jar file is freely available for download from http://www.bioinf.manchester.ac.uk/segminator/

The Evolutionary Analysis of Emerging Low Frequency HIV-1 CXCR4 Using Variants through Time—An Ultra-Deep Approach

Large-scale parallel pyrosequencing produces unprecedented quantities of sequence data. However, when generated from viral populations current mapping software is inadequate for dealing with the high levels of variation present, resulting in the potential for biased data loss. In order to apply the 454 Life Sciences' pyrosequencing system to the study of viral populations, we have developed software for the processing of highly variable sequence data. Here we demonstrate our software by analyzing two temporally sampled HIV-1 intra-patient datasets from a clinical study of maraviroc. This drug binds the CCR5 coreceptor, thus preventing HIV-1 infection of the cell. The objective is to determine viral tropism (CCR5 versus CXCR4 usage) and track the evolution of minority CXCR4-using variants that may limit the response to a maraviroc-containing treatment regimen. Five time points (two prior to treatment) were available from each patient. We first quantify the effects of divergence on initial read k-mer mapping and demonstrate the importance of utilizing population-specific template sequences in relation to the analysis of next-generation sequence data. Then, in conjunction with coreceptor prediction algorithms that infer HIV tropism, our software was used to quantify the viral population structure pre- and post-treatment. In both cases, low frequency CXCR4-using variants (2.5–15%) were detected prior to treatment. Following phylogenetic inference, these variants were observed to exist as distinct lineages that were maintained through time. Our analysis, thus confirms the role of pre-existing CXCR4-using virus in the emergence of maraviroc-insensitive HIV. The software will have utility for the study of intra-host viral diversity and evolution of other fast evolving viruses, and is available from http://www.bioinf.manchester.ac.uk/segminator/

From cassava to gari: Mapping of quality characteristics and end-user preferences in Cameroon and Nigeria

User's preferences of cassava and cassava products along the value chain are supported by specific root quality characteristics that can be linked to root traits. Therefore, providing an evidence base of user preferred characteristics along the value chain, can help in the functional choice of cassava varieties. In this respect, the present paper presents the results from focus group discussions and individual interviews on user preferred quality characteristics of raw cassava roots and the derived product, gari, ‐ one of the major cassava products in Sub Saharan Africa ‐ in major production and consumption areas of Cameroon and Nigeria. Choice of cassava varieties for farming is mainly determined by the multiple end‐uses of the roots, their agricultural yield and the processing determinants of roots that support their major high‐quality characteristics: size, density, low water content, maturity, colour and safety. Processing of cassava roots into gari goes through different technological variants leading to a gari whose high‐quality characteristics are: dryness, colour, shiny/attractive appearance, uniform granules and taste. Eba, the major consumption form of gari in Cameroon and Nigeria is mainly characterized by its textural properties: smoothness, firmness, stickiness, elasticity, mouldability. Recommendations are made, suggesting that breeding will have to start evaluating cassava clones for brightness/shininess, as well as textural properties such as mouldability and elasticity of cassava food products, for the purpose of supporting decision‐making by breeders and the development of high‐throughput selection methods of cassava varieties. Women are identified as important beneficiaries of such initiatives giving their disadvantaged position and their prominent role in cassava processing and marketing of gari

Track A Basic Science

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/138319/1/jia218438.pd

The growth of sociological theory : human nature, knowledge, and social change

xv, 527 p. : il.; 24 cm

Characterizing the diverse mutational pathways associated with r5-tropic maraviroc resistance: hiv-1 that uses the drug-bound ccr5 coreceptor

Entry inhibitors represent a potent class of antiretroviral drugs that target a host cell protein, CCR5, an HIV-1 entry coreceptor, and not viral protein. Lack of sensitivity can occur due to preexisting virus that uses the CXCR4 coreceptor, while true resistance occurs through viral adaptation to use a drug-bound CCR5 coreceptor. To understand this R5 resistance pathway, we analyzed >500 envelope protein sequences and phenotypes from viruses of 20 patients from the clinical trials MOTIVATE 1 and 2, in which treatment-experienced patients received maraviroc plus optimized background therapy. The resistant viral population was phylogenetically distinct and associated with a genetic bottleneck in each patient, consistent with de novo emergence of resistance. Recombination analysis showed that the C2-V3-C3 region tends to genotypically correspond to the recombinant\u27s phenotype, indicating its primary importance in conferring resistance. Between patients, there was a notable lack of commonality in the specific sites conferring resistance, confirming the unusual nature of R5-tropic resistance. We used coevolutionary and positive-selection analyses to characterize the genotypic determinants of resistance and found that (i) there are complicated covariation networks, indicating frequent coevolutionary/compensatory changes in the context of protein structure; (ii) covarying sites under positive selection are enriched in resistant viruses; (iii) CD4 binding sites form part of a unique covariation network independent of the V3 loop; and (iv) the covariation network formed between the V3 loop and other regions of gp120 and gp41 intersects sites involved in glycosylation and protein secretion. These results demonstrate that while envelope sequence mutations are the key to conferring maraviroc resistance, the specific changes involved are context dependent and thus inherently unpredictable. IMPORTANCE The entry inhibitor drug maraviroc makes the cell coreceptor CCR5 unavailable for use by HIV-1 and is now used in combination antiretroviral therapy. Treatment failure with drug-resistant virus is particularly interesting because it tends to be rare, with lack of sensitivity usually associated with the presence of CXCR4-using virus (CXCR4 is the main alternative coreceptor HIV-1 uses, in addition to CD4). We analyzed envelope sequences from HIV-1, obtained from 20 patients who enrolled in maraviroc clinical trials and experienced treatment failure, without detection of CXCR4-using virus. Evolutionary analysis was employed to identify molecular changes that confer maraviroc resistance. We found that in these individuals, resistant viruses form a distinct population that evolved once and was successful as a result of drug pressure. Further evolutionary analysis placed the complex network of interdependent mutational changes into functional groups that help explain the impediments to the emergence of maraviroc-associated R5 drug resistance

Characterizing the diverse mutational pathways associated with r5-tropic maraviroc resistance: hiv-1 that uses the drug-bound ccr5 coreceptor

Entry inhibitors represent a potent class of antiretroviral drugs that target a host cell protein, CCR5, an HIV-1 entry coreceptor, and not viral protein. Lack of sensitivity can occur due to preexisting virus that uses the CXCR4 coreceptor, while true resistance occurs through viral adaptation to use a drug-bound CCR5 coreceptor. To understand this R5 resistance pathway, we analyzed >500 envelope protein sequences and phenotypes from viruses of 20 patients from the clinical trials MOTIVATE 1 and 2, in which treatment-experienced patients received maraviroc plus optimized background therapy. The resistant viral population was phylogenetically distinct and associated with a genetic bottleneck in each patient, consistent with de novo emergence of resistance. Recombination analysis showed that the C2-V3-C3 region tends to genotypically correspond to the recombinant\u27s phenotype, indicating its primary importance in conferring resistance. Between patients, there was a notable lack of commonality in the specific sites conferring resistance, confirming the unusual nature of R5-tropic resistance. We used coevolutionary and positive-selection analyses to characterize the genotypic determinants of resistance and found that (i) there are complicated covariation networks, indicating frequent coevolutionary/compensatory changes in the context of protein structure; (ii) covarying sites under positive selection are enriched in resistant viruses; (iii) CD4 binding sites form part of a unique covariation network independent of the V3 loop; and (iv) the covariation network formed between the V3 loop and other regions of gp120 and gp41 intersects sites involved in glycosylation and protein secretion. These results demonstrate that while envelope sequence mutations are the key to conferring maraviroc resistance, the specific changes involved are context dependent and thus inherently unpredictable. IMPORTANCE The entry inhibitor drug maraviroc makes the cell coreceptor CCR5 unavailable for use by HIV-1 and is now used in combination antiretroviral therapy. Treatment failure with drug-resistant virus is particularly interesting because it tends to be rare, with lack of sensitivity usually associated with the presence of CXCR4-using virus (CXCR4 is the main alternative coreceptor HIV-1 uses, in addition to CD4). We analyzed envelope sequences from HIV-1, obtained from 20 patients who enrolled in maraviroc clinical trials and experienced treatment failure, without detection of CXCR4-using virus. Evolutionary analysis was employed to identify molecular changes that confer maraviroc resistance. We found that in these individuals, resistant viruses form a distinct population that evolved once and was successful as a result of drug pressure. Further evolutionary analysis placed the complex network of interdependent mutational changes into functional groups that help explain the impediments to the emergence of maraviroc-associated R5 drug resistance

Clonal analysis of HIV-1 genotype and function associated with virologic failure in treatment-experienced persons receiving maraviroc: Results from the MOTIVATE phase 3 randomized, placebo-controlled trials.

Detailed clonal phenotypic/genotypic analyses explored viral-escape mechanisms during maraviroc-based therapy in highly treatment-experienced participants from the MOTIVATE trials. To allow real-time assessment of samples while maintaining a blind trial, the first 267 enrolled participants were selected for evaluation. At failure, plasma samples from 20/50 participants (16/20 maraviroc-treated) with CXCR4-using virus and all 38 (13 maraviroc-treated) with CCR5-tropic virus were evaluated. Of those maraviroc-treated participants with CXCR4-using virus at failure, genotypic and phenotypic clonal tropism determinations showed >90% correspondence in 14/16 at Day 1 and 14/16 at failure. Phylogenetic analysis of clonal sequences detected pre-treatment progenitor CXCR4-using virus, or on-treatment virus highly divergent from the Day 1 R5 virus, excluding possible co-receptor switch through maraviroc-mediated evolution. Re-analysis of pre-treatment samples using the enhanced-sensitivity Trofile® assay detected CXCR4-using virus pre-treatment in 16/20 participants failing with CXCR4-using virus. Post-maraviroc reversion of CXCR4-use to CCR5-tropic occurred in 7/8 participants with follow-up, suggesting selective maraviroc inhibition of CCR5-tropic variants in a mixed-tropic viral population, not emergence of de novo mutations in CCR5-tropic virus, as the main virologic escape mechanism. Maraviroc-resistant CCR5-tropic virus was observed in plasma from 5 treated participants with virus displaying reduced maximal percent inhibition (MPI) but no evidence of IC50 change. Env clones with reduced MPI showed 1-5 amino acid changes specific to each V3-loop region of env relative to Day 1. However, transferring on-treatment resistance-associated changes using site-directed mutagenesis did not always establish resistance in Day 1 virus, and key 'signature' mutation patterns associated with reduced susceptibility to maraviroc were not identified. Evolutionary divergence of the CXCR4-using viruses is confirmed, emphasizing natural selection not influenced directly by maraviroc; maraviroc simply unmasks pre-existing lineages by inhibiting the R5 virus. For R5-viral failure, resistance development through drug selection pressure was uncommon and manifested through reduced MPI and with virus strain-specific mutational patterns