This is a test version of Biostars. For the public version, visit https://www.biostars.org.
splice junction information with hisat2

I am newbie to linux and NGS. Can anyone help me out how to get information about splice junctions using HISAT2?

The command I'm using is giving information in single file about alignment in SAM format. The command is as follows:

./hisat2 -p 64 --max-intronlen 10000 -x /data/memona/hisat2-2.1.0/hisat_index -1 /data/memona/SRR959590_A_1P.fq -2 /data/memona/SRR959590_A_2P.fq -S /data/memona/results/hisat_align.sam &
hisat2
• 6,405 views • link
updated by Ram Ram • written by blooming.daisy333 • blooming.daisy333

HI, i am also using same mapping strategy, and getting error "index file does not exist"

(base) macserver@MACs-Mac-Pro hisat2_alingment % hisat2 -p 8 -x /Users/macserver/Desktop/buffalo_refence_genome/buf_rna_index -1 /Volumes/AshishIVRI/OLD/ashishanalysis/Trim_F/Gr_11_R1C.fastq -2 /Volumes/AshishIVRI/OLD/ashishanalysis/Trim_F/Gr_11_R2C.fastq -s /Users/macserver/Desktop/hisat2_alingment/Gr_11.sam    
(ERR): "/Users/macserver/Desktop/buffalo_refence_genome/buf_rna_index" does not exist
Exiting now ...

please suggest possible solution

• 0 views • link
updated by GenoMax GenoMax • written by Ashish • Ashish

What is the output of

ls -lh /Users/macserver/Desktop/buffalo_refence_genome/buf_rna_index*
• 0 views • link
written by GenoMax GenoMax

1 answer

Hi blooming.daisy333,

It is explained in the manual. You have to use the option --novel-splicesite-outfile.
Update:
Be careful there is an error how left splice sites are reported as I mentioned here
The splice sites are correctly reported but the format they use is not necessarily straightforward.

• 1 views • link
written by Juke34 Juke34

Hi Juke,

thanks for the kind guidance. however, im still unable to get the splice junction information despite of using --novel-splicesite-outfile command.

here is the command that i used. kindly point out the mistake and suggest the solution:

./hisat2 --np 0 --pen-noncansplice 10000000 --min-intronlen 20 --max-intronlen 10000 --novel-splicesite-outfile /data/memona/hisat2-2.1.0/result/ --rna-strandness RF --dta -p64 --summary-file -x /data/memona/hisat2-2.1.0/hisat_index -1 /data/memona/Trimmomatic-0.36/SRR959591_E_1P.fq -2 /data/memona/Trimmomatic-0.36/SRR959591_E_2P.fq -S /data/memona/hisat2-2.1.0/result/hisat_align.sam

thank you so much

• 1 views • link
written by blooming.daisy333 • blooming.daisy333

You’re welcome. The problem is you provided a path but no file name. Instead of “/data/memona/hisat2-2.1.0/result/” do “/data/memona/hisat2-2.1.0/result/splice_sites.tsv” and it should be fine.

• 1 views • link
written by Juke34 Juke34

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

• 1 views • link
written by WouterDeCoster WouterDeCoster

thanks for the kind guidance. yes it has produced the output but is giving only 3 fields given below:;

chr1    329728  329839  -
chr1    330066  330757  -
chr1    581256  581357  +

while not mentioning the canonical/non canonical status and doner accepter nucleotides which are important to classify the site.

can you please share the script/tool or command that you used to extract the nucleotide information of doner and accepter site like AT/AC, GT/AT, CT/GC etc...

further is it important to build the index with --ss and --exon options to determine the splice sites...???? or it is OK if the index is built without using these options???

• 1 views • link
written by blooming.daisy333 • blooming.daisy333

I don’t know the use of those options but for sure it is fine without to determine the spicing sites.

For the extraction of spicing sites I have used fasta_domainExtractor_JD.pl from the NBIS/GAAS repository but it is not really adapted for what you want. You can get inspiration from this script to implement what you reallly want.

• 1 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this