Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
mouse or non-human cancer sequence download Peggy Liang Bioinformatics 2 05-13-2016 04:36 PM
Cosmic database analysis m4merg Bioinformatics 1 02-23-2015 11:40 AM
Mouse somatic variants (COSMIC equivalent) id0 Bioinformatics 0 10-31-2013 07:49 AM
Mouse mRNA Database Htao Bioinformatics 1 01-12-2013 03:58 PM
Genomic Sequence Database vs Chromosomal sequence database ? Arupsss Bioinformatics 2 06-20-2012 10:22 AM

Thread Tools
Old 08-05-2016, 03:58 AM   #1
Junior Member
Location: Oslo

Join Date: Mar 2016
Posts: 2
Default Compare mouse cancer sequence data with COSMIC database

I'm working with whole exome sequence data from a mouse breast cancer model. I've done mutation calling using the GATK best practices workflow (BWA mem (mm10) -> Picard deduplicate -> base quality score recalibration -> Mutect2. Annotated .vcfs with SnpEff). I now wanted to compare my variants to the COSMIC data base to see if the orthologous variants have been found in human, and am having trouble figuring out the most appropriate way to do this.

My first thought was to convert my VCF files to .bed files and use the UCSC liftover tool to convert the files from m10 to hg19. I've however also found a publication in which they perform what I'm looking to do, and in their methods they write:

Variants were compared to known human somatic mutations as available via the COSMIC database38. Briefly, the mouse and human sequences for homologous proteins were pairwise aligned using Clustal Omega39 and the human protein position homolo- gous to the mouse mutation was used to query COSMIC for known missense and nonsense mutations at or surrounding this peptide position. Local conservation was determined after sequence alignment using a 6 10 amino acid residue window surrounding the substituted amino acid.
I've looked into Clustal Omega, but I'm having trouble understanding how to implement it into this sort of workflow.

I'm wondering if anyone has done this sort of analysis, and what workflow you've had the most success with?

CLFougner is offline   Reply With Quote
Old 08-07-2016, 11:08 PM   #2
tristan dubos
Location: France

Join Date: Dec 2015
Posts: 39

Hi ,
I looked the Clustal Omega and you can only compared homologous by pair ? and the think is how to automatize it ?

tristan dubos is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:21 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO