This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ensembl gene annotation GRCh37.87

Hi,

I need to quantify gene expression by salmon so I need Ensembl gene annotation GRCh37.87 likely in fasta

I tried ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/ but not working

Do you know know from where I can downlioad such file for salmon?

Thank you

rna-seq ensembl salmon
• 6,729 views • link
updated by GenoMax GenoMax • written by zizigolu zizigolu

If Salmon needs the FASTA sequences for the transcripts, then you can do the follow these steps (http://ccb.jhu.edu/software/stringtie/gff.shtml#gffread_ex) using the GRCh37 reference from GENCODE and getting the GTF file from possibly from ftp://ftp.ensembl.org/pub/grch37/release-87/gtf/homo_sapiens/Homo_sapiens.GRCh37.87.chr.gtf.gz. This assumes your BAM files are alignments to the main assembly and not including the alternative haplotypes or patches.

• 0 views • link
written by jean.elbers jean.elbers

I have bam files I have alignment come from GRCh37_g1k by STAR; Now I need to quantify raw counts; I have done that by featurecounts but I have a lot of strange features so I decided to used Salmon. Thank you anyway but I don't know why I can not open these links

So not I though to convert my bam to fastq

When I used this command in Salmon I obtained this error

salmon quant -t gencode.v30lift37.transcripts.fa -l A -a file.bam -o salmon_quant

If you have access to the genome FASTA and GTF used for alignment
consider generating a transcriptome fasta using a command like:
gffread -w output.fa -g genome.fa genome.gtf
you can find the gffread utility at (http://ccb.jhu.edu/software/stringtie/gff.shtml)

Finally I used gffread but I am getting this error

[fi1d18@cyan01 fi1d18]$ /temp/hgig/fi1d18/gffread-0.11.2.Linux_x86_64/gffread -w transcripts.fa -g hs37d5.fa gencode.v30lift37.annotation.gff3

Warning: couldn't find fasta record for 'chr1'!
Error: no genomic sequence available (check -g option!).
[fi1d18@cyan01 fi1d18]$

[fi1d18@cyan01 fi1d18]$ /temp/hgig/fi1d18/gffread-0.11.2.Linux_x86_64/gffread -w transcripts.fa -g hs37d5.fa gencode.v30lift37.annotation.gtf

Warning: couldn't find fasta record for 'chr1'!
Error: no genomic sequence available (check -g option!).
[fi1d18@cyan01 fi1d18]$
• 0 views • link
written by zizigolu zizigolu

1 answer

Use GRCh37 fasta from GENCODE. If you need annotations then those are available there as well.

That said for new analyses you should stick with current release unless you are trying to reproduce some past analyses.

• 0 views • link
written by GenoMax GenoMax

Thank you, bam files are from GRCh37 so I need that I guess

• 0 views • link
written by zizigolu zizigolu

Sorry @genomax

Where I can find the same genome.fasta and GTF files from GRCh37_1k (1000 genome) ?

I googled, for genome.fasta there were 3 phases but I did not find any corresponding GTF

• 0 views • link
written by zizigolu zizigolu

Log in to answer this question.

More posts like this
  • Mapping RNA seq reads from bed file to genes
    written by jrose • jrose
    I have a bed file with columns: chrom, read_start, read_end, id, strand where every row represents a read from an RNA seq data set (e.g. …
  • VEP 100 and GRCh37
    written by cocchi.e89 cocchi.e89
    I've installed latest VEP distribution (100). It asked me for installing FASTA file and I Saw it installed `homo_sapiens/100_GRCh38` file, I was wondering where I …
  • convert GRCh38 FASTA .fa to .vcf
    written by LTDavid • LTDavid
    How do you convert the GRCh38.fa FASTA file to a VCF file? The file is from Ensembl primary_assembly [ ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/][1] [1]: ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/
  • Which Ensembl chromosome sequence file is used for the gff3 files?
    written by O.rka O.rka
    **I can't which of the following files were used to create the GFF3 files?** ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.dna.toplevel.fa.gz ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.dna.primary_assembly.fa.gz or an entirely different one from the ftp: ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/ …
  • Fasta length does not match Ensembl info
    written by firestar firestar
    I am calculating the total genome length from a fasta file using the following code `zcat genome.fa.gz | grep -v ">" | wc | awk …
  • hg38 Ensembl GTF file
    written by Gregor Rot Gregor Rot
    Dear all, what is the difference between the hg38 Ensembl GTF file: ftp://ftp.ensembl.org/pub/release-90/gtf/homo_sapiens/Homo_sapiens.GRCh38.90.gtf.gz and ftp://ftp.ensembl.org/pub/release-90/gtf/homo_sapiens/Homo_sapiens.GRCh38.90.chr.gtf.gz (note the .chr. addition in the name) Thank you, Gregor
  • No XIST in Emsembl cdna.fa
    written by WouterDeCoster WouterDeCoster
    I noticed that Salmon didn't count XIST in my RNA-seq data, and tracked it down to the transcriptome fasta from Ensembl which was used to …
  • peak annotation chip-seq
    written by ana222 • ana222
    Hello all. I used the Sample file PeakAnalyzer (PeakAnnotator.peaks file) (http://www.ebi.ac.uk/sites/ebi.ac.uk/files/groups/bertone/software/PeakAnalyzer_sample_data.tar.gz) And I downloaded human annotation from the site of Ensembl. (ftp://ftp.ensembl.org/pub/grch37/release-89/gtf/homo_sapiens/Homo_sapiens.GRCh37.87.gtf.gz) . Now running …
  • NCBI vs ENSEMBL FASTA
    written by uki_al • uki_al
    Hi I have a question about the differences between the FASTA files that can be downloaded from the ensembl ftp (ftp://ftp.ensembl.org/pub/release-75/fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa.gz) and the ncbi ftp …
  • Using Ensembl ID History Converter on GRCh37 Transcripts for salmon
    written by stefanos.bamopoulos • stefanos.bamopoulos
    Hello guys, I have a question regarding the use of salmon with the GRCh37 ensembl reference. For my analysis I run salmon for the purpose …