Dataset Information

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data.

ABSTRACT: Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (< 1%) within an infected host, and to assay variants without prior knowledge. Critical to interpreting deep sequence data sets is the ability to distinguish biological variants from process errors with high sensitivity and specificity. To address this challenge, we describe V-Phaser, an algorithm able to recognize rare biological variants in mixed populations. V-Phaser uses covariation (i.e. phasing) between observed variants to increase sensitivity and an expectation maximization algorithm that iteratively recalibrates base quality scores to increase specificity. Overall, V-Phaser achieved > 97% sensitivity and > 97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the -10 kb genome, including 603 rare variants (< 1% frequency) detected only using phase information. V-Phaser identified variants at frequencies down to 0.2%, comparable to the detection threshold of allele-specific PCR, a method that requires prior knowledge of the variants. The high sensitivity and specificity of V-Phaser enables identifying and tracking changes in low frequency variants in mixed populations such as RNA viruses.

SUBMITTER: Macalalad AR

PROVIDER: S-EPMC3305335 | biostudies-literature | 2012

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data.

Macalalad Alexander R AR Zody Michael C MC Charlebois Patrick P Lennon Niall J NJ Newman Ruchi M RM Malboeuf Christine M CM Ryan Elizabeth M EM Boutwell Christian L CL Power Karen A KA Brackney Doug E DE Pesko Kendra N KN Levin Joshua Z JZ Ebel Gregory D GD Allen Todd M TM Birren Bruce W BW Henn Matthew R MR

PLoS computational biology 20120315 3

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing met ...[more]

PMID: 22438797

Similar Datasets

Project description:BackgroundCharacterising genetic diversity through the analysis of massively parallel sequencing (MPS) data offers enormous potential to significantly improve our understanding of the genetic basis for observed phenotypes, including predisposition to and progression of complex human disease. Great challenges remain in resolving genetic variants that are genuine from the millions of artefactual signals.ResultsFAVR is a suite of new methods designed to work with commonly used MPS analysis pipelines to assist in the resolution of some of the issues related to the analysis of the vast amount of resulting data, with a focus on relatively rare genetic variants. To the best of our knowledge, no equivalent method has previously been described. The most important and novel aspect of FAVR is the use of signatures in comparator sequence alignment files during variant filtering, and annotation of variants potentially shared between individuals. The FAVR methods use these signatures to facilitate filtering of (i) platform and/or mapping-specific artefacts, (ii) common genetic variants, and, where relevant, (iii) artefacts derived from imbalanced paired-end sequencing, as well as annotation of genetic variants based on evidence of co-occurrence in individuals. We applied conventional variant calling applied to whole-exome sequencing datasets, produced using both SOLiD and TruSeq chemistries, with or without downstream processing by FAVR methods. We demonstrate a 3-fold smaller rare single nucleotide variant shortlist with no detected reduction in sensitivity. This analysis included Sanger sequencing of rare variant signals not evident in dbSNP131, assessment of known variant signal preservation, and comparison of observed and expected rare variant numbers across a range of first cousin pairs. The principles described herein were applied in our recent publication identifying XRCC2 as a new breast cancer risk gene and have been made publically available as a suite of software tools.ConclusionsFAVR is a platform-agnostic suite of methods that significantly enhances the analysis of large volumes of sequencing data for the study of rare genetic variants and their influence on phenotypes.

Project description:There have been two major eras in the history of gene discovery. The first was the era of linkage analysis, with approximately 1,300 disease-related genes identified by positional cloning by the turn of the millennium. The second era has been powered by two major breakthroughs: the publication of the human genome and the development of massively parallel sequencing (MPS). MPS has greatly accelerated disease gene identification, such that disease genes that would have taken years to map previously can now be determined in a matter of weeks. Additionally, the number of affected families needed to map a causative gene and the size of such families have fallen: de novo mutations, previously intractable by linkage analysis, can be identified through sequencing of the parent-child trio, and genes for recessive disease can be identified through MPS even of a single affected individual. MPS technologies include whole exome sequencing (WES), whole genome sequencing (WGS), and panel sequencing, each with their strengths. While WES has been responsible for most gene discoveries through MPS, WGS is superior in detecting copy number variants, chromosomal rearrangements, and repeat-rich regions. Panels are commonly used for diagnostic purposes as they are extremely cost-effective and generate manageable quantities of data, with no risk of unexpected findings. However, in instances of diagnostic uncertainty, it can be challenging to choose the right panel, and in these circumstances WES has a higher diagnostic yield. MPS has ethical, social, and legal implications, many of which are common to genetic testing generally but amplified due to the magnitude of data (e.g., relationship misattribution, identification of variants of uncertain significance, and genetic discrimination); others are unique to WES and WGS technologies (e.g., incidental or secondary findings). Nonetheless, MPS is rapidly translating into clinical practice as an extremely useful part of the clinical armamentarium.

Project description:BackgroundVirus-associated cell membrane proteins acquired by HIV-1 during budding may give information on the cellular source of circulating virions. In the present study, by applying immunosorting of the virus and of the cells with antibodies targeting monocyte (CD36) and lymphocyte (CD26) markers, it was possible to directly compare HIV-1 quasispecies archived in circulating monocytes and T lymphocytes with that present in plasma virions originated from the same cell types. Five chronically HIV-1 infected patients who underwent therapy interruption after prolonged HAART were enrolled in the study. The analysis was performed by the powerful technology of ultra-deep pyrosequencing after PCR amplification of part of the env gene, coding for the viral glycoprotein (gp) 120, encompassing the tropism-related V3 loop region. V3 amino acid sequences were used to establish heterogeneity parameters, to build phylogenetic trees and to predict co-receptor usage.ResultsThe heterogeneity of proviral and viral genomes derived from monocytes was higher than that of T-lymphocyte origin. Both monocytes and T lymphocytes might contribute to virus rebounding in the circulation after therapy interruptions, but other virus sources might also be involved. In addition, both proviral and circulating viral sequences from monocytes and T lymphocytes were predictive of a predominant R5 coreceptor usage. However, minor variants, segregating from the most frequent quasispecies variants, were present. In particular, in proviral genomes harboured by monocytes, minority variant clusters with a predicted X4 phenotype were found.ConclusionThis study provided the first direct comparison between the HIV-1 quasispecies archived as provirus in circulating monocytes and T lymphocytes with that of plasma virions replicating in the same cell types. Ultra-deep pyrosequencing generated data with some order of magnitude higher than any previously obtained with conventional approaches. Next generation sequencing allowed the analysis of previously inaccessible aspects of HIV-1 quasispecies, such as co-receptor usage of minority variants present in archived proviral sequences and in actually replicating virions, which may have clinical and therapeutic relevance.

Dataset Information

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data.

Publications

Highly sensitive and specific detection of rare variants in mixed viral populations from massively parallel sequence data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets