![]() | Name | Last modified | Size | Description |
---|---|---|---|---|
![]() | Parent Directory | - | ||
![]() | stamp-h1 | 2019-05-16 16:21 | 23 | |
![]() | sparsehash/ | 2018-06-12 18:09 | - | |
![]() | m4/ | 2019-05-03 11:31 | - | |
![]() | lib/ | 2019-05-03 11:31 | - | |
![]() | kmerprint/ | 2019-05-16 16:25 | - | |
![]() | doxygen.conf | 2019-05-03 11:31 | 50K | |
![]() | doc/ | 2019-05-16 16:21 | - | |
![]() | dialign/ | 2019-05-16 16:21 | - | |
![]() | configure.ac | 2019-05-03 11:31 | 9.5K | |
![]() | configure | 2019-05-16 16:15 | 374K | |
![]() | config.status | 2019-05-16 16:21 | 39K | |
![]() | config.log | 2019-05-16 16:21 | 110K | |
![]() | config.h.in | 2019-05-16 16:15 | 9.0K | |
![]() | config.h | 2019-05-16 16:21 | 9.5K | |
![]() | bin/ | 2019-05-16 16:21 | - | |
![]() | autom4te.cache/ | 2019-05-16 16:15 | - | |
![]() | autogen.sh | 2019-05-03 11:31 | 58 | |
![]() | aclocal.m4 | 2019-05-16 16:15 | 34K | |
![]() | _config.yml | 2019-05-03 11:31 | 26 | |
![]() | Unittest/ | 2019-05-16 16:21 | - | |
![]() | SimpleGraph/ | 2019-05-16 16:25 | - | |
![]() | Sealer/ | 2019-05-16 16:25 | - | |
![]() | Scaffold/ | 2019-05-16 16:25 | - | |
![]() | PopBubbles/ | 2019-05-16 16:24 | - | |
![]() | PathOverlap/ | 2019-05-16 16:24 | - | |
![]() | ParseAligns/ | 2019-05-16 16:24 | - | |
![]() | Parallel/ | 2019-05-16 16:21 | - | |
![]() | PairedDBG/ | 2019-05-16 16:24 | - | |
![]() | Overlap/ | 2019-05-16 16:24 | - | |
![]() | Misc/ | 2019-05-16 16:21 | - | |
![]() | MergePaths/ | 2019-05-16 16:24 | - | |
![]() | Map/ | 2019-05-16 16:23 | - | |
![]() | Makefile.in | 2019-05-16 16:15 | 25K | |
![]() | Makefile.am | 2019-05-03 11:31 | 855 | |
![]() | Makefile | 2019-05-16 16:21 | 25K | |
![]() | LogKmerCount/ | 2019-05-16 16:23 | - | |
![]() | Layout/ | 2019-05-16 16:23 | - | |
![]() | LICENSE | 2019-05-03 11:31 | 33K | |
![]() | Konnector/ | 2019-05-16 16:22 | - | |
![]() | KAligner/ | 2019-05-16 16:23 | - | |
![]() | IntegrationTest/ | 2019-05-03 11:31 | - | |
![]() | Graph/ | 2019-05-16 16:21 | - | |
![]() | GapFiller/ | 2019-05-16 16:25 | - | |
![]() | FilterGraph/ | 2019-05-16 16:25 | - | |
![]() | FMIndex/ | 2019-05-16 16:21 | - | |
![]() | Dockerfile | 2019-05-03 11:31 | 786 | |
![]() | DistanceEst/ | 2019-05-16 16:22 | - | |
![]() | DataLayer/ | 2019-05-16 16:21 | - | |
![]() | DataBase/ | 2019-05-16 16:21 | - | |
![]() | DAssembler/ | 2019-05-16 16:22 | - | |
![]() | Consensus/ | 2019-05-16 16:22 | - | |
![]() | Common/ | 2019-05-16 16:21 | - | |
![]() | ChangeLog | 2019-05-03 11:31 | 48K | |
![]() | COPYRIGHT | 2019-05-03 11:31 | 6.3K | |
![]() | CITATION.md | 2019-05-03 11:31 | 1.8K | |
![]() | CITATION.bib | 2019-05-03 11:31 | 556 | |
![]() | BloomDBG/ | 2019-05-16 16:22 | - | |
![]() | Bloom/ | 2019-05-16 16:22 | - | |
![]() | BIRCH-README.txt | 2019-05-16 16:18 | 740 | |
![]() | Assembly/ | 2019-05-16 16:21 | - | |
![]() | Align/ | 2019-05-16 16:21 | - | |
![]() | AdjList/ | 2019-05-16 16:26 | - | |
![]() | ABYSS/ | 2019-05-16 16:22 | - | |
ABySS is a de novo sequence assembler intended for short paired-end reads and large genomes.
Please cite our papers.
3 May 2019
Looking for a fun & worthy challenge? Think you can contribute code to this project? Join our team of developers! We are currently looking for C++ bioinformatics programmers. Inquire for staff, graduate student, and postdoctoral positions.
Install Linuxbrew, and run the command
brew install abyss
Install Homebrew, and run the command
brew install abyss
Install Windows Subsystem for Linux and Linuxbrew, and run the command
brew install abyss
Run the command
sudo apt-get install abyss
wget http://www.bcgsc.ca/platform/bioinfo/software/abyss/releases/1.3.4/test-data.tar.gz
tar xzvf test-data.tar.gz
abyss-pe k=25 name=test \
in='test-data/reads1.fastq test-data/reads2.fastq'
abyss-fac test-unitigs.fa
Dependencies may be installed using the package manager Homebrew on macOS and Linxubrew on Linux and Windows, using Windows Subsystem for Linux.
ABySS requires a C++ compiler that supports OpenMP such as GCC.
ABySS requires the following libraries:
ABySS will receive an error when compiling with Boost 1.51.0 or 1.52.0 since they contain a bug. Later versions of Boost compile without error.
Tigmint to correct assembly errors
brew install brewsci/bio/arcs brewsci/bio/links-scaffolder
zsh for reporting time and memory usage
brew install pigz samtools zsh
When installing ABySS from GitHub source the following tools are required:
To generate the configure script and make files:
./autogen.sh
See "Compiling ABySS from source" for further steps.
To compile and install ABySS in /usr/local
:
./configure
make
sudo make install
To install ABySS in a specified directory:
./configure --prefix=/opt/abyss
make
sudo make install
ABySS uses OpenMP for parallelization, which requires a modern compiler such as GCC 4.2 or greater. If you have an older compiler, it is best to upgrade your compiler if possible. If you have multiple versions of GCC installed, you can specify a different compiler:
./configure CC=gcc-4.6 CXX=g++-4.6
ABySS requires the Boost C++ libraries. Many systems come with Boost installed. If yours does not, you can download Boost. It is not necessary to compile Boost before installing it. The Boost header file directory should be found at /usr/include/boost
, in the ABySS source directory, or its location specified to configure
:
./configure --with-boost=/usr/local/include
If you wish to build the parallel assembler with MPI support, MPI should be found in /usr/include
and /usr/lib
or its location specified to configure
:
./configure --with-mpi=/usr/lib/openmpi
While OpenMPI is assumed by default you can switch to LAM/MPI or MPICH using: ./configure --enable-mpich use MPICH (default is to use Open MPI) ./configure --enable-lammpi use LAM/MPI (default is to use Open MPI)
ABySS should be built using the sparsehash library to reduce memory usage, although it will build without. sparsehash should be found in /usr/include
or its location specified to configure
:
./configure CPPFLAGS=-I/usr/local/include
If the optional dependency SQLite is installed in non-default directories, its location can be specified to configure
:
./configure --with-sqlite=/opt/sqlite3
The default maximum k-mer size is 128 and may be decreased to reduce memory usage or increased at compile time. This value must be a multiple of 32 (i.e. 32, 64, 96, 128, etc):
./configure --enable-maxk=160
If you encounter compiler warnings, you may ignore them like so:
make AM_CXXFLAGS=-Wall
To run ABySS, its executables should be found in your PATH
. If you installed ABySS in /opt/abyss
, add /opt/abyss/bin
to your PATH
:
PATH=/opt/abyss/bin:$PATH
ABySS stores temporary files in TMPDIR
, which is /tmp
by default on most systems. If your default temporary disk volume is too small, set TMPDIR
to a larger volume, such as /var/tmp
or your home directory.
export TMPDIR=/var/tmp
To assemble paired reads in two files named reads1.fa
and reads2.fa
into contigs in a file named ecoli-contigs.fa
, run the command:
abyss-pe name=ecoli k=96 in='reads1.fa reads2.fa'
The parameter in
specifies the input files to read, which may be in FASTA, FASTQ, qseq, export, SRA, SAM or BAM format and compressed with gz, bz2 or xz and may be tarred. The assembled contigs will be stored in ${name}-contigs.fa
.
A pair of reads must be named with the suffixes /1
and /2
to identify the first and second read, or the reads may be named identically. The paired reads may be in separate files or interleaved in a single file.
Reads without mates should be placed in a file specified by the parameter se
(single-end). Reads without mates in the paired-end files will slow down the paired-end assembler considerably during the abyss-fixmate
stage.
The distribution of fragment sizes of each library is calculated empirically by aligning paired reads to the contigs produced by the single-end assembler, and the distribution is stored in a file with the extension .hist
, such as ecoli-3.hist
. The N50 of the single-end assembly must be well over the fragment-size to obtain an accurate empirical distribution.
Here's an example scenario of assembling a data set with two different fragment libraries and single-end reads. Note that the names of the libraries (pea
and peb
) are arbitrary.
pea
has reads in two files, pea_1.fa
and pea_2.fa
.peb
has reads in two files, peb_1.fa
and peb_2.fa
.se1.fa
and se2.fa
.The command line to assemble this example data set is:
abyss-pe k=96 name=ecoli lib='pea peb' \
pea='pea_1.fa pea_2.fa' peb='peb_1.fa peb_2.fa' \
se='se1.fa se2.fa'
The empirical distribution of fragment sizes will be stored in two files named pea-3.hist
and peb-3.hist
. These files may be plotted to check that the empirical distribution agrees with the expected distribution. The assembled contigs will be stored in ${name}-contigs.fa
.
Long-distance mate-pair libraries may be used to scaffold an assembly. Specify the names of the mate-pair libraries using the parameter mp
. The scaffolds will be stored in the file ${name}-scaffolds.fa
. Here's an example of assembling a data set with two paired-end libraries and two mate-pair libraries. Note that the names of the libraries (pea
, peb
, mpa
, mpb
) are arbitrary.
abyss-pe k=96 name=ecoli lib='pea peb' mp='mpc mpd' \
pea='pea_1.fa pea_2.fa' peb='peb_1.fa peb_2.fa' \
mpc='mpc_1.fa mpc_2.fa' mpd='mpd_1.fa mpd_2.fa'
The mate-pair libraries are used only for scaffolding and do not contribute towards the consensus sequence.
ABySS can scaffold using linked reads from 10x Genomics Chromium. The barcodes must first be extracted from the read sequences and added to the BX:Z
tag of the FASTQ header, typically using the longranger basic
command of Long Ranger or EMA preproc. The linked reads are used to correct assembly errors, which requires that Tigmint. The linked reads are also used for scaffolding, which requires ARCS. See Dependencies for installation instructions.
ABySS can combine paired-end, mate-pair, and linked-read libraries. The pe
and lr
libraries will be used to build the de Bruijn graph. The mp
libraries will be used for paired-end/mate-pair scaffolding. The lr
libraries will be used for misassembly correction using Tigmint and scaffolding using ARCS.
abyss-pe k=96 name=hsapiens \
pe='pea' pea='lra.fastq.gz' \
mp='mpa' mpa='lra.fastq.gz' \
lr='lra' lra='lra.fastq.gz'
ABySS performs better with a mixture of paired-end, mate-pair, and linked reads, but it is possible to assemble only linked reads using ABySS, though this mode of operation is experimental.
abyss-pe k=96 name=hsapiens lr='lra' lra='lra.fastq.gz'
Long sequences such as RNA-Seq contigs can be used to rescaffold an assembly. Sequences are aligned using BWA-MEM to the assembled scaffolds. Additional scaffolds are then formed between scaffolds that can be linked unambiguously when considering all BWA-MEM alignments.
Similar to scaffolding, the names of the datasets can be specified with the long
parameter. These scaffolds will be stored in the file ${name}-long-scaffs.fa
. The following is an example of an assembly with PET, MPET and an RNA-Seq assembly. Note that the names of the libraries are arbitrary.
abyss-pe k=96 name=ecoli lib='pe1 pe2' mp='mp1 mp2' long='longa' \
pe1='pe1_1.fa pe1_2.fa' pe2='pe2_1.fa pe2_2.fa' \
mp1='mp1_1.fa mp1_2.fa' mp2='mp2_1.fa mp2_2.fa' \
longa='longa.fa'
Assemblies may be performed using a Bloom filter de Bruijn graph, which typically reduces memory requirements by an order of magnitude. To assemble in Bloom filter mode, the user must specify 3 additional parameters: B
(Bloom filter size in bytes), H
(number of Bloom filter hash functions), and kc
(minimum k-mer count threshold). B
is the overall memory budget for the Bloom filter assembler, and may be specified with unit suffixes 'k' (kilobytes), 'M' (megabytes), 'G' (gigabytes). If no units are specified bytes are assumed. For example, the following will run a E. coli assembly with an overall memory budget of 100 megabytes, 3 hash functions, a minimum k-mer count threshold of 3, with verbose logging enabled:
abyss-pe name=ecoli k=96 in='reads1.fa reads2.fa' B=100M H=3 kc=3 v=-v
At the current time, the user must calculate suitable values for B
and H
on their own, and finding the best value for kc
may require experimentation (optimal values are typically in the range of 2-4). Internally, the Bloom filter assembler divides the memory budget (B
) equally across (kc
+ 1) Bloom filters, where kc
Bloom filters are used for the cascading Bloom filter and one additional Bloom filter is used to track k-mers that have previously been included in contigs. Users are recommended to target a Bloom filter false positive rate (FPR) that is less than 5%, as reported by the assembly log when using the v=-v
option (verbose level 1).
Assemblies may be performed using a paired de Bruijn graph instead of a standard de Bruijn graph. In paired de Bruijn graph mode, ABySS uses k-mer pairs in place of k-mers, where each k-mer pair consists of two equal-size k-mers separated by a fixed distance. A k-mer pair is functionally similar to a large k-mer spanning the breadth of the k-mer pair, but uses less memory because the sequence in the gap is not stored. To assemble using paired de Bruijn graph mode, specify both individual k-mer size (K
) and k-mer pair span (k
). For example, to assemble E. coli with a individual k-mer size of 16 and a k-mer pair span of 96:
abyss-pe name=ecoli K=16 k=96 in='reads1.fa reads2.fa'
In this example, the size of the intervening gap between k-mer pairs is 64 bp (96 - 2*16). Note that the k
parameter takes on a new meaning in paired de Bruijn graph mode. k
indicates kmer pair span in paired de Bruijn graph mode (when K
is set), whereas k
indicates k-mer size in standard de Bruijn graph mode (when K
is not set).
Strand-specific RNA-Seq libraries can be assembled such that the resulting unitigs, contigs and scaffolds are oriented correctly with respect to the original transcripts that were sequenced. In order to run ABySS in strand-specific mode, the SS
parameter must be used as in the following example:
abyss-pe name=SS-RNA k=96 in='reads1.fa reads2.fa' SS=--SS
The expected orientation for the read sequences with respect to the original RNA is RF. i.e. the first read in a read pair is always in reverse orientation.
To find the optimal value of k
, run multiple assemblies and inspect the assembly contiguity statistics. The following shell snippet will assemble for every eighth value of k
from 50 to 90.
for k in `seq 50 8 90`; do
mkdir k$k
abyss-pe -C k$k name=ecoli k=$k in=../reads.fa
done
abyss-fac k*/ecoli-contigs.fa
The default maximum value for k
is 96. This limit may be changed at compile time using the --enable-maxk
option of configure. It may be decreased to 32 to decrease memory usage or increased to larger values.
The np
option of abyss-pe
specifies the number of processes to use for the parallel MPI job. Without any MPI configuration, this will allow you to use multiple cores on a single machine. To use multiple machines for assembly, you must create a hostfile
for mpirun
, which is described in the mpirun
man page.
Do not run mpirun -np 8 abyss-pe
. To run ABySS with 8 threads, use abyss-pe np=8
. The abyss-pe
driver script will start the MPI process, like so: mpirun -np 8 ABYSS-P
.
The paired-end assembly stage is multithreaded, but must run on a single machine. The number of threads to use may be specified with the parameter j
. The default value for j
is the value of np
.
ABySS integrates well with cluster job schedulers, such as:
For example, to submit an array of jobs to assemble every eighth value of k
between 50 and 90 using 64 processes for each job:
qsub -N ecoli -pe openmpi 64 -t 50-90:8 \
<<<'mkdir k$SGE_TASK_ID && abyss-pe -C k$SGE_TASK_ID in=/data/reads.fa'
ABySS supports the use of DIDA (Distributed Indexing Dispatched Alignment), an MPI-based framework for computing sequence alignments in parallel across multiple machines. The DIDA software must be separately downloaded and installed from http://www.bcgsc.ca/platform/bioinfo/software/dida. In comparison to the standard ABySS alignment stages which are constrained to a single machine, DIDA offers improved performance and the ability to scale to larger targets. Please see the DIDA section of the abyss-pe man page (in the doc
subdirectory) for details on usage.
Parameters of the driver script, abyss-pe
a
: maximum number of branches of a bubble [2
]b
: maximum length of a bubble (bp) [""
]B
: Bloom filter size (e.g. "100M")c
: minimum mean k-mer coverage of a unitig [sqrt(median)
]d
: allowable error of a distance estimate (bp) [6
]e
: minimum erosion k-mer coverage [round(sqrt(median))
]E
: minimum erosion k-mer coverage per strand [1 if sqrt(median) > 2
else 0]G
: genome size, used to calculate NG50H
: number of Bloom filter hash functions [1
]j
: number of threads [2
]k
: size of k-mer (when K
is not set) or the span of a k-mer pair (when K
is set)kc
: minimum k-mer count threshold for Bloom filter assembly [2
]K
: the length of a single k-mer in a k-mer pair (bp)l
: minimum alignment length of a read (bp) [40
]m
: minimum overlap of two unitigs (bp) [k-1
]n
: minimum number of pairs required for building contigs [10
]N
: minimum number of pairs required for building scaffolds [n
]np
: number of MPI processes [1
]p
: minimum sequence identity of a bubble [0.9
]q
: minimum base quality [3
]s
: minimum unitig size required for building contigs (bp) [1000
]S
: minimum contig size required for building scaffolds (bp) [1000-10000
]t
: maximum length of blunt contigs to trim [k
]v
: use v=-v
for verbose logging, v=-vv
for extra verbosex
: spaced seed (Bloom filter assembly only)lr_s
: minimum contig size required for building scaffolds with linked reads (bp) [S
]lr_n
: minimum number of barcodes required for building scaffolds with linked reads [10
]Please see the abyss-pe manual page for more information on assembly parameters.
abyss-pe
configuration variables may be set on the command line or from the environment, for example with export k=96
. It can happen that abyss-pe
picks up such variables from your environment that you had not intended, and that can cause trouble. To troubleshoot that situation, use the abyss-pe env
command to print the values of all the abyss-pe
configuration variables:
abyss-pe env [options]
abyss-pe
is a driver script implemented as a Makefile. Any option of make
may be used with abyss-pe
. Particularly useful options are:
-C dir
, --directory=dir
Change to the directory dir
and store the results there.-n
, --dry-run
Print the commands that would be executed, but do not execute them.abyss-pe
uses the following programs, which must be found in your PATH
:
ABYSS
: de Bruijn graph assemblerABYSS-P
: parallel (MPI) de Bruijn graph assemblerAdjList
: find overlapping sequencesDistanceEst
: estimate the distance between sequencesMergeContigs
: merge sequencesMergePaths
: merge overlapping pathsOverlap
: find overlapping sequences using paired-end readsPathConsensus
: find a consensus sequence of ambiguous pathsPathOverlap
: find overlapping pathsPopBubbles
: remove bubbles from the sequence overlap graphSimpleGraph
: find paths through the overlap graphabyss-fac
: calculate assembly contiguity statisticsabyss-filtergraph
: remove shim contigs from the overlap graphabyss-fixmate
: fill the paired-end fields of SAM alignmentsabyss-map
: map reads to a reference sequenceabyss-scaffold
: scaffold contigs using distance estimatesabyss-todot
: convert graph formats and merge graphsThis flowchart shows the ABySS assembly pipeline its intermediate files.
ABySS has a built-in support for SQLite database to export log values into a SQLite file and/or .csv
files at runtime.
Of abyss-pe
: * db
: path to SQLite repository file [$(name).sqlite
] * species
: name of species to archive [ ] * strain
: name of strain to archive [ ] * library
: name of library to archive [ ]
For example, to export data of species 'Ecoli', strain 'O121' and library 'pea' into your SQLite database repository named '/abyss/test.sqlite':
abyss-pe db=/abyss/test.sqlite species=Ecoli strain=O121 library=pea [other options]
Found in your path
:
abyss-db-txt
: create a flat file showing entire repository at a glanceabyss-db-csv
: create .csv
table(s) from the repositoryUsage:
abyss-db-txt /your/repository
abyss-db-csv /your/repository program(s)
For example,
abyss-db-txt repo.sqlite
abyss-db-csv repo.sqlite DistanceEst
abyss-db-csv repo.sqlite DistanceEst abyss-scaffold
abyss-db-csv repo.sqlite --all
Shaun D Jackman, Benjamin P Vandervalk, Hamid Mohamadi, Justin Chu, Sarah Yeo, S Austin Hammond, Golnaz Jahesh, Hamza Khan, Lauren Coombe, René L Warren, and Inanc Birol (2017). ABySS 2.0: Resource-efficient assembly of large genomes using a Bloom filter. Genome research, 27(5), 768-777. doi:10.1101/gr.214346.116
Simpson, Jared T., Kim Wong, Shaun D. Jackman, Jacqueline E. Schein, Steven JM Jones, and Inanc Birol (2009). ABySS: a parallel assembler for short read sequence data. Genome research, 19(6), 1117-1123. doi:10.1101/gr.089532.108
Robertson, Gordon, Jacqueline Schein, Readman Chiu, Richard Corbett, Matthew Field, Shaun D. Jackman, Karen Mungall, et al (2010). De novo assembly and analysis of RNA-seq data. Nature methods, 7(11), 909-912. doi:10.1038/10.1038/nmeth.1517
Nielsen, Cydney B., Shaun D. Jackman, Inanc Birol, and Steven JM Jones (2009). ABySS-Explorer: visualizing genome sequence assemblies. IEEE Transactions on Visualization and Computer Graphics, 15(6), 881-888. doi:10.1109/TVCG.2009.116
Create a new issue on GitHub.
Subscribe to the ABySS mailing list, .
For questions related to transcriptome assembly, contact the Trans-ABySS mailing list, .
Supervised by Dr. Inanc Birol.
Copyright 2016 Canada's Michael Smith Genome Sciences Centre