Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • thedamian
    Member
    • Feb 2012
    • 50

    Find genes in contigs based on homology proteins

    Hello,
    I have contigs from a new organism (its genome isn't known well).
    I know, that this organism produces some proteins that are similar to proteins of other organisms.

    I would like to find an approximate location of a gene that produces a protein similar to other homology protein.

    Do you know some good tools or pipelines for that purpose?
    A have tried tblastn and SWiPS however I am little confused with the result and don't know how to interpert it
    (how to get a linear sequence of supposed gene)

    Does anybody have some good advices or an experience in this area?
  • krobison
    Senior Member
    • Nov 2007
    • 734

    #2
    How many contigs? How big are they? How were they generated?

    BLASTX is a straightforward way to approach this. Make sure you set a low limit on the number of overlapping hits; otherwise a lot of hits in one region will prevent hits from being reported in another region.

    If you have very long contigs & think that frameshifts are rare, you should use a gene finder apropos to the domain of life you are in (for example, Glimmer3 is good for eubacteria & archea but doesn't model introns).

    You should start by reading some genome sequencing papers of similar organisms; this will give you an idea of the sorts of tools that might be applied.

    Comment

    • thedamian
      Member
      • Feb 2012
      • 50

      #3
      Originally posted by krobison View Post
      How many contigs? How big are they? How were they generated?

      BLASTX is a straightforward way to approach this. Make sure you set a low limit on the number of overlapping hits; otherwise a lot of hits in one region will prevent hits from being reported in another region.

      If you have very long contigs & think that frameshifts are rare, you should use a gene finder apropos to the domain of life you are in (for example, Glimmer3 is good for eubacteria & archea but doesn't model introns).

      You should start by reading some genome sequencing papers of similar organisms; this will give you an idea of the sorts of tools that might be applied.
      I have 27050 contigs that have length of 1000 or more bases. The longest is 90691. So lenghts are good. They were generated by De Novo assembler embedded in CLC Genomic Workbench (reads are from MiSeq).
      My organism is roundworm so introns exist.

      I know some basic theory, but I'm struggling with applying it.

      Comment

      • A_Morozov
        Member
        • Feb 2011
        • 40

        #4
        If you only need several genes, you can gather them by hand. Download some sequence browser (I personally prefer artemis), look up hit location, strand and frame in blastx/tblasn output and copypaste aminoacid sequence. Of course, if gene contains introns you will get some uncertainty near exon boundaries, but at least you will have most of the sequence and hopefully active sites.
        If, on the other hand, you want to find ALL genes in a genome, considering all splicing and detecting pseudogenes and whatnot, I advise you to google JGI and TIRG gene finding/annotation pipelines. But that will cost a significant amount of computational power and you will probably want to outsource this part of work to bioinformatician if you don't have one in your lab.

        Comment

        • thedamian
          Member
          • Feb 2012
          • 50

          #5
          I want olny five genes, so this simple approach seems to be fine. I'll check it

          Comment

          Latest Articles

          Collapse

          • seqadmin
            New Genomics Tools and Methods Shared at AGBT 2025
            by seqadmin


            This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

            The Headliner
            The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
            03-03-2025, 01:39 PM
          • seqadmin
            Investigating the Gut Microbiome Through Diet and Spatial Biology
            by seqadmin




            The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
            02-24-2025, 06:31 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-20-2025, 05:03 AM
          0 responses
          17 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-19-2025, 07:27 AM
          0 responses
          18 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-18-2025, 12:50 PM
          0 responses
          19 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-03-2025, 01:15 PM
          0 responses
          186 views
          0 reactions
          Last Post seqadmin  
          Working...