Hi All. Ensembl has the reference genome and annotations that I need separated into files by chromosome. I'd like to combine these into one file with an additional column to specify which chromosome they are found on, so that I can align and generate feature counts for the organisms whole genome, from which I can do differential expression and gene newtwork analysis for the organisms whole genome.
Does anyone know of a way I can do this? Software, scripts?
1 answer
I have downloaded several genomes and annotations from Ensembl, each contained in single genome (fasta) and annotation (gff or gtf) files. It may be a bit confusing at first, with all the .dna.chromosome.1.fa.gz and .chromosome.X.gff3.gz, but you just want the files with no chromosome in their names.
For example, for pig, you will want the Sus_scrofa.Sscrofa11.1.dna.toplevel.fa.gz for reference genome, and Sus_scrofa.Sscrofa11.1.101.gff3.gz for the gff annotation (the gtf directory doesn't split the annotation by chromosome).
Log in to answer this question.
-
RNA-seq analysis, wrong genome build
written by a_bis •
-
How to align and visualize data with .fasta and .gff3 files in IGV?
written by Phuc •
-
Rsubread featureCounts returns no matches
written by devarts •
-
GO and KEGG from DESeq2 counts
written by devarts •
-
read counts per chromosome arm
written by Sara
-
HISAT2 index generation
written by deshpande.neha2 •
-
Warnings while generating STAR indices
written by skhan •
-
RNA-seq analysis: How to output normalized counts or normalized FPKM values for all genes in each c…
written by sclittlerain •
-
Differential expression in edgeR with cuffdiff count output
written by mjg •
-
How to run MISO with ENSEMBL GTF if reads were aligned against UCSC genome
written by Anil Kesarwani
What have you tried? Shell loops in conjuction with
awkcan do this if used well.I didn't have any initial ideas, other than running alignment and feature counting for each chromosome individually and then trying to use deplyr to combine data tables in R, as I was planning on doing the differential expression with EdgeR in R. But, I'll try using awk and shell loops to do this prior to alignment. Thank you. Any additional hints would be greatly appreciated.