Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Script for creating GO annotation file from Interproscan output

    Dear comunity,

    I am working with a non-model organism, thus I have to use alternative approaches for many of the analyses. At the moment I want to run the GO enrichment analysis of my differentially expressed transcripts using BiNGO, and for that I have to load my own annotation file of the transcriptome assembly which I created with Interproscan. The output of IPRS can be found here: https://github.com/ebi-pf-team/inter...example-output

    So basically, using the IPRS output I have to create one association per line of the transcript ID that is in the 1st column with the GO term that is 14th column. The issue is that for many transcripts there are different GO terms associated, while others have none, as for instance:

    transcript_1 ...columns_2-13... GO:0004601|GO:0006979|GO:0020037|GO:0055114
    transcript_1 ...columns_2-13... GO:0004601|GO:0006979|GO:0020037|GO:0055114
    transcript_1 ...columns_2-13...
    transcript_1 ...columns_2-13... GO:0004601|GO:0055114
    transcript_1 ...columns_2-13... GO:0004601|GO:0042744
    transcript_2 ...columns_2-13...
    transcript_2 ...columns_2-13... GO:0055085

    And here is how the the custom annotation file should be:

    transcript_1 = 0004601
    transcript_1 = 0006979
    transcript_1 = 0020037
    transcript_1 = 0055114
    transcript_1 = 0042744
    transcript_2 = 0055085

    Please, can someone help me with that? It wouldn't be a problem if the output of the script generates reduntant lines really, I can remove duplicated values later.

    Best regards,

    Gustavo

  • #2
    Hi Gustavo,

    I've attached a little python script that should do what you want. I had to name it "rearrange_go.txt" because of seqanswers restrictions - just rename it "rearrange_go.py".

    First open the python file in a text editor and replace INPUT_FILE_NAME with the name of your file. Then put the script in the same directory as your file and run the following:

    python rearrange_go.py

    Another file called 'custom_annotation.txt' should be created. The script assumes you have python installed (mac and linux usually do by default). It also assumes that column 14 is always GO terms, and not anything else, though it can be blank.

    Give it a go and let me know if it works!

    Cheers,

    Matt.
    Attached Files

    Comment


    • #3
      Hi Matt,

      It worked, thanks a lot!!! You saved my day xD

      Cheers,

      Gustavo

      Comment


      • #4
        Excellent! Glad it worked

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Essential Discoveries and Tools in Epitranscriptomics
          by seqadmin


          The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
          Yesterday, 07:01 AM
        • seqadmin
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        39 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        41 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        35 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        55 views
        0 likes
        Last Post seqadmin  
        Working...
        X