Tweet
I wanted to try tophat 1.1 with a UCSC supplied GTF file, but the binary (linux-x86_64) I downloaded keeps asking for a GFF file. Has anyone had success running this version with a GTF file?
-
-
I have had succes with the 1.1.0 binary. Are you sure you are running the updated TopHat and not an old one left on your system? You can check using the --version option when you run TopHat ($ tophat --version). -
It's working! I had some deprecated paths in my scripts which were causing the problem. Thanks for your help!Comment
-
where can I get reliable GTF annotation files?
I'm trying to get the GTF files for hg18, preferably at isoform level, from UCSC browser portal, in order to run with Tophat and cufflinks, but apparently I can't find such files there.
So far I've managed to download a table from here http://genome.ucsc.edu/cgi-bin/hgTables?command=start but I'm not quite sure if it's the right way to do it. How did you guys get the file?
thanksComment
-
You are correct in using the table browser. To download a GTF file of a track, you need to select GTF in the output format dropdown menu, type a name for the output file, for instance hg18.UCSCknowngene.isoforms.gtf, then click get output. That should get you a GTF file.Comment
-
yes that what I did, but the problem is that there are loads of options for group/track/table and I don't find any single combination to look significantly more appealing. I'm looking for a detailed annotation for hg18, currently I have taken, the GTF file with:
Group: Genes and Gene prediction Tracks
Track: RefSeq Genes
Table: refGene
but I dunno how sensible/standard choice it is! and I guess it does not contain isoform level annotation.Comment
-
RefSeq gene does include a lot of isoforms (any that have RefSeq mRNA entries), but there are certainly isoforms expected to be missing.Comment
-
So you can supply TopHat with a GTF file of annotated transcripts, which, using the --GTF option, will be the first place where reads are mapped, followed by the whole genome, with or without novel junction discovery in this second stage. As I understand it, this is after TopHat 1.4.
I'm curious to know how t was before 1.4. I think you could already give TopHat a GTF file, but it used it second. Am I right? If so, what is the difference between using it [the GTF file] first and using it second after the genome?Comment
-
Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.
[Article Coming Soon!]...Today, 08:07 AM -
Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...09-23-2024, 06:35 AM -
During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.
Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...09-09-2024, 10:59 AM
Comment