14 research outputs found
Basic components of the iterative workflow as compared to a standard NGS whole genome analysis.
<p>Basic components of the iterative workflow as compared to a standard NGS whole genome analysis.</p
Concordance of SNP data with variants from standard and iterative workflows for sample NA12878.
<p>Concordance of SNP data with variants from standard and iterative workflows for sample NA12878.</p
Evaluation of SNVs and Indels called by the iterative and standard workflow.
<p>Evaluation of SNVs and Indels called by the iterative and standard workflow.</p
SoftSearch: Integration of Multiple Sequence Features to Identify Breakpoints of Structural Variations
<div><p>Background</p><p>Structural variation (SV) represents a significant, yet poorly understood contribution to an individual’s genetic makeup. Advanced next-generation sequencing technologies are widely used to discover such variations, but there is no single detection tool that is considered a community standard. In an attempt to fulfil this need, we developed an algorithm, SoftSearch, for discovering structural variant breakpoints in Illumina paired-end next-generation sequencing data. SoftSearch combines multiple strategies for detecting SV including split-read, discordant read-pair, and unmated pairs. Co-localized split-reads and discordant read pairs are used to refine the breakpoints. </p> <p>Results</p><p>We developed and validated SoftSearch using real and synthetic datasets. SoftSearch’s key features are 1) not requiring secondary (or exhaustive primary) alignment, 2) portability into established sequencing workflows, and 3) is applicable to any DNA-sequencing experiment (e.g. whole genome, exome, custom capture, etc.). SoftSearch identifies breakpoints from a small number of soft-clipped bases from split reads and a few discordant read-pairs which on their own would not be sufficient to make an SV call. </p> <p>Conclusions</p><p>We show that SoftSearch can identify more true SVs by combining multiple sequence features. SoftSearch was able to call clinically relevant SVs in the BRCA2 gene not reported by other tools while offering significantly improved overall performance.</p> </div
Example IGV screenshot of a 71bp tandem duplication in the BRCA2 gene identified by SoftSearch.
<p>Discordant reads are blue (plus strand) or red (minus strand). Soft clipped bases appear as multicolour “rainbows”.</p
Overlap of true positive calls for the NA12878 and NA18507 datasets.
<p>Overlap of true positive calls for the NA12878 and NA18507 datasets.</p
The general strategy for SoftSearch.
<p>A) Left clipped reads are defined as where the clipped portion of the read is at a smaller genome coordinate than the opposite end (opposite for right clipping). For a left clipped read located on the “+” strand, SoftSearch looks upstream for a discordant read pair where the read is oriented in the “-” direction. The orientation and location of the mate is where SoftSearch links the first region to. To increase the likelihood of exactly detecting the breakpoint, it then looks upstream for a right clipped read cluster. If none is found, then the default breakpoint location is the discordant read mate location; otherwise it is the position of soft clipping at the right clipped read. B) SoftSearch determines discordant read pairs by their insert size and orientation and places them in a temporary BAM file. It also reads the input BAM file for soft clipped reads and converts them to a BED file. Overlapping soft clip locations are counted to identify putative breakpoints, and then queried against the discordant read pair bam file for properly oriented supporting reads, which are then output in VCF format.</p
Visualization of single nucleotide variant validation.
<p>Sanger sequence validation of highly expressed novel somatic SNVs for <i>MRPL3</i> variant in the BCT40 HER2 tumor sample. RNA-Seq sequence reads shown above Sanger sequencing tracing, with mutation shown by an arrow.</p
Visualization of single nucleotide variant validation.
<p>Sanger sequencing validation of MPG eSNV in HER2 tumor. RNA-Seq reads shown over Sanger sequence tracing with mutation indicated by an arrow.</p
HER2-positive tumor interactome.
<p>HER2 tumor interactome network developed using cytoscape and functional reactome of 1055 genes obtained from integration analyses of genomic features. Different functional modules within the network are color coded.</p