Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • chickenmcfu
    replied
    Yes, the description of Deseq is easy to follow, even with no R experience.

    Unfortunately I have run in other problems. Deseq analysis gives me 125 differentially regulated genes. The commercial sequencing service gives me a cuffdiff result of 200 genes (mixed up annotation though). So now I went through the process of cuffdiff and got 129 regulated genes.

    As I now have conducted nearly every possible mapping (bowtie2, tophat2), whole genome or with options for tophat -G and -T and then for cuffdiff -M with the rtRNA.gtf, I could bring it up to 136.

    The only thing in which my analysis differs from them is, that they align straight to only the cds and ncRNAs, completely disregarding the RNAs, but then give cuffdiff a full gtf. They afterwards did a second mapping with the fastqs to rtRNAs, so I know that up to 4.5% map to them.

    Is this right? Can this difference within the alignment give such a huge difference in the results?

    Leave a comment:


  • crazyhottommy
    replied
    I just went through the same process as what you did
    see my post here http://crazyhottommy.blogspot.com/20...-sort-and.html

    I am following the protocol from Simon Anders http://www.nature.com/nprot/journal/....2013.099.html
    Count-based differential expression analysis of RNA sequencing data using R and Bioconductor






    Originally posted by chickenmcfu View Post
    Hey this is my first try to analyse a rna-seq project. Since the company we worked with is not able to give me a usefull annotated differential expression table...

    I just want to know for sure, if its half-way right what I do.

    My samples: 2 conditions, 2 replicats from each condition, 50bp single-end, not strand-specific. I got 23m-50m reads per library.

    I want to know differential expressed genes between conditions.

    I first wanted to use bowtie2 for alignment and that worked pretty well until I noticed that no NH tag for htseqcount is written.

    So I switched to tophat and there it got complicated:

    In default tophat2 finds lesser alignments than bowtie2. Why? As I understand tophat2 uses bowtie2 for aligment.

    As I dont want to find novel junctions, as there is no splicing my bacterium, the final command I used after several attempts is:

    tophat2 -G file.gtf --no-novel-juncs --no-coverage-search --library-typ fr-unstranded index file.fastq

    With every attempt (first: --no-coverage-search; second: added -G; third: added --no-novel juncs) the count of aligned reads dropped a little bit. Why it droppend between this 3 modes?

    With the last mode I got an alignment rate of 65-75%.

    I finally got my count tables with
    samtools view file.bam | htseq-count -t gene -s no - file.gtf > counts.txt

    Now I will use deseq for differential expression.

    Everything ok so far?

    Leave a comment:


  • chickenmcfu
    started a topic Tophat2 Bowtie2 Htseq-count for bacteria

    Tophat2 Bowtie2 Htseq-count for bacteria

    Hey this is my first try to analyse a rna-seq project. Since the company we worked with is not able to give me a usefull annotated differential expression table...

    I just want to know for sure, if its half-way right what I do.

    My samples: 2 conditions, 2 replicats from each condition, 50bp single-end, not strand-specific. I got 23m-50m reads per library.

    I want to know differential expressed genes between conditions.

    I first wanted to use bowtie2 for alignment and that worked pretty well until I noticed that no NH tag for htseqcount is written.

    So I switched to tophat and there it got complicated:

    In default tophat2 finds lesser alignments than bowtie2. Why? As I understand tophat2 uses bowtie2 for aligment.

    As I dont want to find novel junctions, as there is no splicing my bacterium, the final command I used after several attempts is:

    tophat2 -G file.gtf --no-novel-juncs --no-coverage-search --library-typ fr-unstranded index file.fastq

    With every attempt (first: --no-coverage-search; second: added -G; third: added --no-novel juncs) the count of aligned reads dropped a little bit. Why it droppend between this 3 modes?

    With the last mode I got an alignment rate of 65-75%.

    I finally got my count tables with
    samtools view file.bam | htseq-count -t gene -s no - file.gtf > counts.txt

    Now I will use deseq for differential expression.

    Everything ok so far?

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 06-14-2024, 07:24 AM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-13-2024, 08:58 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-12-2024, 02:20 PM
0 responses
17 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-07-2024, 06:58 AM
0 responses
186 views
0 likes
Last Post seqadmin  
Working...
X