Dataset Information

Designing deep sequencing experiments: detecting structural variation and estimating transcript abundance.

ABSTRACT:

Background

Massively parallel DNA sequencing technologies have enabled the sequencing of several individual human genomes. These technologies are also being used in novel ways for mRNA expression profiling, genome-wide discovery of transcription-factor binding sites, small RNA discovery, etc. The multitude of sequencing platforms, each with their unique characteristics, pose a number of design challenges, regarding the technology to be used and the depth of sequencing required for a particular sequencing application. Here we describe a number of analytical and empirical results to address design questions for two applications: detection of structural variations from paired-end sequencing and estimating mRNA transcript abundance.

Results

For structural variation, our results provide explicit trade-offs between the detection and resolution of rearrangement breakpoints, and the optimal mix of paired-read insert lengths. Specifically, we prove that optimal detection and resolution of breakpoints is achieved using a mix of exactly two insert library lengths. Furthermore, we derive explicit formulae to determine these insert length combinations, enabling a 15% improvement in breakpoint detection at the same experimental cost. On empirical short read data, these predictions show good concordance with Illumina 200 bp and 2 Kbp insert length libraries. For transcriptome sequencing, we determine the sequencing depth needed to detect rare transcripts from a small pilot study. With only 1 Million reads, we derive corrections that enable almost perfect prediction of the underlying expression probability distribution, and use this to predict the sequencing depth required to detect low expressed genes with greater than 95% probability.

Conclusions

Together, our results form a generic framework for many design considerations related to high-throughput sequencing. We provide software tools http://bix.ucsd.edu/projects/NGS-DesignTools to derive platform independent guidelines for designing sequencing experiments (amount of sequencing, choice of insert length, mix of libraries) for novel applications of next generation sequencing.

SUBMITTER: Bashir A

PROVIDER: S-EPMC3091630 | biostudies-literature | 2010 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Designing deep sequencing experiments: detecting structural variation and estimating transcript abundance.

Bashir Ali A Bansal Vikas V Bafna Vineet V

BMC genomics 20100618

<h4>Background</h4>Massively parallel DNA sequencing technologies have enabled the sequencing of several individual human genomes. These technologies are also being used in novel ways for mRNA expression profiling, genome-wide discovery of transcription-factor binding sites, small RNA discovery, etc. The multitude of sequencing platforms, each with their unique characteristics, pose a number of design challenges, regarding the technology to be used and the depth of sequencing required for a part ...[more]

PMID: 20565853

Similar Datasets

Project description:Deep-sequencing allows for an in-depth characterization of sequence variation in complex populations. However, technology associated errors may impede a powerful assessment of low-frequency mutations. Fortunately, base calls are complemented with quality scores which are derived from a quadruplet of intensities, one channel for each nucleotide type for Illumina sequencing. The highest intensity of the four channels determines the base that is called. Mismatch bases can often be corrected by the second best base, i.e. the base with the second highest intensity in the quadruplet. A virus variant model-based clustering method, ViVaMBC, is presented that explores quality scores and second best base calls for identifying and quantifying viral variants. ViVaMBC is optimized to call variants at the codon level (nucleotide triplets) which enables immediate biological interpretation of the variants with respect to their antiviral drug responses.Using mixtures of HCV plasmids we show that our method accurately estimates frequencies down to 0.5%. The estimates are unbiased when average coverages of 25,000 are reached. A comparison with the SNP-callers V-Phaser2, ShoRAH, and LoFreq shows that ViVaMBC has a superb sensitivity and specificity for variants with frequencies above 0.4%. Unlike the competitors, ViVaMBC reports a higher number of false-positive findings with frequencies below 0.4% which might partially originate from picking up artificial variants introduced by errors in the sample and library preparation step.ViVaMBC is the first method to call viral variants directly at the codon level. The strength of the approach lies in modeling the error probabilities based on the quality scores. Although the use of second best base calls appeared very promising in our data exploration phase, their utility was limited. They provided a slight increase in sensitivity, which however does not warrant the additional computational cost of running the offline base caller. Apparently a lot of information is already contained in the quality scores enabling the model based clustering procedure to adjust the majority of the sequencing errors. Overall the sensitivity of ViVaMBC is such that technical constraints like PCR errors start to form the bottleneck for low frequency variant detection.

Project description:Phakopsora pachyrhizi, the causal agent responsible for soybean rust, is among the top hundred most virulent plant pathogens and can cause soybean yield losses of up to 80% when appropriate conditions are met. We used mRNA-Seq by Illumina to analyze pathogen transcript abundance at 15 seconds (s), 7 hours (h), 48 h, and 10 days (d) after inoculation (ai) of susceptible soybean leaves with P. pachyrhizi to gain new insights into transcript abundance in soybean and the pathogen at specific time-points during the infection including the uredinial stage.Over three million five hundred thousand sequences were obtained for each time-point. Energy, nucleotide metabolism, and protein synthesis are major priorities for the fungus during infection and development as indicated by our transcript abundance studies. At all time-points, energy production is a necessity for P. pachyrhizi, as indicated by expression of many transcripts encoding enzymes involved in oxidative phosphorylation and carbohydrate metabolism (glycolysis, glyoxylate and dicarboxylate, pentose phosphate, pyruvate). However, at 15 sai, transcripts encoding enzymes involved in ATP production were highly abundant in order to provide enough energy for the spore to germinate, as observed by the expression of many transcripts encoding proteins involved in electron transport. At this early time-point, transcripts encoding proteins involved in RNA synthesis were also highly abundant, more so than transcripts encoding genes involved in DNA and protein synthesis. At 7 hai, shortly after germination during tube elongation and penetration, transcripts encoding enzymes involved in deoxyribonucleotide and DNA synthesis were highly abundant. At 48 hai, transcripts encoding enzymes involved in amino acid metabolism were highly abundant to provide for increased protein synthesis during haustoria maturation. During sporulation at 10 dai, the fungus still required carbohydrate metabolism, but there also was increased expression of transcripts encoding enzymes involved in fatty acid metabolism.This information provides insight into molecular events and their timing throughout the life cycle of the P. pachyrhizi, and it may be useful in the development of new methods of broadening resistance of soybean to soybean rust.

Dataset Information

Designing deep sequencing experiments: detecting structural variation and estimating transcript abundance.

Background

Results

Conclusions

Publications

Designing deep sequencing experiments: detecting structural variation and estimating transcript abundance.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure