Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • synonymous snps from vcf

    Hi

    Anyone has any ideas how would one find out (programmatically) synonymous and non-synonymous snps from vcf files? I have used mpileup on several hundred bacterial genomes to get the vcf file.

    Thanks

  • #2
    Well, you either write your own tool to do that or try annovar

    Comment


    • #3
      Or Ensembl's VEP (http://www.ensembl.org/tools.html) or snpEff (http://snpeff.sourceforge.net/) or...

      Comment


      • #4
        Thanks guys but....

        I am trying to program it myself and I thought I can get some leads into how to do this from a vcf file.

        What do you think of this quick way of doing this:

        1- get the nucleotide sequence of the CDS that has the SNP
        2- perform 6-frame translation
        3- compare with reference translated sequence
        4- if the sequences are different then the SNP at point (1) is non-syn if they are the same then its syn.

        Not accurate but will give you an idea. What do you guys think?

        Comment


        • #5
          You can try comparing the coordinates in the variant VCF with the coding region start/ends in refseq to see where your variant falls in and make a determination based on that.

          Comment


          • #6
            What I've done is using the coordiante from the vcf to get the sequence around and including the SNP. Then I blastx those sequences against a database of the proteins from that bacterium. Then I parse the blastx to find out which changes cause amino acid differences.

            But yes, annovar is easier, if you can get a file for annovar to use to compare to.

            Comment


            • #7
              Originally posted by bioinfun View Post
              Hi

              Anyone has any ideas how would one find out (programmatically) synonymous and non-synonymous snps from vcf files? I have used mpileup on several hundred bacterial genomes to get the vcf file.

              Thanks
              bioinfun, I have a similar problem. Are there any solutions now? Thanks.

              Comment


              • #8
                Look at this pub. "De novo Transcriptome Assembly and SNP Discovery in the Wing Polymorphic Salt Marsh Beetle Pogonus chalceus (Coleoptera, Carabidae)"

                I now provide a quote from the primary author, reference their paper if you use the script

                "The script for finding amino acid changes uses several data files.

                - I searched the ORFs in the unigenes with this program: http://proteomics.ysu.edu/tools/OrfPredictor.html

                è Output: a CDS file (DNA sequences of the ORFs) and a PEP file (AA sequences of the ORFs, and also contains START, STOP and READINGFRAME of the ORFs)



                - SNP calling with SAMtools

                è Output: VCF file (SNP and positions of SNP)



                - Perl script (SNP_in_ORF_nonsyn.pl) infers whether SNPs are located within an ORF and whether the SNP results in an amino acid change. The script gets the SNP position from the VCF file, mutates the position in the original sequence in the unigene fasta file, then translates that sequence according its ORF (from PEP file) and then checks whether the original sequence differs from the mutated sequence. The script uses bioperl.

                è Output: each line in the VCF file that contains a nonsynonymous SNP. At the end, the number of synonymous and nonsynonymous is also outputted.



                I made the script and data available here: http://users.ugent.be/~slvbelle/NGS/

                (I added an example PEP and VCF file which should work)



                The script should be used as follows:

                ./SNP_in_ORF_nonsyn.pl Trinity_GC018ALL_unique.fasta PEP.fasta SNP.vcf > output"

                Comment


                • #9
                  Originally posted by JackieBadger View Post
                  Look at this pub. "De novo Transcriptome Assembly and SNP Discovery in the Wing Polymorphic Salt Marsh Beetle Pogonus chalceus (Coleoptera, Carabidae)"

                  I now provide a quote from the primary author, reference their paper if you use the script

                  "The script for finding amino acid changes uses several data files.

                  - I searched the ORFs in the unigenes with this program: http://proteomics.ysu.edu/tools/OrfPredictor.html

                  è Output: a CDS file (DNA sequences of the ORFs) and a PEP file (AA sequences of the ORFs, and also contains START, STOP and READINGFRAME of the ORFs)



                  - SNP calling with SAMtools

                  è Output: VCF file (SNP and positions of SNP)



                  - Perl script (SNP_in_ORF_nonsyn.pl) infers whether SNPs are located within an ORF and whether the SNP results in an amino acid change. The script gets the SNP position from the VCF file, mutates the position in the original sequence in the unigene fasta file, then translates that sequence according its ORF (from PEP file) and then checks whether the original sequence differs from the mutated sequence. The script uses bioperl.

                  è Output: each line in the VCF file that contains a nonsynonymous SNP. At the end, the number of synonymous and nonsynonymous is also outputted.



                  I made the script and data available here: http://users.ugent.be/~slvbelle/NGS/

                  (I added an example PEP and VCF file which should work)



                  The script should be used as follows:

                  ./SNP_in_ORF_nonsyn.pl Trinity_GC018ALL_unique.fasta PEP.fasta SNP.vcf > output"
                  I haven’t tried it but I think your script on the end of the post should work (definitely cite that PLoS paper). I also wonder if there are alternative ways because my case is much simple. I sequenced a long and heterogeneous viral ORF using HiSeq 2000. Thus ORF prediction is unnecessary. My destination is to calculate the number of dS and dN over the viral ORF through a sliding window. Only tool I am aware is CLC’s SNP analysis tool from a publication (http://www.ncbi.nlm.nih.gov/pubmed/22278255). There may other facilities to be able to do this job too. Thanks in advance.

                  Comment


                  • #10
                    JackieBadger, I tried the script. It came with:

                    Use of uninitialized value $countSyn in concatenation (.) or string at SNP_in_ORF_nonsyn.pl line 101, <GEN0> line 39393.
                    Use of uninitialized value $countNonSyn in concatenation (.) or string at SNP_in_ORF_nonsyn.pl line 102, <GEN0> line 39393.

                    any advice? pls.

                    Comment


                    • #11
                      SNPdat

                      SNPdat can be used for this

                      Background Single nucleotide polymorphisms (SNPs) are the most abundant genetic variant found in vertebrates and invertebrates. SNP discovery has become a highly automated, robust and relatively inexpensive process allowing the identification of many thousands of mutations for model and non-model organisms. Annotating large numbers of SNPs can be a difficult and complex process. Many tools available are optimised for use with organisms densely sampled for SNPs, such as humans. There are currently few tools available that are species non-specific or support non-model organism data. Results Here we present SNPdat, a high throughput analysis tool that can provide a comprehensive annotation of both novel and known SNPs for any organism with a draft sequence and annotation. Using a dataset of 4,566 SNPs identified in cattle using high-throughput DNA sequencing we demonstrate the annotations performed and the statistics that can be generated by SNPdat. Conclusions SNPdat provides users with a simple tool for annotation of genomes that are either not supported by other tools or have a small number of annotated SNPs available. SNPdat can also be used to analyse datasets from organisms which are densely sampled for SNPs. As a command line tool it can easily be incorporated into existing SNP discovery pipelines and fills a niche for analyses involving non-model organisms that are not supported by many available SNP annotation tools. SNPdat will be of great interest to scientists involved in SNP discovery and analysis projects, particularly those with limited bioinformatics experience.




                      (there is also a short tutorial in the downloads section)

                      You only need a VCF for input, annotation file (GTF) and reference sequence (Fasta file). The annotation and sequence information can be from your own assembly and dont require any preprocessing.

                      Comment


                      • #12
                        Originally posted by fanx View Post
                        JackieBadger, I tried the script. It came with:

                        Use of uninitialized value $countSyn in concatenation (.) or string at SNP_in_ORF_nonsyn.pl line 101, <GEN0> line 39393.
                        Use of uninitialized value $countNonSyn in concatenation (.) or string at SNP_in_ORF_nonsyn.pl line 102, <GEN0> line 39393.

                        any advice? pls.
                        Please check again: http://users.ugent.be/~slvbelle/NGS/

                        I made some modifications, it should work now...

                        Comment


                        • #13
                          JackieBadger, I tried the script. It came with:
                          Error, Reference nucleotide does not equal the one in the original sequence at ./SNP_in_ORF_nonsyn_multiSNP.pl line 85, <GEN0> line 6.

                          any suggestions, please?

                          Maoshigua

                          Comment


                          • #14
                            Hi maoshigua,

                            can you send me ([email protected]) a sample of your data? I will try to fix it.

                            Cheers,
                            Steven

                            Comment


                            • #15
                              Hi Steven,
                              i send you those three input files. thanks a lot.

                              Maoshigua

                              Comment

                              Latest Articles

                              Collapse

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

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 03-14-2024, 06:13 AM
                              0 responses
                              32 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-08-2024, 08:03 AM
                              0 responses
                              71 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-07-2024, 08:13 AM
                              0 responses
                              80 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-06-2024, 09:51 AM
                              0 responses
                              68 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X