Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Options for PGM Variant Analysis

    I am preparing to order barcoded primers for a sequencing project on the PGM, and would appreciate some feedback on the options for data analysis before I make this purchase. We are sequencing 4 loci in approximately 65 individuals (diploid, mammal, WGS available). The amplicons will be ~350 bp. We were told this could be done on a single 314 chip.

    We are going through a sequencing center, so I do not have access to the Ion Torrent server. If I understand correctly, I could pay to use the Variant Caller through Amazon though. What are the general costs? Does this work when not using the human genome as a reference?

    Does anyone have experience with any free programs for variant analysis of PGM data? How about GATK? I hear so much about the need to stick with platform-specific analyses to deal with their particular types of sequencing errors, and just wanted to know what the general thoughts were on that topic.

  • #2
    I think Nils Homer develop the TMAP software, which is a variant of his BFAST short read aligner that is geared toward the Ion Torrent sequencing error profile, and post-sequencing data analysis for variant calls.

    Regards,

    -Tom Blomquist

    Comment


    • #3
      You could probably suggest to the sequencing center to have them install a plug in called PathogenAnalyzer. It's on the community. It looks to me to be designed specifically for viral/bacterial amplicons and viral whole genomes. However, I read the manual where you can put in different reference sequences (unlimited?) that it will map against and pick best fit. It also does variant calling. I assume you could fool it into using small regions of the human genome as applicable. I haven't used it myself, so can't help much more.

      Comment


      • #4
        The center should be able to run the variant caller that's part of the server. There is no need for amazon web services. All the torrent server tools are are open source and posted to github.

        Comment


        • #5
          Thank you everyone for your suggestions, and for the tip that everything's available on github. I watched a webinar yesterday for the Ion Reporter Software, in which the instructor stated that only the human genome can be used as a reference. I will be working with non-primate species, so that may rule out using their variant caller.

          I looked into TMAP and found that it can be combined with the Ion Variant Hunter while specifying your own reference, so that may be an option. PathogenAnalyzer also sounds like it could be useful, so I will explore these further. Thanks!

          Comment


          • #6
            TMAP is a software suite. It utilizes 4 algorithms. By default on the torrent server tmap4 is run, which is a variant of bwa-mem. The others are:
            tmap1 - bwa-short
            tmap2- bwa long read (bwa-sw)
            tmap3 - variant of SSAHA

            You can pick which ones you want to run if you compile it yourself OR you can run any combination of the the 4 with tmap mapall command.

            That said, I've had a wicked time trying to track down who actually upkeeps tmap at lifetech right now. Where did you hear that Nils Homer developed it? I think it was someone internal at Ion Torrent (life tech).

            Comment


            • #7
              Email communication from dr. Homer to myself in Oct. 2012.


              "Check out TMAP on github for non SOLID reads. It is my next software. N"


              I took it to mean it was his development. This email was in response to an email chain regarding the gearing of BFAST index seeds for the ion torrent error profile. It may not actually be a BFAST code base, but I think it is Nils' work based on my email strings.

              Btw. Did I mention how much I love the customizability of BFAST? :-)

              -Tom Blomquist

              Comment


              • #8
                Interesting.... I've never used it actually. Just bwa or bowtie2. Did they get BFAST to work better with Ion Torrent data?

                Comment


                • #9
                  Originally posted by jtackney View Post
                  Interesting.... I've never used it actually. Just bwa or bowtie2. Did they get BFAST to work better with Ion Torrent data?
                  Honestly, wouldnt know the difference between tmap and bfast. I ended up sticking with BFAST because of how reproducible my sequencing alignment calls were. I work primarily with my own custom amplicon libraries and parsing my reads into subregions being analyzed by unique/custom index masks geared for each fasta reference database, I can get my error rate for alignments into the 1:100,000 and 10 million reads processed on an i7 intel chip with 32gb ram and 8 cores in 1-2 hours. I'm working toward developing some acute care panels so speed, accuracy and reproducibility are of essence.

                  -Tom

                  Comment


                  • #10
                    Originally posted by thomasblomquist View Post
                    Honestly, wouldnt know the difference between tmap and bfast. I ended up sticking with BFAST because of how reproducible my sequencing alignment calls were. I work primarily with my own custom amplicon libraries and parsing my reads into subregions being analyzed by unique/custom index masks geared for each fasta reference database, I can get my error rate for alignments into the 1:100,000 and 10 million reads processed on an i7 intel chip with 32gb ram and 8 cores in 1-2 hours. I'm working toward developing some acute care panels so speed, accuracy and reproducibility are of essence.

                    -Tom
                    BFAST and TMAP are not related, as the latter is derived from BWA and customized for Ion Torrent data, while the former was primarly used for SOLiD data alignment and the early days of Ion Torrent alignments. I was the primary developer of both bfast and the extension of BWA (Heng Li is the author of BWA.

                    URLs:

                    Torrent Mapping Alignment Program. Contribute to iontorrent/TMAP development by creating an account on GitHub.


                    Originally posted by jtackney View Post
                    TMAP is a software suite. It utilizes 4 algorithms. By default on the torrent server tmap4 is run, which is a variant of bwa-mem. The others are:
                    tmap1 - bwa-short
                    tmap2- bwa long read (bwa-sw)
                    tmap3 - variant of SSAHA

                    You can pick which ones you want to run if you compile it yourself OR you can run any combination of the the 4 with tmap mapall command.

                    That said, I've had a wicked time trying to track down who actually upkeeps tmap at lifetech right now. Where did you hear that Nils Homer developed it? I think it was someone internal at Ion Torrent (life tech).
                    I did develop it, although I am no longer working at Life Technologies.

                    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
                    49 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