Question: StringTie's merge function
0
gravatar for Ammu
11 days ago by
Ammu0
Ammu0 wrote:

Is the stringtie's merge used to merge the transcripts fronm biological replicates prior to differential expression analysis?

rna-seq • 144 views
ADD COMMENTlink modified 9 days ago by Kevin Blighe7.3k • written 11 days ago by Ammu0
0
gravatar for Kevin Blighe
9 days ago by
Kevin Blighe7.3k
Republic of Ireland (√Čire)
Kevin Blighe7.3k wrote:

It can be used for that if you want - yes.

I believe that it was designed for situations where you want to just create a consensus transcriptome reference GTF from multiple samples, whatever they may have been. Take a look at the release notes:

1/3/2016 - v1.2.0 release

  • new feature: a "transcripts merge" usage mode of StringTie which is triggered by the new --merge option; with this option StringTie expects as input a list of GTF files and merges/assembles all these transcripts into a non-redundant set of transcripts. This performs a
    function similar to the CuffMerge script in the Cufflinks/Tuxedo suite -- please see the updated protocol in the manual.

  • an improved "estimation only" usage mode (-e option); in this mode only, StringTie now assigns reads to partially covered reference
    transcripts as well, therefore relaxing the previous requirement that reference transcripts need to be covered end-to-end in order to be
    found as expressed.

  • improvements in dealing with abnormally large bundles by filtering of likely spurious spliced alignments produced by some aligners (e.g. STAR); note that HISAT2, used with the --dta option, is now the recommended aligner to use for StringTie.

[source: https://ccb.jhu.edu/software/stringtie/]

ADD COMMENTlink modified 9 days ago • written 9 days ago by Kevin Blighe7.3k

This term of multiple samples is what I am not clear about. Does it mean the biological replicates or different genotypes ? In my case, I have biological replicates, different genotypes with different tretment to identify differential gene expression.

ADD REPLYlink modified 9 days ago • written 9 days ago by Ammu0
1

It can be absolutely anything that you want. The logic of why you are merging/comparing the transcriptomes needs to make sense. StringTie is just providing functionality for you to use in whichever way(s) you please.

Technically, you could merge the replicates independently with each other, and then further perform a merge of all 'replicate transcriptomes'

ADD REPLYlink written 9 days ago by Kevin Blighe7.3k
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 2.3.0
Traffic: 963 users visited in the last hour