Question: Which reference to map against for RNA-Seq using HISAT2
0
gravatar for skhan
9 months ago by
skhan40
skhan40 wrote:

Hi all,

I am trying to determine differential expression between two test conditions - in rat. The first step would be to map the sequencing data to a rat reference of some sort, and I would like to use HISAT2 (locally installed on my machine) along with the Ensemble rat genome.

Which reference file(s) do I map against?

1) Do I use the hisat2 index for rn6 as is: ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/data/rn6.tar.gz ?

Or

2) Do I need to generate my own hisat2 index from the Ensemble reference genome (ftp://ftp.ensembl.org/pub/release-86/fasta/rattus_norvegicus/dna/Rattus_norvegicus.Rnor_6.0.dna_sm.toplevel.fa.gz) and annotations file (ftp://ftp.ensembl.org/pub/release-86/gtf/rattus_norvegicus/Rattus_norvegicus.Rnor_6.0.86.gtf.gz)

Thanks, skhan

rna-seq hisat2 • 606 views
ADD COMMENTlink modified 9 months ago by Jennifer Hillman Jackson25k • written 9 months ago by skhan40
0
gravatar for Jennifer Hillman Jackson
9 months ago by
United States
Jennifer Hillman Jackson25k wrote:

Hi,

The rn6 reference genome is sourced from UCSC so has different chromosome identifier formats than the Ensembl version, so it would require a different reference annotation file. UCSC has an Ensemble track based on the rn6 build. Check the date on it to see if it is current enough for your needs - it is probably older than the one you linked to.

Using the Ensembl genome/annotation versions are also an option.

For either, install the genome and create the indexes with the Data Managers (fetch genome, sam, picard, hisat2). These are in the tool shed. See the DM tutorial here for how-to: https://galaxyproject.org/admin/#data-resources

The Chrom identifier mismatch FAQ here explains why it is important to make sure that data is all from the same build when doing an analysis (plus how to check/adjust for mismatches, when possible): https://galaxyproject.org/support/chrom-identifiers/

Hope that helps! Jen, Galaxy team

ADD COMMENTlink written 9 months ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 122 users visited in the last hour