Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to detect a replacement gene/region

    Illumina PE sequences.

    I am having trouble finding a program that detects a replacement gene/region. Either I need a different program than the ones have been running (Pindel, Breakdancer, indelMiner) or better command line parameters. The programs detect deletions and short insertions just fine but don't tell me "this gene has been deleted and replaced with a different gene from outside the organism".

    It is easy to see the region in IGV or other visualization program. Visually the region looks like a deletion however unlike a normal deletion where it is possible to see the two PE reads map to either side of the break and see a single read get broke, in the replacement case one side of the PE reads maps while the side does not.

    Any suggestion as to a program to use?

    Thanks,

  • #2
    As an example of what I mean, here are two IGV screenshots.

    The first ("deletion_found") show the reads mapped to the genome. You can see paired reads spanning the deletion (the thin lines between the reads). There are also some reads that are split apart. Not that coverage goes down on either side of the deletion breakpoint.

    The second screenshot shows a similarly sized region (~4500 bases) and has a similar coverage pattern. However note that there are no pairs spanning and that broken reads then to have mismatched bases near the breakpoint. The lack of paired reads is because the pair is part of the replacement gene and since said gene is not part of the reference then no mapping occurs.

    Pindel, for one, picks up the first region (deletion) just fine but skips over the second region (replacement)

    I am using Samtools to pull out places where there are a bunch of mapped but mate unmapped reads. This might work for my purposes. I was just hoping that there would be some program that can pick up replacement genes for me.
    Attached Files

    Comment


    • #3
      In the end I used a combination of cortex as well as pindel for determining indel/replacements. Neither program seems to be perfect and I am not totally satisified with the results in that I have the feeling that I may have missed something. Annotation was done with snpEff.

      BTW: I really like the cortex manual. Well done.

      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
      30 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
      86 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-20-2024, 08:57 AM
      0 responses
      70 views
      0 likes
      Last Post seqadmin  
      Working...
      X