StringTie-1.3.3_merge
Rationale and background:
Â
StringTie enables improved reconstruction of a transcriptome from RNA-seq reads
Mihaela Pertea,   Geo M Pertea,   Corina M Antonescu,   Tsung-Cheng Chang,   Joshua T Mendell   & Steven L Salzberg
 doi:10.1038/nbt.3122
Â
StringTie is a fast and highly efficient assembler of RNA-Seq alignments into potential transcripts. It uses a novel network flow algorithm as well as an optional de novo assembly step to assemble and quantitate full-length transcripts representing multiple splice variants for each gene locus. Its input can include not only the alignments of raw reads used by other transcript assemblers, but also alignments longer sequences that have been assembled from those reads.In order to identify differentially expressed genes between experiments, StringTie's output can be processed by specialized software like Ballgown, Cuffdiff or other programs (DESeq2, edgeR, etc.)
Version: 1.3.3
If StringTie is run with the --merge option, it takes as input a list of GTF/GFF files and merges/assembles these transcripts into a non-redundant set of transcripts. This step creates a uniform set of transcripts for all samples to facilitate the downstream calculation of differentially expressed levels for all transcripts among the different experimental conditions. Output is a merged GTF file with all merged gene models, but without any numeric results on coverage, FPKM, and TPM. Then, with this merged GTF, StringTie can re-estimate abundances by running it again with the -e option on the original set of alignment filesÂ
Pre-Requisites
- A CyVerse account. (Register for an CyVerse account here -Â user.cyverse.org)
- Mandatory arguments -
- Input gtf files to merge:Â (in gtf format)
- Optional arguments
- Annotation: provide gtf file or select from the list (a reference annotation file in GTF/GFF3 format can be provided to StringTie)
- output file name for the merged transcripts GTF:merged.out.gtf(output file name for the merged transcripts GTF)
- minimum isoform fraction:0.1(Sets the minimum length allowed for the predicted transcripts)
- minimum assembled transcript length:50(minimum input transcript length to include in the merge)
- minimum reads per bp coverage to consider for transcript assembly:0 (minimum input transcript coverage to include in the merge)
- minimum input transcript FPKM to include in the merge: 0 (minimum input transcript FPKM to include in the merge)
- minimum input transcript TPM to include in the merge : 0 (minimum input transcript TPM to include in the merge)
- keep merged transcripts with retained introns:uncheck (keep merged transcripts with retained introns)
The following test data are provided for testing StringTie1.3.3_merge in here - /iplant/home/shared/iplantcollaborative/example_data/StringTie/StringTie1.3.3_merge:
- reference gtf file- reference.gtf
- Directory of gtf files in (.bam format sorted by their genomic location) -
- Â sample1.gtf
- sample2.gtf
Run StringTie1.3.3_merge on gtf files using reference file.
ResultsÂ
Successful execution of the StringTie1.3.3_merge will contain several files and directories:
- merged.out.gtf: merged GTF file with all merged gene models
More information on the tool can be found here - https://ccb.jhu.edu/software/stringtie/index.shtml?t=manual