Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Using TopHat output files with UCSC genome browser

    Hi all,

    Recently, I ran TopHat with 76bp reads data and got the results (sam, bed, and wig files).

    Actual a few lines of my input (fasta file) are:
    >HWUSI-EAS366:4:1:4:624#0/1:
    CTCNGGATGGAGTACAGTGGTGTGATCATGGCTCACTGTAGNNNNNANCN CNTGGGCGCAAGCNNNNNNNNNCTAN
    >HWUSI-EAS366:4:1:4:243#0/1:
    CGGNGCCGTTGCTGGTTCTCACACCTTTTAGGTCTGTTCTCNNNNNCNGN TNCGACTCTCTCTNNNNNANNNCCGN
    >HWUSI-EAS366:4:1:4:1373#0/1:
    GAAAAAACCACCCAGCGGTGATGGCAGCGCGCGTGGGTCCCNNNGNGNGN GGGGCGGGTCGCGCNNNNGNNNCGAN
    >HWUSI-EAS366:4:1:4:1672#0/1:
    GGGCAGGAAAAAAAGGGAAGANAAAATACTGGGGAAGAAAANNNANCNCN GTTTGGCAGCTCTTNNNNGNNNCAGN


    And a few lines of junctions.bed file are:

    track name=junctions description="TopHat junctions"
    gi|29823169|ref|NT_025004.13|Hs18_25160 9690 19656 JUNC00000001 1 + 9690 19656 255,0,0 2 37,38 0,9928
    gi|29823169|ref|NT_025004.13|Hs18_25160 14260 19654 JUNC00000002 2 + 14260 19654 255,0,0 2 57,36 0,5358
    gi|29823169|ref|NT_025004.13|Hs18_25160 19701 160104 JUNC00000003 3 + 19701 160104 255,0,0 2 32,66 0,140337


    A few lines of coverage.wig file are:

    track type=bedGraph name="TopHat - read coverage"
    gi|29823169|ref|NT_025004.13|Hs18_25160 0 9580 0
    gi|29823169|ref|NT_025004.13|Hs18_25160 9580 9655 1
    gi|29823169|ref|NT_025004.13|Hs18_25160 9655 9690 0


    Here is the problem.

    When I copied and pasted the results (either bed file or wig file), I always got an error and when I change the gi|29823169|ref|NT... part to something like chromosome name, it works.

    As you can see from my input file, I don't have gi|29823169|ref|NT... part. I am not sure where the TopHat find such label or reference.

    Can someone tell me what gi|29823169|ref|NT... part means and how I can convert these files into the one that UCSC genome brower understands. I think I need to get the actual chromosome names.

    Thank you,
    Statsteam

  • #2
    The gi lines you see are the fasta file headers from the NCBI human assembly. Each chromosome in that assembly comes in a separate file and it is the accession codes for those separate files that you are seeing.

    The full header for the first accession you found is:

    >gi|29823169|ref|NT_025004.13|Hs18_25160 Homo sapiens chromosome 18 genomic contig, reference assembly

    You therefore need to find all the accessions for the different chromosomes and replace them with the corresponding chromosome name.

    Alternatively you could edit the original fasta files and change the first lines to just contain a chromosome name, eg:

    >chr18

    ..and then reindex the genome and run tophat again. This should put usable chromosome names into your output files.

    Comment


    • #3
      Thank you simon.
      I just started bowtie-build with fasta files containing only chromosome names.

      Statsteam

      Comment


      • #4
        as the output above
        :A few lines of coverage.wig file are:

        track type=bedGraph name="TopHat - read coverage"
        gi|29823169|ref|NT_025004.13|Hs18_25160 0 9580 0
        gi|29823169|ref|NT_025004.13|Hs18_25160 9580 9655 1
        then 9580 has 1 or 0 hit??

        Comment


        • #5
          Originally posted by melody View Post
          as the output above
          :A few lines of coverage.wig file are:

          track type=bedGraph name="TopHat - read coverage"
          gi|29823169|ref|NT_025004.13|Hs18_25160 0 9580 0
          gi|29823169|ref|NT_025004.13|Hs18_25160 9580 9655 1
          then 9580 has 1 or 0 hit??

          No, that is a data column because the output is in bedGraph format.
          When you copy and paste with correct chromosome name, it will draw a bedGraph based on the value of the data column.

          In this example, it will draw 0 for chr18:0-9580 then draw 1 for chr18:9580-9655.

          -Statsteam

          Comment


          • #6
            just to add to this discussion, i found when using sequencing data from mice it worked best for all of my source references to come from UCSC. i used FASTA files for each chromosome downloaded from UCSC's downloads area to build my Bowtie index and I also used UCSC's table browser to produce the GTF file (which i converted to GFF3 using scripts from seq ontology). only when I had built everything from those sources did i have reliable output files that work straight away with the UCSC browser. in fact, when I used the NCBI reference (and swapped the chromosome names out with UCSC's names) the output from Tophat didn't even align with the genome.
            /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
            Salk Institute for Biological Studies, La Jolla, CA, USA */

            Comment


            • #7
              Originally posted by sdriscoll View Post
              just to add to this discussion, i found when using sequencing data from mice it worked best for all of my source references to come from UCSC. i used FASTA files for each chromosome downloaded from UCSC's downloads area to build my Bowtie index and I also used UCSC's table browser to produce the GTF file (which i converted to GFF3 using scripts from seq ontology). only when I had built everything from those sources did i have reliable output files that work straight away with the UCSC browser. in fact, when I used the NCBI reference (and swapped the chromosome names out with UCSC's names) the output from Tophat didn't even align with the genome.
              What is the config file needed to use the Seq ontology script? I can't find the documentation for it.

              Comment


              • #8
                Originally posted by sdriscoll View Post
                just to add to this discussion, i found when using sequencing data from mice it worked best for all of my source references to come from UCSC. i used FASTA files for each chromosome downloaded from UCSC's downloads area to build my Bowtie index and I also used UCSC's table browser to produce the GTF file (which i converted to GFF3 using scripts from seq ontology). only when I had built everything from those sources did i have reliable output files that work straight away with the UCSC browser. in fact, when I used the NCBI reference (and swapped the chromosome names out with UCSC's names) the output from Tophat didn't even align with the genome.
                I have that exact problem but is there a way to fix it if all I have is either the raw file or the bam. or bam.bai files? Do I need to ask my core personnel to realign using the UCSC files? Any help would be greatly appreciated..

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Strategies for Sequencing Challenging Samples
                  by seqadmin


                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                  03-22-2024, 06:39 AM
                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-27-2024, 06:37 PM
                0 responses
                13 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-27-2024, 06:07 PM
                0 responses
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                53 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                69 views
                0 likes
                Last Post seqadmin  
                Working...
                X