Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Genome Annotation Pipeline Help Required

    First post as a user here, so please go easy on me for lack of due diligence

    We have a genome assembled from Illumina data. There is a reference genome of a closely related species (same genus). I downloaded the proteins from this reference genome and sought to map them to our genome using local tblastn, as a homology-based annotation (we also have predicted transcripts from MAKER as an ab initio annotation method).

    I have seen this method used in the literature, but all the method descriptions skip an important step - actually physically mapping the best tblastn hits (from whatever criteria) to the genome.

    I assume there is some way to convert the blast xml output to an annotation file (GFF or similar) - one that conserves the info from the blast (especially protein name and function). I tried looking into Biopython and BioPerl but could not lay hands on the proper method of doing this.

    Can someone please point me in the right direction?

  • #2
    Originally posted by marct View Post
    First post as a user here, so please go easy on me for lack of due diligence

    We have a genome assembled from Illumina data. There is a reference genome of a closely related species (same genus). I downloaded the proteins from this reference genome and sought to map them to our genome using local tblastn, as a homology-based annotation (we also have predicted transcripts from MAKER as an ab initio annotation method).

    I have seen this method used in the literature, but all the method descriptions skip an important step - actually physically mapping the best tblastn hits (from whatever criteria) to the genome.

    I assume there is some way to convert the blast xml output to an annotation file (GFF or similar) - one that conserves the info from the blast (especially protein name and function). I tried looking into Biopython and BioPerl but could not lay hands on the proper method of doing this.

    Can someone please point me in the right direction?
    I'm currently doing something similar using gmap.. I have a transcriptmome that I'm mapping to a genome though.. (we assembled the transcriptome and then the genome of a related species was subsequently released) also have ests that I'm mapping to a genome with gmap. Exonerate I believe does a similar job and has a protein matching mode.. Anyway both these programs will output in gff format

    Comment


    • #3
      Thanks for the reply. As of this moment I am running genBlast and exonerate as well as the tblastn. In all of these cases, I am using the protein database of the model species as the query and my genomic sequence as the target (or database). I'll let you know how it goes.

      Comment


      • #4
        How come the functions of MAKER were not sufficient for your analysis?

        Comment


        • #5
          I am joining this project somewhat in the middle of the process.

          From my understanding, the initial Maker run was ab initio only, we do not have ESTs or RNA-seq data to add to the pipeline. So while important, the SNAP/Augustus etc gene calls from Maker should constitute one line of evidence for our annotations, while direct alignment of homologous proteins coupled with splice-site detection (a la exonerate) should constitute another, homology-based line of evidence.

          Stop me if I'm wrong.

          Comment


          • #6
            Hang on..

            I thought MAKER was just for ab initio. You can use that to bring together ESTs and RNA-Seq data too?

            That's what I have. I'm just mapping my transcripts to the ESTs at the moment

            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, 06:54 AM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 06-14-2024, 07:24 AM
            0 responses
            16 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 06-13-2024, 08:58 AM
            0 responses
            15 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 06-12-2024, 02:20 PM
            0 responses
            17 views
            0 likes
            Last Post seqadmin  
            Working...
            X