Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to obtain full length RNA transcript sequence

    Hi everyone,
    i'm new to this kind of tasks so, please be patient!
    I'm trying to create a blast DB using the RNAseq data from ENCODE.
    I've downloaded both the FASTQ reads and the .bam/bai files.
    I need the fasta sequences of all the full length transcripts: is it possible to extract/obtain them from the BAM file?
    Alternatively should i try to do a de novo assembly using Trinity?
    Thanx a lot.
    Regards,

    Davide

  • #2
    I thought this task would have been easy or at least possible since i have the reads aligned to the ref genome (homo sapiens)
    Anyone can help?
    Thanx

    Comment


    • #3
      Which BAM files are you talking about? ENCODE has many.

      Why do you want to make your BLAST data base from RNA-Seq reads rather than simply from, say, the cDNA FASTA file from Ensembl?

      Comment


      • #4
        Originally posted by Simon Anders View Post
        Which BAM files are you talking about? ENCODE has many.

        Why do you want to make your BLAST data base from RNA-Seq reads rather than simply from, say, the cDNA FASTA file from Ensembl?
        I'm talking about BAM file of the human total RNA extract from CSHL Long RNA seq.

        I don't use Ensembl data because cDNA FASTA from Ensembl does not contain all the transcript (i guess) but only "known, novel and pseudogenes" as stated on their website

        Moreover i will probably repeat this task using RNAseq data from cell in particular conditions

        Thanx a lot.

        Davide

        Comment


        • #5
          What you want to do is called reference-based (as opposed to: de-novo) transcript assembly. A tool commonly used for this purpose is cufflinks:

          Roberts, Pimentel, Trapnell, and Pachter:
          Identification of novel transcripts in annotated genomes using RNA-Seq
          Bioinformatics (2011) 27 (17): 2325-2329.
          doi:10.1093/bioinformatics/btr355

          However, before doing this yourself, you may want to check whether the ENCODE people have not already done this analysis. It seems obvious that they would do this.

          I still wonder what you would need a database of all transcripts for. Instead of blasting against it, you can always blast against the genome.

          Comment


          • #6
            Originally posted by Simon Anders View Post
            What you want to do is called reference-based (as opposed to: de-novo) transcript assembly. A tool commonly used for this purpose is cufflinks:

            Roberts, Pimentel, Trapnell, and Pachter:
            Identification of novel transcripts in annotated genomes using RNA-Seq
            Bioinformatics (2011) 27 (17): 2325-2329.
            doi:10.1093/bioinformatics/btr355

            However, before doing this yourself, you may want to check whether the ENCODE people have not already done this analysis. It seems obvious that they would do this.

            I still wonder what you would need a database of all transcripts for. Instead of blasting against it, you can always blast against the genome.
            Thank you.
            I've taken a look at cufflinks, just the manual, but i did not find th FASTA file of the transcript as an output file of some task. Cufflinks instead talk about gtf file as an output (that does not contain the FASTA sequence of the transcript). I'll take a better look to the program.
            I've also read just yesterday this interesting article:
            "Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks" nature protocol

            If i blast in the genome i lose informations that are in the RNA sequence and not in genome (ex. sequences in the transposable element that are not integrated in the genome...)

            Thanx again.

            Comment


            • #7
              You use the GTF file to produce the cDNA FASTA file from the reference FASTA file. This is a simple exercise in script programming.

              Comment


              • #8
                I've taken a look at GTF specs.
                Yes, it is.
                Thanx

                Comment


                • #9
                  Hi i saw the thread..
                  Can i get the logic for the program to create transcripts from genome file.
                  how it differ based on orientation. i mean reads which have positive and negative orientation.
                  Thank you.
                  Deepak

                  Comment


                  • #10
                    Refer to my earlier post to get fasta from Cufflinks GTF.
                    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

                    Comment

                    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, Today, 08:58 AM
                    0 responses
                    8 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, Yesterday, 02:20 PM
                    0 responses
                    14 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 06-07-2024, 06:58 AM
                    0 responses
                    181 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 06-06-2024, 08:18 AM
                    0 responses
                    231 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X