- Technical Note
SeqBench: Integrated solution for the management and analysis of exome sequencing data
BMC Research Notesvolume 7, Article number: 43 (2014)
The rapid development of next generation sequencing technologies, including the recently introduced benchtop sequencers, made sequencing affordable for smaller research institutions. A widely applied method to identify causing mutations of diseases is exome sequencing, which proved to be cost-effective and time-saving.
SeqBench, a web-based application, combines management and analysis of exome sequencing data into one solution. It provides a user friendly data acquisition module to facilitate comprehensive and intuitive data handling. SeqBench provides direct access to the analysis pipeline SIMPLEX, which can be configured to run locally, on a cluster, or in the cloud. Identified genomic variants are presented along with several functional annotations and can be interpreted in a family context.
The web-based application SeqBench supports the management and analysis of exome sequencing data, is open-source and available athttp://www.icbi.at/SeqBench.
The introduction of next-generation sequencing (NGS) technologies has sparked a major revolution in biomedical research. Modern benchtop sequencers - such as the MiSeq (Illumina), PGM (life technologies), and 454 GS Junior (Roche) - generate only modest set-up and maintenance costs, and as a consequence genomic sequencing data are now produced at an unprecedented rate. At the same time, they are confronted with enormous challenges in data management, computer hardware requirements, and data analysis.
Exome sequencing is a widely used approach to elucidate the causing mutations of various diseases, including Mendelian disorders and cancer[1–3]. A plethora of tools for analyzing exome sequencing data are available, making efficient and standardized data analysis challenging for inexperienced users.
SIMPLEX is a cloud-enabled analysis pipeline, which autonomously performs the complete workflow of exome sequencing data analysis. The pipeline is able to perform the following steps: (a) initial quality control, (b) intelligent data filtering and pre-processing, (c) sequence alignment to a reference genome, (d) SNP and INDEL detection, (e) functional annotation of variants, and (f) detailed report generation during various stages of the workflow.
SIMPLEX combines several established analysis applications (such as BWA and GATK) with in-house developed analysis tools (e.g. the quality report) into a tailored pipeline. Underlying applications can be configured using a variety of exposed parameters. The pipeline annotates identified variants using GATK, ANNOVAR, and SnpEff as well as a variety of attributes derived from the databases RefSeq, GO, KEGG, dbSNP, and the 1000 genomes project. Additionally, annotation scores from SIFT, PolyPhen2, PhyloP, MutationTaster, and LRT are added to the variants.
However, to the best of our knowledge there is currently no tool available which supports cloud computing and combines data management, data analysis, and result visualization into a single solution. Such an integrated solution would be of great help for labs without dedicated computational infrastructure. Additionally, the system should be easy to install and therefore be provided as a virtual machine (VM) image.
The new web-based application SeqBench embeds the established analysis pipeline SIMPLEX and provides a data acquisition module supporting data derived from Illumina and SOLiD platforms. This module guides the user through all stages of data entry. In addition, SeqBench provides a module to monitor ongoing SIMPLEX analyses and offers an intuitive interface to interpret computed results.
The preconfigured out of the box solution is very easy to install and presents itself as the ideal solution for inexperienced users. Yet, as SeqBench is completely open-source, it can be adapted and modified by experienced users to answer specific biological questions. Furthermore, SeqBench offers multiple customization options and can be configured to run the computationally expensive data analysis pipeline either in the cloud, on a high-performance computing infrastructure, or on a powerful local desktop machine using the provided VM image.
SeqBench is organized in different modules as depicted in Figure1. After logging in, a dashboard shows currently running SIMPLEX processes, results of finished analyses, and a list of available projects. New and existing projects can be added and modified using an intuitive data acquisition wizard. Users are guided through all necessary steps, making sure that project information and complementary data is entered in a consistent and coherent way. SeqBench allows users to upload and analyze raw sequencing data in FASTQ (Illumina) or CSFASTA/QUAL (SOLiD) format. As the analysis pipeline offers a wide variety of parameters to tailor the analysis to a specific biological problem, SeqBench offers a flexible system to set those parameters.
A well chosen set of standard values is based on the suggestions by the used applications of SIMPLEX. SeqBench passes the uploaded raw files and selected parameter values on to SIMPLEX thus starting the analysis process. As output of the pipeline, files for several steps are created, SeqBench constantly monitors those output files and displays the progress in the result view. At the end of a pipeline run a summary of result files is generated, whereby selected files can be directly accessed via SeqBench.
The result visualization module in SeqBench allows users to access and download all generated analysis files. Furthermore, a table of identified SNPs and INDELs is displayed, including several functional annotations with outgoing links. Furthermore, all annotations found by SIMPLEX can be displayed as columns of this table. This table features dynamic column selection, real-time filtering, and sorting of mutations. Descriptive statistics about the filtered variants are illustrated in charts as depicted in Figure2, which are updated dynamically after filtering. As SIMPLEX supports the analysis of multiple samples at once, SeqBench provides a combined view of identified variants to interpret results in their family context. Therefore, after uploading all files into SeqBench, one SIMPLEX analysis run can be started including all samples. Next, identified variants can be loaded into the variant table view, where they can be filtered and investigated in a combined way.
Access to the system is secured by an authorization and authentication system (AAS), which can be configured to use either a simple XML file or a web-based application for the usermanagement. The AAS allows the assignment of different user roles with dedicated data access permissions, which are enforced within SeqBench.
SeqBench is open-source and distributed as a VM image that ships with a preconfigured AAS and a fully configured SIMPLEX pipeline. A detailed user manual provides information about installation, configuration, and usage of SeqBench as well as the installation process using VirtualBox as the virtual machine software which runs on Windows, Linux and Mac OS X (https://seqbench.i-med.ac.at/SeqBench/doc/userManual.jsf). Details about how to use SIMPLEX in the cloud is provided in the SIMPLEX user manual (http://icbi.at/software/simplex/simplex.shtml).
SeqBench is a Java EE 6 based web application running on a JBoss application server. It uses a relational database (PostgreSQL) accessed via an object relational mapper (Hibernate). JBoss Seam in combination with "Contexts and Dependency Injection" (CDI) is used as underlying framework. The presentation layer is based on Java Server Faces (JSF) and the JBoss Richfaces component library, which supports internationalization and input validation.
SeqBench combines data management, data analysis, and result visualization into a single web-based platform. The complete software suite is distributed as a VM image providing several benefits to users such as easy installation, replication of analyses, or taking snapshots for archiving results (see recently published review). To the best of our knowledge it is the only integrated solution that provides, amongst other features, SOLiD as well as cloud support.
SeqBench provides an intuitive user interface and facilitates comprehensive data acquisition through a wizard system and can be easily installed as an out of the box solution. The application covers the complete workflow from data acquisition to result interpretation and supports Illumina as well as SOLiD data. In conclusion, SeqBench is a light-weight and effective solution for the streamlined integration of data management and analyses of exome sequencing data.
Availability and requirements
AD, SP, MS, and MF developed SeqBench. AD designed the database scheme. SP and GS implemented the front-end of the application. MS developed the results report view. MF integrated SIMPLEX into SeqBench. GS and ZT conceived the project and suppervised the development. All authors participated in debugging and testing of SeqBench, and assisted in preparation of the article. All authors read and approved the final manuscript.
Authorization and Authentication System
Affero General Public License
Kyoto Encyclopedia of Genes and Genomes
Schossig A, Wolf NI, Fischer C, Fischer M, Stocker G, Pabinger S, Dander A, Steiner B, Tönz O, Kotzot D, Haberlandt E, Amberger A, Burwinkel B, Wimmer K, Fauth C, Grond-Ginsbach C, Koch MJ, Deichmann A, von Kalle C, Bartram CR, Kohlschütter A, Trajanoski Z, Zschocke J: Mutations in ROGDI Cause Kohlschütter-Tönz Syndrome. Am J Hum Genet. 2012, 90 (4): 701-707. 10.1016/j.ajhg.2012.02.012. [PMID: 22424600]
Ishida S, Picard F, Rudolf G, Noé E, Achaz G, Thomas P, Genton P, Mundwiller E, Wolff M, Marescaux C, Miles R, Baulac M, Hirsch E, Leguern E, Baulac S: Mutations of DEPDC5 cause autosomal dominant focal epilepsies. Nat Genet. 2013, 45: 552-555. 10.1038/ng.2601. [PMID: 23542701]
Dulak AM, Stojanov P, Peng S, Lawrence MS, Fox C, Stewart C, Bandla S, Imamura Y, Schumacher SE, Shefler E, McKenna A, Carter SL, Cibulskis K, Sivachenko A, Saksena G, Voet D, Ramos AH, Auclair D, Thompson K, Sougnez C, Onofrio RC, Guiducci C, Beroukhim R, Zhou Z, Lin L, Lin J, Reddy R, Chang A, Landrenau R, Pennathur A, et al: Exome and whole-genome sequencing of esophageal adenocarcinoma identifies recurrent driver events and mutational complexity. Nat Genet. 2013, 45: 478-486. 10.1038/ng.2591. [PMID: 23525077]
Pabinger S, Dander A, Fischer M, Snajder R, Sperk M, Efremova M, Krabichler B, Speicher MR, Zschocke J, Trajanoski Z: A survey of tools for variant analysis of next-generation genome sequencing data. Brief Bioinform. 2013, [PMID: 23341494]
Fischer M, Snajder R, Pabinger S, Dander A, Schossig A, Zschocke J, Trajanoski Z, Stocker G: SIMPLEX: cloud-enabled pipeline for the comprehensive analysis of exome sequencing data. PloS One. 2012, 7 (8): e41948-10.1371/journal.pone.0041948. [PMID: 22870267]
Mariette J, Escudié F, Allias N, Salin G, Noirot C, Thomas S, Klopp C: NG6: Integrated next generation sequencing storage and processing environment. BMC Genomics. 2012, 13: 462-10.1186/1471-2164-13-462. [PMID: 22958229]
Scholtalbers J, Rößler J, Sorn P, Graaf Jd, Boisguérin V, Castle J, Sahin U: Galaxy LIMS for next-generation sequencing. Bioinformatics. 2013, 29 (9): 1233-1234. 10.1093/bioinformatics/btt115. [PMID: 23479349]
Nix DA, Di Sera TL, Dalley BK, Milash BA, Cundick RM, Quinn KS, Courdy SJ: Next generation tools for genomic data generation, distribution, and visualization. BMC Bioinformatics. 2010, 11: 455-10.1186/1471-2105-11-455. [PMID: 20828407]
Bauch A, Adamczyk I, Buczek P, Elmer FJ, Enimanev K, Glyzewski P, Kohler M, Pylak T, Quandt A, Ramakrishnan C, Beisel C, Malmström L, Aebersold R, Rinn B: openBIS: a flexible framework for managing and analyzing complex data in biology research. BMC Bioinformatics. 2011, 12: 468-10.1186/1471-2105-12-468. [PMID: 22151573]
Pabinger S, Rader R, Agren R, Nielsen J, Trajanoski Z: MEMOSys: Bioinformatics platform for genome-scale metabolic models. BMC Syst Biol. 2011, 5: 20-10.1186/1752-0509-5-20. [PMID: 21276275]
Nocq J, Celton M, Gendron P, Lemieux S, Wilhelm BT: Harnessing virtual machines to simplify next generation DNA sequencing analysis. Bioinformatics. 2013, 29 (17): 2075-2083. 10.1093/bioinformatics/btt352. [PMID: 23786767]
We would like to thank Rene Snajder and Mark Bessem for their contributions to this project.
Funding: This work was supported by ONCOTYROL (FFG, Austrian Federal Ministries BMVIT/BMWFJ), the Austrian Drug Screening Institute, and a grant from the Standortagentur Tirol.
The authors declare that they have no competing interests.
Andreas Dander, Stephan Pabinger contributed equally to this work.