Described in:
Song, L., Sabunciyan, S., Yang, G. and Florea, L. A multi-sample approach increases the accuracy of transcript assembly. Nat Commun 10, 5000 (2019)
Copyright (C) 2018- and GNU GPL by Li Song, Liliana Florea
Includes portions copyright from:
samtools - Copyright (C) 2008-, Genome Research Ltd, Heng Li
Commands, scripts and supporting data for the paper can be found here.
PsiCLASS is a reference-based transcriptome assembler for single or multiple RNA-seq samples. Unlike conventional methods that analyze each sample separately and then merge the outcomes to create a unified set of meta-annotations, PsiCLASS takes a multi-sample approach, simultaneously analyzing all RNA-seq data sets in an experiment. PsiCLASS is both a transcript assembler and a meta-assembler, producing separate transcript sets for the individual samples and a unified set of meta-annotations. The algorithmic underpinnings of PsiCLASS include using a global subexon splice graph, statistical cross-sample feature (intron, subexon) selection methods, and an efficient dynamic programming algorithm to select a subset of transcripts from among those encoded in the graph, based on the read support in each sample. Lastly, the set of meta-annotations is selected from among the transcripts generated for individual samples by voting. While PsiCLASS is highly accurate and efficient for medium-to-large collections of RNA-seq data, its accuracy is equally high for small RNA-seq data sets (2-10 samples) and is competitive to reference methods for single samples. Additionally, its performance is robust with the aggregation method used, including the built-in voting and assembly-based approaches such as StringTie-merge and TACO. Therefore, it can be effectively used as a multi-sample and as a single-sample assembler, as well as in conventional assemble-and-merge protocols.
- Clone the GitHub repo, e.g. with
git clone https://github.com/splicebox/psiclass.git - Run
makein the repo directory
You will find the executable files in the downloaded directory. If you want to run PsiCLASS without specifying the directory, you can either add the directory of PsiCLASS to the environment variable PATH or create a soft link ("ln -s") of the file "psiclass" to a directory in PATH.
PsiCLASS depends on pthreads and samtools depends on zlib.
Usage: ./psiclass [OPTIONS]
Required:
-b STRING: paths to the alignment BAM files; use comma to separate multiple BAM files
or
--lb STRING: path to the file listing the alignment BAM files
Optional:
-s STRING: path to the trusted splice file (default: not used)
-o STRING: prefix of output files (default: ./psiclass)
-p INT: number of threads (default: 1)
-c FLOAT: only use the subexons with classifier score <= than the given number. (default: 0.05)
--sa FLOAT: the minimum average number of supported read for retained introns (default: 0.5)
--vd FLOAT : the minimum average coverage depth of a transcript to be reported in voting (defaults: 1.0)
--stranded STRING: un/rf/fr for library unstranded/fr-firstand/fr-secondstrand (default: not used)
--maxDpConstraintSize: the number of subexons a constraint can cover in DP. (default: 7. -1 for inf)
--primaryParalog: use primary alignment to retain paralog genes (default: use unique alignments)
--tssTesQuantile FLOAT: the quantile for transcription start/end sites in subexon graph (default: 0.5; 1.0 for longest exon margin)
--version: print version and exit
--stage INT: (default: 0)
0-start from the beginning - building the splice site file for each sample
1-start from building the subexon file for each samples
2-start from combining the subexon files across samples
3-start from assembling the transcripts for each sample
4-start from voting the consensus transcripts across samples
Alignment compatibility. PsiCLASS has been tuned to run on alignments generated with the tools HISAT and STAR.
When running PsiCLASS with STAR alignments, run STAR with the option --outSAMstrandField intronMotif, which will include the XS field indicating the strand in the BAM alignments. Further, when including alignments with non-canonical splice sites, use the provided addXS executable to add the XS field:
samtools view -h in.bam | ./addXS reference_genome.fa | samtools view -bS - > out.bam
Trusted introns from other sources. By default, PsiCLASS determines a set of trusted introns from the input spliced alignments, to use in building the global subexon graph. Alternatively, the user can supply an external set of trusted introns, for instance extracted from the GENCODE gene annotations or judiciously selected from the input data using a tool like JULIP. This file must contain three columns:
chr_name start_site end_site
Voting optimization. The default parameters for voting have been calibrated and perform near-optimally for a wide variety of data, including with varying levels of coverage and different library construction protocols. However, if further optimization is desired, to determine a better cutoff value one can run the voting stage (see Usage above) with different parameter values, and assess the performance against a reference set of gene annotations, such as GENCODE. The program 'grader', included in the package, can be used for this purpose. Note that the per sample sets of transcripts will remain unchanged.
Add gene name. For many applications, it would be desirable to associate the known (annotated) gene name with each transcript. PsiCLASS provides the program "add-genename" for such purpose. "add-genename" takes as input a GTF file containing a reference set of gene annotations and a file listing the raw GTF files, and generates a new GTF file for each input raw GTF file by appending the annotated gene names. If a gene is not found in the annotation, "add-genename" will use "novel_INT" to represent its gene name. The program can be run as:
./add-genename annotation.gtf gtflist
The primary input to PsiCLASS is a set of BAM alignment files, one for each RNA-seq sample in the analysis. The program calculates a set of subexon files and a set of splice (intron) files, for the individual samples. (Optionally, one may specify a path to an external file of trusted introns as explained above.) The output consists of one GTF file of transcripts for each sample, and the GTF file of meta-annotations produced by voting, stored in the output directory:
Sample-wise GTF files: (psiclass)_sample_{0,1,...,n-1}.gtf
Meta-assembly GTF file: (psiclass)_vote.gtf
where indices 0,1,...,n-1 match the order of the input BAM files.
Subexon and splice (intron) files, and other auxiliary files, are in the subdirectories:
Intron files: splice/*
Subexon graph files: subexon/*
Log file: (psiclass)_classes.log
The directory './example' in this distribution contains two BAM files, along with an example of a BAM list file. Run PsiCLASS with:
./psiclass -b example/s1.bam,example/s2.bam
or
./psiclass --lb example/slist
The run will generate the files 'psiclass_sample_0.gtf' for 's1.bam', 'psiclass_sample_1.gtf' for 's2.bam', and the file 'psiclass_vote.gtf' containing the meta-assemblies.
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.
This program 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. See the GNU General Public License for more details.
You should have received (LICENSE.txt) a copy of the GNU General Public License along with this program; if not, you can obtain one from http://www.gnu.org/licenses/gpl.txt or by writing to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
Create a GitHub issue.