NAME

Bio::ViennaNGS - A Perl distribution for Next-Generation Sequencing (NGS) data analysis

DESCRIPTION

Bio::ViennaNGS is a distribution of Perl modules and utilities for building efficient Next-Generation Sequencing (NGS) analysis pipelines. It covers various aspects of NGS data analysis, including (but not limited to) conversion of sequence annotation, evaluation of mapped data, expression quantification and visualization.

The main Bio::ViennaNGS module is shipped with a complementary set of (sub)modules:

Bio::ViennaNGS::AnnoC: A Moose interface for storage and conversion of sequence annotation data.
Bio::ViennaNGS::Bam: Routines for high-level manipulation of BAM files.
Bio::ViennaNGS::BamStat: A Moose based class for collecting mapping statistics.
Bio::ViennaNGS::BamStatSummary: A Moose interface for processing Bio::ViennaNGS::BamStat objects on a set of BAM files.
Bio::ViennaNGS::Bed: A Moose interface for manipulation of genomic interval data in BED format.
Bio::ViennaNGS::BedGraphEntry: A Moose interface for storing genomic cvoerage and interval data in bedGraph format.
Bio::ViennaNGS::Expression: An object oriented interface for read-count based gene expression analysis.
Bio::ViennaNGS::ExtFeature: A Moose wrapper for extended BED6 entries.
Bio::ViennaNGS::Fasta: Routines for accessing genomic sequences implemented through a Moose interface to Bio::DB::Fasta.
Bio::ViennaNGS::Feature: A Moose based BED6 wrapper.
Bio::ViennaNGS::FeatureChain: Yet another Moose class for chaining gene annotation features.
Bio::ViennaNGS::FeatureInterval: A Moose interface for handling elementary genomic intervals, corresponding to BED3.
Bio::ViennaNGS::FeatureIO: A Moose interface for efficient input/output handling of genomic annotation formats.
Bio::ViennaNGS::FeatureLine: An abstract Moose class for combining several Bio::ViennaNGS::FeatureChain objects.
Bio::ViennaNGS::MinimalFeature: A Moose interface for handling elementary gene annotation, corresponding to BED4.
Bio::ViennaNGS::Peak: A Moose interface for identification and characterization of peaks/enriched regions in RNA-seq data.
Bio::ViennaNGS::SpliceJunc: A collection of routines for alternative splicing analysis.
Bio::ViennaNGS::Tutorial: A comprehensive tutorial of the Bio::ViennaNGS core routines with real-world NGS data.
Bio::ViennaNGS::UCSC: Routines for visualization of genomics data with the UCSC genome browser.
Bio::ViennaNGS::Util: A collection of wrapper routines for commonly used third-party NGS utilities, code for normalization of gene expression values based on read count data and a set of utility functions.

UTILITIES

Bio::ViennaNGS comes with a collection of command line utilities for accomplishing routine tasks often required in NGS data processing. These utilities serve as reference implementation of the routines implemented throughout the modules and can readily be used for atomic tasks in NGS data processing:

assembly_hub_constructor.pl: The UCSC genome browser offers the possibility to visualize any organism (including organisms that are not included in the standard UCSC browser bundle) through hso called 'Assembly Hubs'. This script constructs Assembly Hubs from genomic sequence and annotation data.
bam_split.pl: Split (paired-end and single-end) BAM alignment files by strand and compute statistics. Optionally create BED output, as well as normalized bedGraph and bigWig files for coverage visualization in genome browsers (see dependencies on third-patry tools below).
bam_to_bigWig.pl: Produce bigWig coverage profiles from (aligned) BAM files, explicitly considering strandedness. The most natural use case of this tool is to create strand-aware coverage profiles in bigWig format for genome browser visualization.
bam_uniq.pl: Extract unique and multi mapping reads from BAM alignment files and create a separate BAM file for both uniqe (.uniq.) and multi (.mult.) mappers.
bed2bedGraph.pl: Convert BED files to (strand specific) bedGraph files, allowing additional annotation and automatic generation of bedGraph files which can easily be converted to big-type files for easy UCSC visualization.
bed2nt2aa.pl: Provide nucleotide and amino acid sequences for BED6 intervals.
extend_bed.pl: Extend genomic features in BED files by a certain number of nucleotides, either on both sides or specifically at the 5' or 3' end, respectively.
gff2bed.pl: Convert RefSeq GFF3 annotation files to BED12 format. Individual BED12 files are created for each feature type (CDS/tRNA/rRNA/etc.). Tested with RefSeq bacterial GFF3 annotation.
kmer_analysis.pl: Count k-mers of predefined length in FastQ and Fasta files
MEME_XML_motif_extractor.pl: Compute simple statistics from MEME XML output and return a list of found motifs with the number of sequences containing those motifs as well as nice ggplot graphs.
newUCSCdb.pl: Create a new genome database to a locally installed instance of the UCSC genome browser in order to add a novel organism for visualization. Based on this Genomewiki article.
normalize_multicov.pl: Compute normalized expression data in RPKM and TPM from (raw) read counts in bedtools multicov format. TPM reference: Wagner et al, Theory Biosci. 131(4), pp 281-85 (2012).
rnaseq_peakfinder.pl: Find and characterize peaks/enriched regions of certain size and coverage in RNA-seq data.
sj_visualizer.pl: Convert splice junctions from mapped RNA-seq data in segemehl BED6 splice junction format to BED12 for easy visualization in genome Browsers.
splice_site_summary.pl: Identify and characterize splice junctions from RNA-seq data by intersecting them with annotated splice junctions.
trim_fastq.pl: Trim sequence and quality string fields in a Fastq file by user defined length.

DEPENDENCIES

The Bio::ViennaNGS modules and classes depend on a set of Perl modules, some of which are part of the Perl core distribution:

Bio::Perl >= 1.00690001
Bio::DB::Sam >= 1.37
Bio::DB::Fasta
Bio::Tools::GFF
File::Basename
File::Share
File::Slurp
File::Temp
List::Util
Path::Class
IPC::Cmd
Carp
Template
Moose
Moose::Util::TypeConstraints
namespace::autoclean
MooseX::Clone
MooseX::InstanceTracking
Tie::Hash::Indexed
Test::Files
Test::File::Contents

In addition the following modules are required by the Bio::ViennaNGS utilities:

PerlIO::gzip
Math::Round
XML::Simple
Statistics::R

Bio::ViennaNGS depends on a set of third-party tools and libraries which are required for specific filtering and file format conversion tasks as well as for building internally used Perl modules:

bedtools2
bedGraphToBigWig, fetchChromSizes, faToTwoBit from the UCSC Genome Browser applications
R
samtools <=v0.1.19 for building Bio::DB::Sam.

Please ensure that all third-party utilities are available on your system, and that hey can be found and executed by the Perl interpreter.

SOURCE AVAILABILITY

Source code for this distribution is available from the ViennaNGS Github repository.

PAPERS

If the Bio::ViennaNGS suite is useful for your work or if you use any component of the distribution in a custom pipeline, please cite the following publication:

"ViennaNGS - A toolbox for building efficient next-generation sequencing analysis pipelines"

Michael T. Wolfinger, Joerg Fallmann, Florian Eggenhofer and Fabian Amman

F1000Research 2015, 4:50 (doi: 10.12688/f1000research.6157.2)

NOTES

The Bio::ViennaNGS suite is actively developed and tested on different flavours of Linux and Mac OS X. We have taken care of writing platform-independent code that should run out of the box on most UNIX-based systems, however we do not have access to machines running Microsoft Windows. As such, we have not tested and will not test Windows compatibility.

AUTHORS

Michael T. Wolfinger <michael@wolfinger.eu>
Joerg Fallmann <fall@tbi.univie.ac.at>
Florian Eggenhofer <florian.eggenhofer@tbi.univie.ac.at>
Fabian Amman <fabian@tbi.univie.ac.at<gt>

COPYRIGHT AND LICENSE

This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself, either Perl version 5.10.0 or, at your option, any later version of Perl 5 you may have available.

This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

To install Bio::ViennaNGS, copy and paste the appropriate command in to your terminal.

cpanm

cpanm Bio::ViennaNGS

CPAN shell

perl -MCPAN -e shell
install Bio::ViennaNGS

For more information on module installation, please visit the detailed CPAN module installation guide.

	Global
`s`	Focus search bar
`?`	Bring up this help dialog

	GitHub
`g` `p`	Go to pull requests
`g` `i`	go to github issues (only if github is preferred repository)

	POD
`g` `a`	Go to author
`g` `c`	Go to changes
`g` `i`	Go to issues
`g` `d`	Go to dist
`g` `r`	Go to repository/SCM
`g` `s`	Go to source
`g` `b`	Go to file browse

	Search terms
module: (e.g. module:Plugin)
distribution: (e.g. distribution:Dancer auth)
author: (e.g. author:SONGMU Redis)
version: (e.g. version:1.00)