BACKGROUND
next-generation sequencing  has become the method of choice for molecular genetic analyses such as transcriptome profiling , chromatin immunoprecipitation followed by sequencing  as well as resequencing of complete genomes.
numerous software packages designed to analyze massive amounts of ngs data are now publicly available.
preprocessing of ngs data typically takes advantage of a complex hardware architecture composed of, for instance, a parallel compute cluster, a database server and a web server.
as this requires specialized it skills, the widespread access to ngs technology is still hampered by technical barriers.
the primary data analysis is therefore often centralized into core laboratories which face the challenge of using a reasonable selection out of the available software packages to process a growing flow of newly produced data.
we introduce a new framework named quickngs which can be operated by it staff in bioinformatics core labs to process vast amounts of data provided by their end users, typically experimental scientists.
quickngs enables a major decrease in time and effort put into the primary analysis of ngs data, thus contributing to the further evolution of ngs into a standard tool with high accessibility to researchers.
implementation
quickngs enables rapid and professional data analysis for the aforementioned major applications of ngs in a batch-like operation mode.
the core of the quickngs workflow is formed by a comprehensive mysql database used to organize sample metadata such as species, library type and batch information as well as the analysis results.
the system operators can use the quickngs database interface  to upload metadata and monitor the status of quickngs analysis runs from any location .
on the other hand, the purpose of the quickngs web interface is to deliver web reports on the analysis results to the end users who can use their personal login and password to browse the results plus a lot of useful on-top information .
the core results are provided as spreadsheet files which the user can download to a local computer.fig.
<dig>
a the quickngs database contains meta information on samples  and sample groups .
both can be efficiently organized by an intuitive web interface.
new samples and sample groups can be inserted  by following the „ + “button.
b the status page on the quickngs database monitors time, user information and current status of each quickngs module on a clearly arranged website, enabling password-protected interrogation of the current status from any working location, including mobile access
fig.
<dig> results report for the test run on  <dig> drosophila rna-seq samples.
after login, the interface provides links to the main database export files, qc statistics as well as visualisation plots and the user’s personal password-protected ucsc genome browser track hub
as the sample metadata in the quickngs database are used to completely control the overall workflow, these have to be provided to the quickngs database before any analysis can be started.
to achieve this, the file locations of the raw data first need to be saved into a text file .
this file is then uploaded into the quickngs web interface  together with information on the library type, ngs application type, species of origin and the laboratory in which the input material has been generated.
for each sample listed in the text file, the user is then asked to provide a human-readable sample label as well as batch information for the case that samples have been processed with different library preparation protocols or in different ngs runs .
finally, samples can be assigned to pairs and groups for comparative analysis, e.g.
differential gene expression .
subsequently, the raw data received from the sequencing center are linked into the quickngs stack directory which is then processed fully automatically.
thus, the only user action needed is providing sample information and linking the files to the stack directory which both are trivial amounts of effort.
the results produced by the workflow are saved back into the quickngs database and made accessible through a report on the web interface.
this report comprises standard qc metrics  and results on typical data-related research questions, for instance which genes are differentially expressed or differentially spliced, which genomic variations are unique to a sample compared to a control, which transcription factor binding motifs are enriched in a chip-seq data set etc.. the results and report are generated fully automatically without any additional user action.
the analysis relies on widely adopted ngs analysis packages which are listed in table  <dig>  for the core analysis of the raw data, we have carefully selected the most appropriate previously published software programs.
the selection criteria were  performance in published and in-house benchmarking studies,  comprehensiveness of the analysis output,  quality of the implementation and steadiness of maintenance, and  popularity in the community.
our choice of bioinformatics software follows these criteria as far as possible.
the code for qc and visualisation as well as for data management and the workflow itself is unique to quickngs.
as a reference to genomic sequence and annotation, the system uses the mirbase  <cit>  for the mirna-seq workflow and the ensembl database  <cit>  for all other applications.
for instance, rna-seq or chip-seq analyses can thus be carried out on data from any arbitrary organism listed on either ensembl  or ensemblgenomes .
the reference files are downloaded to the local system and updated automatically.
the same applies to genomic annotation data which are retrieved using biomart  <cit> .table  <dig> algorithms and software tools used by quickngs, version  <dig> .
<dig> .
the selection may be modified, updated or extended in future releases of quickngs, however, an up-to-date version of this table will be kept available online at the quickngs website
RESULTS
test run on previously published rna-seq data
to illustrate the practical use of our software, we have re-analyzed  <dig> rna-seq samples from a study on transcription factors in drosophila where  <dig> mutant conditions were compared to a control with two biological replicates  <cit> .
after feeding the sample metadata into the quickngs database, we have linked the  <dig> fastq files into the quickngs stack directory.
these preparing steps took us an overall time of 2 min.
while waiting for the subsequent pipeline run to finish, we were able to monitor the current status of the respective modules using the status page on the quickngs database interface .
the rna-seq workflow comprises an initial quality check using fastqc plus some software which is unique to quickngs.
the basic data processing consists of a splicing-aware alignment using tophat <dig>  <cit>  followed by reference-guided transcriptome reassembly with cufflinks <dig>  <cit> .
differential gene expression and differential exon usage are analyzed with deseq <dig>  and dexseq .
after the processing finished overnight, we logged in to the quickngs user interface and found a report which summarizes all results of the quickngs workflow .
from the initial quality check, we received some basic read statistics  as well as standard qc plots, a heat map  and a plot from a principle component analysis  for the  <dig> samples.
the results of the core analysis for the comparison of atf <dig> mutants  against controls  are provided as additional files.
at thresholds  <dig> and  <dig>  for fold-change and p-value, we get a set of  <dig> differentially expressed genes  and a set of  <dig> differentially used exons .
additional file  <dig> reports the p-values and fold-changes for differential gene expression  together with those for the comparisons of the remaining three mutant conditions to control .
on the web interface, the same three spreadsheet files are given also for these comparisons.
all tables contain a comprehensive selection of genomic and functional annotation.
visualisation of the rna-seq wiggle files on the ucsc genome browser can be accessed by a hyperlink which uses a local password-protected track hub for the browser.
the fastq files for these test data are available from the ncbi short read archive  at accession number srp <dig> table  <dig> reads statistics on the test data from  <cit>  – read counts are given in multiples of  <dig> .
duplicate removal was not performed because this was a single-read analysis.
two samples  were treated with ribo-zero, whereas for the remaining samples, there is a significant degree of contamination with ribosomal rna.
for all samples, about the half of the reads map to the original strand because all data origin from unstranded libraries
fig.
<dig>
a heatmap on the  <dig> rna-seq test data sets: the replicates of each genotype do not perfectly cluster together in distinct subclusters.
this is likely to be cause by a combined effect of ribosomal contamination and batch effects.
b the principle component analysis confirms that two samples  which were processed in a separate batch and with ribozero treatment cluster distantly from the remaining samples
description of other quickngs workflows
although the current quickngs release also comprises workflows for mirna sequencing, chip-seq and whole-genome resequencing, we gave above a detailed description only for the rna-seq workflow.
however, the same level of efficiency and automation is also achieved in all other quickngs workflows.
the mirna-seq workflow comprises quantification and differential profiling of 3p and 5p mature mirnas using mirdeep  <cit>  as well as statistics on mirna families.
differential mirna expression is profiled with the deseq <dig> package  <cit> .
the chip-seq workflow takes advantage of bwa  <cit>  for genomic alignment of the reads and uses macs <dig>  <cit>  for peak calling.
furthermore, quickngs identifies all genes which are 2000 bp up- or downstream from the macs <dig> peaks.
the peak sequences are analyzed for enrichment of transcription factor binding motifs using meme-chip  <cit> .
the results comprise lists of significant peaks and reports for motif enrichment.
for the whole-genome resequencing workflow, finally, the software uses bwa for genomic alignment and calls single nucleotide polymorphisms with samtools  <cit>  and structural variations with delly  <cit> .
the results are annotated and functionally classified with snpeff  <cit> .
basic qc statistics and password-protected track hubs for the ucsc genome browser with direct hyperlinks for visualisation are part of all workflows.
the quickngs database comes with ready-made metadata for additional test data which are available from the sra at ncbi at accession numbers srp <dig> , srp <dig>  and srp <dig> .
additional modules dedicated to cancer genomics and more recent ngs applications such as clip-seq  are currently under development.
features of quickngs compared to other ngs workflow systems
in order to elaborate how quickngs performs in comparison to other ngs workflow systems, we discuss the features that are unique to our solution as well as its limitations.
the degree of automation in quickngs is much higher than, for instance, that of an appropriate workflow in popular data analysis frameworks like galaxy  <cit> , genepattern  <cit>  or chipster  <cit> .
this makes the system more efficient for the typical standard analyses, but also less flexible to modifications.
in particular, our system enables an extreme reduction of the hands-on  time that staff have to spend for the basic ngs data analyses.
data processing for tens or hundreds of samples can be initiated in less than 10 min.
while the subsequent analyses completely run in the background, they can be monitored on the status website and, once finished, the results are ready for immediate access by any scientist without specific it skills.
in contrast to all other systems, the quickngs database is capable of organizing sample meta information along with the analysis results, enabling a high degree of reproducibility and documentation of what analyses have been done.
this is essential whenever large numbers of samples are processed.
our software is also the only one to summarize all analysis results into user accounts with ready-to-deliver web reports.
an overview of the features of several ngs workflow systems compared to quickngs is given in table  <dig>
CONCLUSIONS
we have contributed quickngs, a software which extremely reduces the efforts put into basic data processing for next-generation sequencing.
quickngs takes advantage of powerful hardware architectures and a comprehensive database to control the overall workflow.
as a result, our approach enables laboratories with a high throughput of ngs data analyses to now accomplish their basic bioinformatics work on next-generation sequencing data essentially in zero time.
availability and requirements
project home page: http://bifacility.uni-koeln.de/quickngs/web
operating system: linux
programming languages: bash scripting, perl, r
other requirements: mysql server, apache web server on separate machine
license: gnu gpl version 3
additional files
additional file 1:
list of genes differentially expressed between the atf3a mutants and the yw controls in our test dataset.
the fold-changes and p-values are produced from the deseq <dig> package, whereas the fpkm values are taken from cufflinks output.
the file contains multiple tabs representing lists cut at particular thresholds for fold-change and p-value.
additional file 2:
list of exons differentially used between the atf3a mutants and the yw controls according to the dexseq package.
the file contains multiple tabs representing lists cut at particular thresholds for fold-change and p-value.
additional file 3:
list containing the same p-values and fold-changes as additional file
1
plus p-values and fold-changes for all other comparisons performed in the current quickngs runs.
in this example, these are the comparisons of the three other mutant condition against wild type.
this table facilitates comparisons between the genes differentially regulated under one condition and other conditions.
additional file 4:
procedure to upload sample meta data into the quickngs database.
first, the file locations of the raw data need to be saved into a text file.
together with information on library type, ngs application, species and laboratory, this file can be uploaded into the quickngs web interface.
human-readable sample labels as well as batch information can be provided for each sample listed in the text file.
pairs of samples or sample groups can be defined for comparative analysis within the workflow.
abbreviations
ngsnext-generation sequencing
dbdatabase
qcquality check
itinformation technology
sqlstructured query language
chipchromatin immunoprecipitation
competing interests
the authors declare that they have no competing interests.
authors’ contributions
pw, mn and pf developed the code.
pw tested and debugged the code in high-throughput production mode.
pf and mn wrote the manuscript.
all authors read and approved the final manuscript.

 next-generation sequencing has become the method of choice for molecular genetic analyses such as transcriptome profiling, chromatin immunoprecipitation followed by sequencing as well as resequencing of complete genomes. numerous software packages designed to analyze massive amounts of ngs data are now publicly available. preprocessing of ngs data typically takes advantage of a complex hardware architecture composed of, for instance, a parallel compute cluster, a database server and a web server. as this requires specialized it skills, the widespread access to ngs technology is still hampered by technical barriers. the primary data analysis is therefore often centralized into core laboratories which face the challenge of using a reasonable selection out of the available software packages to process a growing flow of newly produced data. we introduce a new framework named quickngs which can be operated by it staff in bioinformatics core labs to process vast amounts of data provided by their end users, typically experimental scientists. quickngs enables a major decrease in time and effort put into the primary analysis of ngs data, thus contributing to the further evolution of ngs into a standard tool with high accessibility to researchers. implementation quickngs enables rapid and professional data analysis for the aforementioned major applications of ngs in a batch-like operation mode. the core of the quickngs workflow is formed by a comprehensive mysql database used to organize sample metadata such as species, library type and batch information as well as the analysis results. the system operators can use the quickngs database interface to upload metadata and monitor the status of quickngs analysis runs from any location. on the other hand, the purpose of the quickngs web interface is to deliver web reports on the analysis results to the end users who can use their personal login and password to browse the results plus a lot of useful on-top information. the core results are provided as spreadsheet files which the user can download to a local computer.fig. a the quickngs database contains meta information on samples and sample groups. both can be efficiently organized by an intuitive web interface. new samples and sample groups can be inserted by following the „ + “button. b the status page on the quickngs database monitors time, user information and current status of each quickngs module on a clearly arranged website, enabling password-protected interrogation of the current status from any working location, including mobile access fig. results report for the test run on drosophila rna-seq samples. after login, the interface provides links to the main database export files, qc statistics as well as visualisation plots and the user’s personal password-protected ucsc genome browser track hub as the sample metadata in the quickngs database are used to completely control the overall workflow, these have to be provided to the quickngs database before any analysis can be started. to achieve this, the file locations of the raw data first need to be saved into a text file. this file is then uploaded into the quickngs web interface together with information on the library type, ngs application type, species of origin and the laboratory in which the input material has been generated. for each sample listed in the text file, the user is then asked to provide a human-readable sample label as well as batch information for the case that samples have been processed with different library preparation protocols or in different ngs runs. finally, samples can be assigned to pairs and groups for comparative analysis, e.g. differential gene expression. subsequently, the raw data received from the sequencing center are linked into the quickngs stack directory which is then processed fully automatically. thus, the only user action needed is providing sample information and linking the files to the stack directory which both are trivial amounts of effort. the results produced by the workflow are saved back into the quickngs database and made accessible through a report on the web interface. this report comprises standard qc metrics and results on typical data-related research questions, for instance which genes are differentially expressed or differentially spliced, which genomic variations are unique to a sample compared to a control, which transcription factor binding motifs are enriched in a chip-seq data set etc.. the results and report are generated fully automatically without any additional user action. the analysis relies on widely adopted ngs analysis packages which are listed in table  for the core analysis of the raw data, we have carefully selected the most appropriate previously published software programs. the selection criteria were performance in published and in-house benchmarking studies, comprehensiveness of the analysis output, quality of the implementation and steadiness of maintenance, and popularity in the community. our choice of bioinformatics software follows these criteria as far as possible. the code for qc and visualisation as well as for data management and the workflow itself is unique to quickngs. as a reference to genomic sequence and annotation, the system uses the mirbase for the mirna-seq workflow and the ensembl database for all other applications. for instance, rna-seq or chip-seq analyses can thus be carried out on data from any arbitrary organism listed on either ensembl or ensemblgenomes. the reference files are downloaded to the local system and updated automatically. the same applies to genomic annotation data which are retrieved using biomart.table algorithms and software tools used by quickngs, version.. the selection may be modified, updated or extended in future releases of quickngs, however, an up-to-date version of this table will be kept available online at the quickngs website  test run on previously published rna-seq data to illustrate the practical use of our software, we have re-analyzed rna-seq samples from a study on transcription factors in drosophila where mutant conditions were compared to a control with two biological replicates. after feeding the sample metadata into the quickngs database, we have linked the fastq files into the quickngs stack directory. these preparing steps took us an overall time of 2 min. while waiting for the subsequent pipeline run to finish, we were able to monitor the current status of the respective modules using the status page on the quickngs database interface. the rna-seq workflow comprises an initial quality check using fastqc plus some software which is unique to quickngs. the basic data processing consists of a splicing-aware alignment using tophat followed by reference-guided transcriptome reassembly with cufflinks. differential gene expression and differential exon usage are analyzed with deseq and dexseq. after the processing finished overnight, we logged in to the quickngs user interface and found a report which summarizes all results of the quickngs workflow. from the initial quality check, we received some basic read statistics as well as standard qc plots, a heat map and a plot from a principle component analysis for the samples. the results of the core analysis for the comparison of atf mutants against controls are provided as additional files. at thresholds and for fold-change and p-value, we get a set of differentially expressed genes and a set of differentially used exons. additional file reports the p-values and fold-changes for differential gene expression together with those for the comparisons of the remaining three mutant conditions to control. on the web interface, the same three spreadsheet files are given also for these comparisons. all tables contain a comprehensive selection of genomic and functional annotation. visualisation of the rna-seq wiggle files on the ucsc genome browser can be accessed by a hyperlink which uses a local password-protected track hub for the browser. the fastq files for these test data are available from the ncbi short read archive at accession number srp table reads statistics on the test data from – read counts are given in multiples of. duplicate removal was not performed because this was a single-read analysis. two samples were treated with ribo-zero, whereas for the remaining samples, there is a significant degree of contamination with ribosomal rna. for all samples, about the half of the reads map to the original strand because all data origin from unstranded libraries fig. a heatmap on the rna-seq test data sets: the replicates of each genotype do not perfectly cluster together in distinct subclusters. this is likely to be cause by a combined effect of ribosomal contamination and batch effects. b the principle component analysis confirms that two samples which were processed in a separate batch and with ribozero treatment cluster distantly from the remaining samples description of other quickngs workflows although the current quickngs release also comprises workflows for mirna sequencing, chip-seq and whole-genome resequencing, we gave above a detailed description only for the rna-seq workflow. however, the same level of efficiency and automation is also achieved in all other quickngs workflows. the mirna-seq workflow comprises quantification and differential profiling of 3p and 5p mature mirnas using mirdeep as well as statistics on mirna families. differential mirna expression is profiled with the deseq package. the chip-seq workflow takes advantage of bwa for genomic alignment of the reads and uses macs for peak calling. furthermore, quickngs identifies all genes which are 2000 bp up- or downstream from the macs peaks. the peak sequences are analyzed for enrichment of transcription factor binding motifs using meme-chip. the results comprise lists of significant peaks and reports for motif enrichment. for the whole-genome resequencing workflow, finally, the software uses bwa for genomic alignment and calls single nucleotide polymorphisms with samtools and structural variations with delly. the results are annotated and functionally classified with snpeff. basic qc statistics and password-protected track hubs for the ucsc genome browser with direct hyperlinks for visualisation are part of all workflows. the quickngs database comes with ready-made metadata for additional test data which are available from the sra at ncbi at accession numbers srp, srp and srp. additional modules dedicated to cancer genomics and more recent ngs applications such as clip-seq are currently under development. features of quickngs compared to other ngs workflow systems in order to elaborate how quickngs performs in comparison to other ngs workflow systems, we discuss the features that are unique to our solution as well as its limitations. the degree of automation in quickngs is much higher than, for instance, that of an appropriate workflow in popular data analysis frameworks like galaxy , genepattern  or chipster . this makes the system more efficient for the typical standard analyses, but also less flexible to modifications. in particular, our system enables an extreme reduction of the hands-on time that staff have to spend for the basic ngs data analyses. data processing for tens or hundreds of samples can be initiated in less than 10 min. while the subsequent analyses completely run in the background, they can be monitored on the status website and, once finished, the results are ready for immediate access by any scientist without specific it skills. in contrast to all other systems, the quickngs database is capable of organizing sample meta information along with the analysis results, enabling a high degree of reproducibility and documentation of what analyses have been done. this is essential whenever large numbers of samples are processed. our software is also the only one to summarize all analysis results into user accounts with ready-to-deliver web reports. an overview of the features of several ngs workflow systems compared to quickngs is given in table  we have contributed quickngs, a software which extremely reduces the efforts put into basic data processing for next-generation sequencing. quickngs takes advantage of powerful hardware architectures and a comprehensive database to control the overall workflow. as a result, our approach enables laboratories with a high throughput of ngs data analyses to now accomplish their basic bioinformatics work on next-generation sequencing data essentially in zero time. availability and requirements project home page: http://bifacility.uni-koeln.de/quickngs/web operating system: linux programming languages: bash scripting, perl, r other requirements: mysql server, apache web server on separate machine license: gnu gpl version 3 additional files additional file 1: list of genes differentially expressed between the atf3a mutants and the yw controls in our test dataset. the fold-changes and p-values are produced from the deseq package, whereas the fpkm values are taken from cufflinks output. the file contains multiple tabs representing lists cut at particular thresholds for fold-change and p-value. additional file 2: list of exons differentially used between the atf3a mutants and the yw controls according to the dexseq package. the file contains multiple tabs representing lists cut at particular thresholds for fold-change and p-value. additional file 3: list containing the same p-values and fold-changes as additional file 1 plus p-values and fold-changes for all other comparisons performed in the current quickngs runs. in this example, these are the comparisons of the three other mutant condition against wild type. this table facilitates comparisons between the genes differentially regulated under one condition and other conditions. additional file 4: procedure to upload sample meta data into the quickngs database. first, the file locations of the raw data need to be saved into a text file. together with information on library type, ngs application, species and laboratory, this file can be uploaded into the quickngs web interface. human-readable sample labels as well as batch information can be provided for each sample listed in the text file. pairs of samples or sample groups can be defined for comparative analysis within the workflow. abbreviations ngsnext-generation sequencing dbdatabase qcquality check itinformation technology sqlstructured query language chipchromatin immunoprecipitation competing interests the authors declare that they have no competing interests. authors’ contributions pw, mn and pf developed the code. pw tested and debugged the code in high-throughput production mode. pf and mn wrote the manuscript. all authors read and approved the final manuscript. 
we demonstrate the efficiency of our new software by a comprehensive analysis of rna-seq samples which we can finish in only a few minutes of hands-on time. since multiple algorithmic approaches for the basic data analysis have been developed, there is now an increasing need to efficiently use these tools to obtain results in reasonable time. while time and cost for the sequencing itself are decreasing, the analysis of the massive amounts of data remains challenging.
 next-generation sequencing has emerged as a widely used tool in molecular biology. while time and cost for the sequencing itself are decreasing, the analysis of the massive amounts of data remains challenging. since multiple algorithmic approaches for the basic data analysis have been developed, there is now an increasing need to efficiently use these tools to obtain results in reasonable time.  we have developed quickngs, a new workflow system for laboratories with the need to analyze data from multiple ngs projects at a time. quickngs takes advantage of parallel computing resources, a comprehensive back-end database, and a careful selection of previously published algorithmic approaches to build fully automated data analysis workflows. we demonstrate the efficiency of our new software by a comprehensive analysis of rna-seq samples which we can finish in only a few minutes of hands-on time. the approach we have taken is suitable to process even much larger numbers of samples and multiple projects at a time.  our approach considerably reduces the barriers that still limit the usability of the powerful ngs technology and finally decreases the time to be spent before proceeding to further downstream analysis and interpretation of the data. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords next-generation sequencingbatch processingdata managementhigh-performance computinganalysis workflowissue-copyright-statement© the author 2015 