Indel and Carryforward Correction (ICC): a new analysis approach for processing 454 pyrosequencing data.
Ontology highlight
ABSTRACT:
SUBMITTER: Deng W
PROVIDER: S-EPMC3777116 | biostudies-literature | 2013 Oct
REPOSITORIES: biostudies-literature
Publications
Indel and Carryforward Correction (ICC): a new analysis approach for processing 454 pyrosequencing data.
Deng Wenjie W Maust Brandon S BS Westfall Dylan H DH Chen Lennie L Zhao Hong H Larsen Brendan B BB Iyer Shyamala S Liu Yi Y Mullins James I JI
Bioinformatics (Oxford, England) 20130729 19
<h4>Motivation</h4>Pyrosequencing technology provides an important new approach to more extensively characterize diverse sequence populations and detect low frequency variants. However, the promise of this technology has been difficult to realize, as careful correction of sequencing errors is crucial to distinguish rare variants (∼1%) in an infected host with high sensitivity and specificity.<h4>Results</h4>We developed a new approach, referred to as Indel and Carryforward Correction (ICC), to c ...[more]
Similar Datasets
Project description:MOTIVATION: Throughout the recent years, 454 pyrosequencing has emerged as an efficient alternative to traditional Sanger sequencing and is widely used in both de novo whole-genome sequencing and metagenomics. Especially the latter application is extremely sensitive to sequencing errors and artificially duplicated reads. Both are common in 454 pyrosequencing and can create a strong bias in the estimation of diversity and composition of a sample. To date, there are several tools that aim to remove both sequencing noise and duplicates. Nevertheless, duplicate removal is often based on nucleotide sequences rather than on the underlying flow values, which contain additional information. RESULTS: With the novel tool JATAC, we present an approach towards a more accurate duplicate removal by analysing flow values directly. Making use of previous findings on 454 flow data characteristics, we combine read clustering with Bayesian distance measures. Finally, we provide a benchmark with an existing algorithm. AVAILABILITY: JATAC is freely available under the General Public License from http://malde.org/ketil/jatac/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
| S-EPMC3605598 | biostudies-literature
Project description:Pyrosequencing of the 16S ribosomal RNA gene (16S) has become one of the most popular methods to assess microbial diversity. Pyrosequencing reads containing ambiguous bases (Ns) are generally discarded based on the assumptions of their non-sequence-dependent formation and high error rates. However, taxonomic composition differed by removal of reads with Ns. We determined whether Ns from pyrosequencing occur in a sequence-dependent manner. Our reads and the corresponding flow value data revealed occurrence of sequence-specific N errors with a common sequential pattern (a homopolymer + a few nucleotides with bases other than the homopolymer + N) and revealed that the nucleotide base of the homopolymer is the true base for the following N. Using an algorithm reflecting this sequence-dependent pattern, we corrected the Ns in the 16S (86.54%), bphD (81.37%) and nifH (81.55%) amplicon reads from a mock community with high precisions of 95.4, 96.9 and 100%, respectively. The new N correction method was applicable for determining most of Ns in amplicon reads from a soil sample, resulting in reducing taxonomic biases associated with N errors and in shotgun sequencing reads from public metagenome data. The method improves the accuracy and precision of microbial community analysis and genome sequencing using 454 pyrosequencing.
| S-EPMC3985643 | biostudies-literature
Project description:BACKGROUND: Many potentially life-threatening infectious viruses are highly mutable in nature. Characterizing the fittest variants within a quasispecies from infected patients is expected to allow unprecedented opportunities to investigate the relationship between quasispecies diversity and disease epidemiology. The advent of next-generation sequencing technologies has allowed the study of virus diversity with high-throughput sequencing, although these methods come with higher rates of errors which can artificially increase diversity. RESULTS: Here we introduce a novel computational approach that incorporates base quality scores from next-generation sequencers for reconstructing viral genome sequences that simultaneously infers the number of variants within a quasispecies that are present. Comparisons on simulated and clinical data on dengue virus suggest that the novel approach provides a more accurate inference of the underlying number of variants within the quasispecies, which is vital for clinical efforts in mapping the within-host viral diversity. Sequence alignments generated by our approach are also found to exhibit lower rates of error. CONCLUSIONS: The ability to infer the viral quasispecies colony that is present within a human host provides the potential for a more accurate classification of the viral phenotype. Understanding the genomics of viruses will be relevant not just to studying how to control or even eradicate these viral infectious diseases, but also in learning about the innate protection in the human host against the viruses.
| S-EPMC4234478 | biostudies-literature
Project description:There has been a rapid proliferation of approaches for processing and manipulating second generation DNA sequence data. However, users are often left with uncertainties about how the choice of processing methods may impact biological interpretation of data. In this report, we probe differences in output between two different processing pipelines: a de-noising approach using the AmpliconNoise algorithm for error correction, and a standard approach using quality filtering and preclustering to reduce error. There was a large overlap in reads culled by each method, although AmpliconNoise removed a greater net number of reads. Most OTUs produced by one method had a clearly corresponding partner in the other. Although each method resulted in OTUs consisting entirely of reads that were culled by the other method, there were many more such OTUs formed in the standard pipeline. Total OTU richness was reduced by AmpliconNoise processing, but per-sample OTU richness, diversity and evenness were increased. Increases in per-sample richness and diversity may be a result of AmpliconNoise processing producing a more even OTU rank-abundance distribution. Because communities were randomly subsampled to equalize sample size across communities, and because rare sequence variants are less likely to be selected during subsampling, fewer OTUs were lost from individual communities when subsampling AmpliconNoise-processed data. In contrast to taxon-based diversity estimates, phylogenetic diversity was reduced even on a per-sample basis by de-noising, and samples switched widely in diversity rankings. This work illustrates the significant impacts of processing pipelines on the biological interpretations that can be made from pyrosequencing surveys. This study provides important cautions for analyses of contemporary data, for requisite data archiving (processed vs. non-processed data), and for drawing comparisons among studies performed using distinct data processing pipelines.
| S-EPMC3431371 | biostudies-literature
Project description:The feline gastrointestinal microbiota have direct influence on feline health and also human health as a reservoir for potential zoonotic pathogens and antibiotic resistant bacterial strains. In order to describe the feline gastrointestinal microbial diversity, fecal samples from cats have been characterized using both culture-dependent and culture-independent methods. However, data correlating total microbial composition and their functions are lacking. Present descriptive study evaluated both phylogenetic and metabolic diversity of the feline intestinal microbiota using GS Junior titanium shotgun pyrosequencing. A total of 152,494 pyrosequencing reads (5405 assembled contigs) were generated and classified into both phylogenetic and metabolic profiles of the feline intestinal microbiota. The Bacteroides/Chlorobi group was the most predominant bacterial phylum comprising ~68% of total classified diversity, followed by Firmicutes (~13%) and Proteobacteria (~6%) respectively. Archaea, fungi and viruses made up the minor communities in the overall microbial diversity. Interestingly, this study also identified a range of potential enteric zoonotic pathogens (0.02-1.25%) and genes involved in antimicrobial resistance (0.02-0.7%) in feline fecal materials. Based on clustering among nine gastrointestinal metagenomes from five different monogastric hosts (dog, human, mice, cat and chicken), the cat metagenome clustered closely together with chicken in both phylogenetic and metabolic level (>80%). Future studies are required to provide deeper understandings on both intrinsic and extrinsic effects such as impact of age, genetics and dietary interventions on the composition of the feline gastrointestinal microbiome.
| S-EPMC7114314 | biostudies-literature
Project description:BackgroundRNA viruses infecting a host usually exist as a set of closely related sequences, referred to as quasispecies. The genomic diversity of viral quasispecies is a subject of great interest, particularly for chronic infections, since it can lead to resistance to existing therapies. High-throughput sequencing is a promising approach to characterizing viral diversity, but unfortunately standard assembly software was originally designed for single genome assembly and cannot be used to simultaneously assemble and estimate the abundance of multiple closely related quasispecies sequences.ResultsIn this paper, we introduce a new Viral Spectrum Assembler (ViSpA) method for quasispecies spectrum reconstruction and compare it with the state-of-the-art ShoRAH tool on both simulated and real 454 pyrosequencing shotgun reads from HCV and HIV quasispecies. Experimental results show that ViSpA outperforms ShoRAH on simulated error-free reads, correctly assembling 10 out of 10 quasispecies and 29 sequences out of 40 quasispecies. While ShoRAH has a significant advantage over ViSpA on reads simulated with sequencing errors due to its advanced error correction algorithm, ViSpA is better at assembling the simulated reads after they have been corrected by ShoRAH. ViSpA also outperforms ShoRAH on real 454 reads. Indeed, 7 most frequent sequences reconstructed by ViSpA from a real HCV dataset are viable (do not contain internal stop codons), and the most frequent sequence was within 1% of the actual open reading frame obtained by cloning and Sanger sequencing. In contrast, only one of the sequences reconstructed by ShoRAH is viable. On a real HIV dataset, ShoRAH correctly inferred only 2 quasispecies sequences with at most 4 mismatches whereas ViSpA correctly reconstructed 5 quasispecies with at most 2 mismatches, and 2 out of 5 sequences were inferred without any mismatches. ViSpA source code is available at http://alla.cs.gsu.edu/~software/VISPA/vispa.html.ConclusionsViSpA enables accurate viral quasispecies spectrum reconstruction from 454 pyrosequencing reads. We are currently exploring extensions applicable to the analysis of high-throughput sequencing data from bacterial metagenomic samples and ecological samples of eukaryote populations.
| S-EPMC3194189 | biostudies-literature
Project description:BACKGROUND: The Himalayan or Indian Mayapple (Podophyllum hexandrum Royle) produces podophyllotoxin, which is used in the production of semisynthetic anticancer drugs. High throughput transcriptome sequences or genomic sequence data from the Indian Mayapple are essential for further understanding of the podophyllotoxin biosynthetic pathway. RESULTS: 454 pyrosequencing of a P. hexandrum cell culture normalized cDNA library generated 2,667,207 raw reads and 1,503,232 high quality reads, with an average read length of 138 bp. The denovo assembly was performed by Newbler using default and optimized parameters. The optimized parameter generated 40, 380 assembled sequences, comprising 12,940 contigs and 27,440 singlets which resulted in better assembly as compared to default parameters. BLASTX analysis resulted in the annotation of 40,380 contigs/singlet using a cut-off value of ? 1E-03. High similarity to Medicago truncatula using optimized parameters and to Populus trichocarpa using default parameters was noted. The Kyoto encyclopedia of genes and genomes (KEGG) analysis using KEGG Automatic Annotation Server (KAAS) combined with domain analysis of the assembled transcripts revealed putative members of secondary metabolism pathways that may be involved in podophyllotoxin biosynthesis. A proposed schematic pathway for phenylpropanoids and podophyllotoxin biosynthesis was generated. Expression profiling was carried out based on fragments per kilobase of exon per million fragments (FPKM). 1036 simple sequence repeats were predicted in the P. hexandrum sequences. Sixty-nine transcripts were mapped to 99 mature and precursor microRNAs from the plant microRNA database. Around 961 transcripts containing transcription factor domains were noted. High performance liquid chromatography analysis showed the peak accumulation of podophyllotoxin in 12-day cell suspension cultures. A comparative qRT-PCR analysis of phenylpropanoid pathway genes identified in the present data was performed to analyze their expression patterns in 12-day cell culture, callus and rhizome. CONCLUSIONS: The present data will help the identification of the potential genes and transcription factors involved in podophyllotoxin biosynthesis in P. hexandrum. The assembled transcripts could serve as potential candidates for marker discovery and conservation, which should form the foundations for future endeavors.
| S-EPMC3840631 | biostudies-literature
Project description:BackgroundWeevils of the genus Otiorhynchus are regarded as devastating pests in a wide variety of horticultural crops worldwide. So far, little is known on the presence of endosymbionts in Otiorhynchus spp.. Investigation of endosymbiosis in this genus may help to understand the evolution of different reproductive strategies in these weevils (parthenogenesis or sexual reproduction), host-symbiont interactions, and may provide a future basis for novel pest management strategy development. Here, we used a multitag 454 pyrosequencing approach to assess the bacterial endosymbiont diversity in larvae of four economically important Otiorhynchus species.ResultsHigh-throughput tag-encoded FLX amplicon pyrosequencing of a bacterial 16S rDNA fragment was used to characterise bacterial communities associated with different Otiorhynchus spp. larvae. By sequencing a total of ~48,000 PCR amplicons, we identified 49 different operational taxonomic units (OTUs) as bacterial endosymbionts in the four studied Otiorhynchus species. More than 90% of all sequence reads belonged either to the genus Rickettsia or showed homology to the phylogenetic group of "Candidatus Blochmannia" and to endosymbionts of the lice Pedicinus obtusus and P. badii. By using specific primers for the genera Rickettsia and "Candidatus Blochmannia", we identified a new phylogenetic clade of Rickettsia as well as "Candidatus Nardonella" endosymbionts in Otiorhynchus spp. which are closely related to "Candidatus Blochmannia" bacteria.ConclusionsHere, we used multitag 454 pyrosequencing for assessment of insect endosymbiotic communities in weevils. As 454 pyrosequencing generates only quite short sequences, results of such studies can be regarded as a first step towards identifying respective endosymbiotic species in insects. In the second step of our study, we analysed sequences of specific gene regions for a more detailed phylogeny of selected endosymbiont genera. As a result we identified the presence of Rickettsia and "Candidatus Nardonella" endosymbionts in Otiorhynchus spp.. This knowledge is an important step in exploring bacteria-insect associations for potential use in insect pest control.
| S-EPMC3287517 | biostudies-literature
Project description:Viruses, the most abundant biological entities on the planet, are capable of infecting organisms from all three branches of life, although the majority infect bacteria where the greatest degree of cellular diversity lies. However, the characterization and assessment of viral diversity in natural environments is only beginning to become a possibility. Through the development of a novel technique for the harvest of viral DNA and the application of 454 pyrosequencing, a snapshot of the diversity of the DNA viruses harvested from a standing pond on a cattle farm has been obtained. A high abundance of viral genotypes (785) were present within the virome. The absolute numbers of lambdoid and Shiga toxin (Stx) encoding phages detected suggested that the depth of sequencing had enabled recovery of only ca. 8% of the total virus population, numbers that agreed within less than an order of magnitude with predictions made by rarefaction analysis. The most abundant viral genotypes in the pond were bacteriophages (93.7%). The predominant viral genotypes infecting higher life forms found in association with the farm were pathogens that cause disease in cattle and humans, e.g. members of the Herpesviridae. The techniques and analysis described here provide a fresh approach to the monitoring of viral populations in the aquatic environment, with the potential to become integral to the development of risk analysis tools for monitoring the dissemination of viral agents of animal, plant and human diseases.
| S-EPMC3954088 | biostudies-literature
Project description:BackgroundThe Antarctic krill Euphausia superba is a keystone species in the Antarctic food chain. Not only is it a significant grazer of phytoplankton, but it is also a major food item for charismatic megafauna such as whales and seals and an important Southern Ocean fisheries crop. Ecological data suggest that this species is being affected by climate change and this will have considerable consequences for the balance of the Southern Ocean ecosystem. Hence, understanding how this organism functions is a priority area and will provide fundamental data for life history studies, energy budget calculations and food web models.Methodology/principal findingsThe assembly of the 454 transcriptome of E. superba resulted in 22,177 contigs with an average size of 492bp (ranging between 137 and 8515bp). In depth analysis of the data revealed an extensive catalogue of the cellular chaperone systems and the major antioxidant proteins. Full length sequences were characterised for the chaperones HSP70, HSP90 and the super-oxide dismutase antioxidants, with the discovery of potentially novel duplications of these genes. The sequence data contained 41,470 microsatellites and 17,776 Single Nucleotide Polymorphisms (SNPs/INDELS), providing a resource for population and also gene function studies.ConclusionsThis paper details the first 454 generated data for a pelagic Antarctic species or any pelagic crustacean globally. The classical "stress proteins", such as HSP70, HSP90, ferritin and GST were all highly expressed. These genes were shown to be over expressed in the transcriptomes of Antarctic notothenioid fish and hypothesized as adaptations to living in the cold, with the associated problems of decreased protein folding efficiency and increased vulnerability to damage by reactive oxygen species. Hence, these data will provide a major resource for future physiological work on krill, but in particular a suite of "stress" genes for studies understanding marine ectotherms' capacities to cope with environmental change.
| S-EPMC3017093 | biostudies-literature