We use the tophat-cufflinks-cuffmerge-cuffdiff pipeline for RNA Seq analysis. I am wondering if there is a way to have "cufflinks" - the command not the package - just do the transcript assembly since we never use the fpkm outputs from it anyway?
-
-
Just run Cufflinks alone. Cufflinks provides a transcript.gtf file which consists of all possible transcripts. You could use gffread then and provide your reference genome and this gtf to obtain a transcript fasta file. For your case, you just wouldn't continue with any of the other modules they have for downstream analysis. Is this what you were asking for? -
@cedance: Thanks for your response. My post perhaps was not clear. After the cufflinks run, we use cuffmerge to generate merged gtf and then run cuffdiff on it on each of our bam files. But, I see that cufflinks also spends time doing abundance estimates using the GTF file provided for guidance assembly (RABT) - I dont think these abundance results are used by any downstream application. I was wondering if there was a way to avoid this.Comment
-
NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...Yesterday, 06:46 AM -
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM
Comment