Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Unannotated sequence

    Hello all,

    I have a mystery sequence that I was hoping y'all could help me with.
    I've assembled a transcriptome set from both illumina and 454 datasets of the same plant sample. My approach to functional annotation has been to use blast searches against nr, refseq and hmmer searches vs interpro.
    After performing expression analysis, I found that my 3rd highest expressed contig a long with a couple more very highly expressed ones were unannotated.
    This 1200 bp sequence is found near identical in both my illumina and 454 assemblies. It has no hits vs nt, nr, refseq, cdd, interpro... it also seems to exist as the linking part of an assembly chimeric of 2600bp
    (kinase-unannotated mystery seq-transferase).
    Could anyone suggest ideas? I'm at a loss. Thanks very much!!

    >c1
    gttttatattagaattgacaaaaacataataataaaaaggttgtgtaactaaagatggcacttattgaaacaacattggctccgaatattactaaccataaccacaaccacggcgtttaagtggtgaacacagtataggtagaaacaaaatccataacatagcactaggaaaccctagaaaaacagggagacagagatgatgatcctctcctctccaaaagcattaccatccacactccagtcaccatcggtgctcaaaattattttactttctttcwsttaattgcctttgttcatcctcaactctctctttcctcccttttcttgcaattcaaattagtattccaatggccactgctgaggttgtatctgcagcgactgcattgcaagagaaagacacaactcatgaggaattgaacaagagcccagttgttgatgagaccaaggaagagaagccaacagaagaagtggtgacaccaccacccacatcagaagaggtcaaggaagacaaggctgatgcttcaatagaggaaccagtagccactacagatcaagcagaagccactgctgaagaagagaaagcagaggaggcacaagttgaggaggtaaaggaaacaaaggattcagttgaagaggagaaagcagtggttgaggagactaaagaagaggaatcaaaagaagataaggttagtactcctgaaccagtagcacctgaagagaagacccacgaaactacacctactactactaaagatgttagtgagagtactgttgaagcagaagagaaagttgttcaatcagagataccagttgaggaagccaaagcaacagaagcagaagagaaagttgttgccgcagagactactccagtagttgagaaggctgaggagtagatttgatctgtgttccccaggattctgattgtctggccagctagtagtgctggatgtttgtgtagtgggtgttttatatagaagcttgcatgttaagagttgatgagttttaatgtagtaaaagaatctatgtaggttatgaggcttggaagttagtatattaacttggtattttcttgggcccaccaagtaccattctgccatgtggcaaaggtccaatgccccaaaaatatgtttcaaattcccaaagcttttgtttggtggagggacactctctctattgcagtgtggagggttgtgtactggttggtctcgtgtactgtggatvdwcatsrsagttgcagcwacatyyraagccaaagagtaatcacagattttaaaagtaagtgttgggccctg

  • #2
    A functional but non-protein coding RNA? Maybe something in the ENCODE database could help you.

    Comment


    • #3
      Thanks for your suggestion Jeremy. I'll inquire with ENCODE.
      I would have thought that non coding rna would have nt annotations at least.

      Comment


      • #4
        Hello all,

        I have some 454 data and in that there are approx 138 sequences that I am not able to found in Mouse Genome. the are unplaced or uncharacterized.

        So, Please help me to find them.??????//

        Comment


        • #5
          Originally posted by ikim View Post
          Hello all,

          I have a mystery sequence that I was hoping y'all could help me with.
          I've assembled a transcriptome set from both illumina and 454 datasets of the same plant sample. My approach to functional annotation has been to use blast searches against nr, refseq and hmmer searches vs interpro.
          After performing expression analysis, I found that my 3rd highest expressed contig a long with a couple more very highly expressed ones were unannotated.
          This 1200 bp sequence is found near identical in both my illumina and 454 assemblies. It has no hits vs nt, nr, refseq, cdd, interpro... it also seems to exist as the linking part of an assembly chimeric of 2600bp
          (kinase-unannotated mystery seq-transferase).
          Could anyone suggest ideas? I'm at a loss. Thanks very much!!

          Code:
          >c1
          gttttatattagaattgacaaaaacataataataaaaaggttgtgtaactaaagatggcacttattgaaacaacattggctccgaatattactaaccataaccacaaccacggcgtttaagtggtgaacacagtataggtagaaacaaaatccataacatagcactaggaaaccctagaaaaacagggagacagagatgatgatcctctcctctccaaaagcattaccatccacactccagtcaccatcggtgctcaaaattattttactttctttcwsttaattgcctttgttcatcctcaactctctctttcctcccttttcttgcaattcaaattagtattccaatggccactgctgaggttgtatctgcagcgactgcattgcaagagaaagacacaactcatgaggaattgaacaagagcccagttgttgatgagaccaaggaagagaagccaacagaagaagtggtgacaccaccacccacatcagaagaggtcaaggaagacaaggctgatgcttcaatagaggaaccagtagccactacagatcaagcagaagccactgctgaagaagagaaagcagaggaggcacaagttgaggaggtaaaggaaacaaaggattcagttgaagaggagaaagcagtggttgaggagactaaagaagaggaatcaaaagaagataaggttagtactcctgaaccagtagcacctgaagagaagacccacgaaactacacctactactactaaagatgttagtgagagtactgttgaagcagaagagaaagttgttcaatcagagataccagttgaggaagccaaagcaacagaagcagaagagaaagttgttgccgcagagactactccagtagttgagaaggctgaggagtagatttgatctgtgttccccaggattctgattgtctggccagctagtagtgctggatgtttgtgtagtgggtgttttatatagaagcttgcatgttaagagttgatgagttttaatgtagtaaaagaatctatgtaggttatgaggcttggaagttagtatattaacttggtattttcttgggcccaccaagtaccattctgccatgtggcaaaggtccaatgccccaaaaatatgtttcaaattcccaaagcttttgtttggtggagggacactctctctattgcagtgtggagggttgtgtactggttggtctcgtgtactgtggatvdwcatsrsagttgcagcwacatyyraagccaaagagtaatcacagattttaaaagtaagtgttgggccctg
          By eye one can see that this sequence very low complexity (repetitive). Did you run your BLAST search with the default parameters? The default parameters for BLAST mask low complexity regions prior to the search to speed it up by avoiding low information hsps. Using the low complexity filter for this sequence is means very little of the sequence will actually be searched.

          Run a BLASTX of this sequence against the NCBI nr protein database with the low complexity filter disabled. You will see a very large number of hits to proteins in nr. Now most of those are annotated as hypothetical or predicted but at least it's a place to start.

          Comment


          • #6
            Thanks very much for your advice Kmcarr and bringing up the low complexity issue. That makes a lot of sense, DUST or Seg was on during my initial searches.
            I'd run repeatmasker on this sequence though with settings to determine low complexity regions but it did not report anything. And unfortunately, turning off filtering on the NCBI blast server did not return any additional results for me.

            Comment


            • #7
              Originally posted by ikim View Post
              And unfortunately, turning off filtering on the NCBI blast server did not return any additional results for me.
              Sorry, forgot to mention that you also need to disable Composition Based Statistics. If you are using the NCBI web BLASTX select "No adjustment" from the pop-up for Compositional adjustments. If you are running blastx on the command line use "-comp_based_stats 0" (that's a zero, not the letter 'O').

              Disabling low complexity filtering and composition based statistics your sequence returned > 100 significant hits (E ≤ 1e-10) vs. the nr protein database. What is being matched is a general pattern of very glutamate rich proteins.

              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, Yesterday, 06:37 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X