LongISLND: in silico sequencing of lengthy and noisy datatypes.

LongISLND is a software package designed to simulate sequencing data according to the characteristics of third generation, single-molecule sequencing technologies. The general software architecture is easily extendable, as demonstrated by the emulation of Pacific Biosciences (PacBio) multi-pass sequencing with P5 and P6 chemistries, producing data in FASTQ, H5, and the latest PacBio BAM format. We demonstrate its utility by downstream processing with consensus building and variant calling. AVAILABILITY AND IMPLEMENTATION: LongISLND is implemented in Java and available at http://bioinform.github.io/longislnd CONTACT: hugo.lam@roche.comSupplementary information: Supplementary data are available at Bioinformatics online.

An ensemble approach to accurately detect somatic mutations using SomaticSeq.

SomaticSeq is an accurate somatic mutation detection pipeline implementing a stochastic boosting algorithm to produce highly accurate somatic mutation calls for both single nucleotide variants and small insertions and deletions. The workflow currently incorporates five state-of-the-art somatic mutation callers, and extracts over 70 individual genomic and sequencing features for each candidate site. A training set is provided to an adaptively boosted decision tree learner to create a classifier for predicting mutation statuses. We validate our results with both synthetic and real data. We report that SomaticSeq is able to achieve better overall accuracy than any individual tool incorporated.

Leveraging long read sequencing from a single individual to provide a comprehensive resource for benchmarking variant calling methods.

A high-confidence, comprehensive human variant set is critical in assessing accuracy of sequencing algorithms, which are crucial in precision medicine based on high-throughput sequencing. Although recent works have attempted to provide such a resource, they still do not encompass all major types of variants including structural variants (SVs). Thus, we leveraged the massive high-quality Sanger sequences from the HuRef genome to construct by far the most comprehensive gold set of a single individual, which was cross validated with deep Illumina sequencing, population datasets, and well-established algorithms. It was a necessary effort to completely reanalyze the HuRef genome as its previously published variants were mostly reported five years ago, suffering from compatibility, organization, and accuracy issues that prevent their direct use in benchmarking. Our extensive analysis and validation resulted in a gold set with high specificity and sensitivity. In contrast to the current gold sets of the NA12878 or HS1011 genomes, our gold set is the first that includes small variants, deletion SVs and insertion SVs up to a hundred thousand base-pairs. We demonstrate the utility of our HuRef gold set to benchmark several published SV detection tools.

MetaSV: an accurate and integrative structural-variant caller for next generation sequencing.

UNLABELLED: Structural variations (SVs) are large genomic rearrangements that vary significantly in size, making them challenging to detect with the relatively short reads from next-generation sequencing (NGS). Different SV detection methods have been developed; however, each is limited to specific kinds of SVs with varying accuracy and resolution. Previous works have attempted to combine different methods, but they still suffer from poor accuracy particularly for insertions. We propose MetaSV, an integrated SV caller which leverages multiple orthogonal SV signals for high accuracy and resolution. MetaSV proceeds by merging SVs from multiple tools for all types of SVs. It also analyzes soft-clipped reads from alignment to detect insertions accurately since existing tools underestimate insertion SVs. Local assembly in combination with dynamic programming is used to improve breakpoint resolution. Paired-end and coverage information is used to predict SV genotypes. Using simulation and experimental data, we demonstrate the effectiveness of MetaSV across various SV types and sizes. AVAILABILITY AND IMPLEMENTATION: Code in Python is at http://bioinform.github.io/metasv/. CONTACT: rd@bina.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VarSim: a high-fidelity simulation and validation framework for high-throughput genome sequencing with cancer applications.

SUMMARY: VarSim is a framework for assessing alignment and variant calling accuracy in high-throughput genome sequencing through simulation or real data. In contrast to simulating a random mutation spectrum, it synthesizes diploid genomes with germline and somatic mutations based on a realistic model. This model leverages information such as previously reported mutations to make the synthetic genomes biologically relevant. VarSim simulates and validates a wide range of variants, including single nucleotide variants, small indels and large structural variants. It is an automated, comprehensive compute framework supporting parallel computation and multiple read simulators. Furthermore, we developed a novel map data structure to validate read alignments, a strategy to compare variants binned in size ranges and a lightweight, interactive, graphical report to visualize validation results with detailed statistics. Thus far, it is the most comprehensive validation tool for secondary analysis in next generation sequencing. AVAILABILITY AND IMPLEMENTATION: Code in Java and Python along with instructions to download the reads and variants is at http://bioinform.github.io/varsim. CONTACT: rd@bina.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploring the physiologic role of human gastroesophageal reflux by analyzing time-series data from 24-h gastric and esophageal pH recordings.

Our previous finding of a fractal pattern for gastric pH and esophageal pH plus the statistical association of sequential pH values for up to 2 h led to our hypothesis that the fractal pattern encodes information regarding gastric acidity and that depending on the value of gastric acidity, the esophagus can signal the stomach to alter gastric acidity by influencing gastric secretion of acid or bicarbonate. Under our hypothesis values of gastric pH should provide information regarding values of esophageal pH and vice versa. We used vector autoregression, a theory-free set of inter-related linear regressions used to measure relationships that can change over time, to analyze data from 24-h recordings of gastric pH and esophageal pH. We found that in pH records from normal subjects, as well as from subjects with gastroesophageal reflux disease alone and after treatment with a proton pump inhibitor, gastric pH values provided important information regarding subsequent values of esophageal pH and values of esophageal pH provided important information regarding subsequent values of gastric pH. The ability of gastric pH and esophageal pH to provide information regarding subsequent values of each other was reduced in subjects with gastroesophageal reflux disease compared to normal subjects. Our findings are consistent with the hypothesis that depending on the value of gastric acidity, the esophagus can signal the stomach to alter gastric acidity, and that this ability is impaired in subjects with gastroesophageal reflux disease.

Fast and accurate read alignment for resequencing.

MOTIVATION: Next-generation sequence analysis has become an important task both in laboratory and clinical settings. A key stage in the majority sequence analysis workflows, such as resequencing, is the alignment of genomic reads to a reference genome. The accurate alignment of reads with large indels is a computationally challenging task for researchers. RESULTS: We introduce SeqAlto as a new algorithm for read alignment. For reads longer than or equal to 100 bp, SeqAlto is up to 10 × faster than existing algorithms, while retaining high accuracy and the ability to align reads with large (up to 50 bp) indels. This improvement in efficiency is particularly important in the analysis of future sequencing data where the number of reads approaches many billions. Furthermore, SeqAlto uses less than 8 GB of memory to align against the human genome. SeqAlto is benchmarked against several existing tools with both real and simulated data. AVAILABILITY: Linux and Mac OS X binaries free for academic use are available at http://www.stanford.edu/group/wonglab/seqalto CONTACT: whwong@stanford.edu.

Mutation of a U2 snRNA gene causes global disruption of alternative splicing and neurodegeneration.

Although uridine-rich small nuclear RNAs (U-snRNAs) are essential for pre-mRNA splicing, little is known regarding their function in the regulation of alternative splicing or of the biological consequences of their dysfunction in mammals. Here, we demonstrate that mutation of Rnu2-8, one of the mouse multicopy U2 snRNA genes, causes ataxia and neurodegeneration. Coincident with the observed pathology, the level of mutant U2 RNAs was highest in the cerebellum and increased after granule neuron maturation. Furthermore, neuron loss was strongly dependent on the dosage of mutant and wild-type snRNA genes. Comprehensive transcriptome analysis identified a group of alternative splicing events, including the splicing of small introns, which were disrupted in the mutant cerebellum. Our results suggest that the expression of mammalian U2 snRNA genes, previously presumed to be ubiquitous, is spatially and temporally regulated, and dysfunction of a single U2 snRNA causes neuron degeneration through distortion of pre-mRNA splicing.