Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • shorter time for building bowtie index

    Hi:
    I am not sure if this question has been asked before.

    I have 120 RNA-Seq fasta files. I am using tophat2 for aligning to human genome.

    I am queing this job on a 70 node cluster. I am using bowtie genome files downloaded from Illumina genome website.

    Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/genome

    I see for every fastq file alignment, tophat building index file from genes.fa..

    It takes more than 30 minutes for every file.

    Can this time consuming process be avoided since we already have index files for bowtie.

    Following is my command line:

    tophat2 -p 16 \
    -G /Refs/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/genes.gtf \
    -o tpht \
    /Refs/Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/genome \

    myfile_1.fastq myfile_2.fastq

    Thanks

  • #2
    I don't know the answer to your question. Probably not.

    You can however get faster TopHat runs simply by adding the option --no-novel-juncs.

    Comment


    • #3
      Yes, in the most recent version(s) of tophat, you can pre-generate the indexes with the 'transcriptome-index' option:

      Comment


      • #4
        @gringer

        Very interesting, and useful. Thanks for the information.

        Comment


        • #5
          Thank you so much! that is helpful.

          However when I use the following I don't see files - known.gff, known.fa, known.fa.tlst, known.fa.ver and the known.* Bowtie index files in the directory.

          following is my command:

          tophat2 -p 16 \
          -G /Refs/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/genes.gtf \
          -o tpht \
          /Refs/Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/genome \
          --transcriptome-index=known
          myfile_1.fastq myfile_2.fastq

          Is there something wrong. Where would the indexed known files will be deposited?
          thanks
          Adrian

          Comment


          • #6
            To prepare transcriptome indexes you should not include the -o tpht or the sequence files. This is a special one time tophat run (Check the example in the link gringer included)

            Your command would be:
            Code:
            $ tophat2  \
            -G /Refs/Homo_sapiens/Ensembl/GRCh37/Annotation/Genes/genes.gtf \
            --transcriptome-index=IF_you_want_a_directory_name/known \
            /Refs/Homo_sapiens/Ensembl/GRCh37/Sequence/Bowtie2Index/genome
            This needs to be a non-threaded run so you have to omit the -p 16 as well.
            Last edited by GenoMax; 06-30-2014, 09:59 AM. Reason: Corrected error in the genome_base_name

            Comment

            Latest Articles

            Collapse

            • seqadmin
              The Impact of AI in Genomic Medicine
              by seqadmin



              Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
              02-26-2024, 02:07 PM
            • seqadmin
              Multiomics Techniques Advancing Disease Research
              by seqadmin


              New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

              A major leap in the field has
              ...
              02-08-2024, 06:33 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 02-28-2024, 06:12 AM
            0 responses
            27 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-23-2024, 04:11 PM
            0 responses
            74 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-21-2024, 08:52 AM
            0 responses
            82 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-20-2024, 08:57 AM
            0 responses
            69 views
            0 likes
            Last Post seqadmin  
            Working...
            X