Question: Compute Genetic Map
6
gravatar for Pierre
8.7 years ago by
Pierre480
Spain
Pierre480 wrote:

Dear all,

I would like to get the position in cM for a set of SNPs (SNPs from 1000Genomes Project). What I have is:

  • list of SNPs with their physical position
  • Genotype (or infered genotype) for each SNP for around 50 individuals

Some of these snps (around 25%) are among the HapMap II SNPs used to compute the genetic map available on HapMap webpage (http://ftp.hapmap.org/recombination/2008-03_rel22_B36/rates). But is it enough information to calculate the cM position for 1000genomes SNPs from HapMap II genetic map?

What do you suggest?

Thanks for your help.

Yours truly Pierre

genetics map distance • 5.6k views
ADD COMMENTlink modified 18 months ago by Biostar ♦♦ 20 • written 8.7 years ago by Pierre480

FYI, I could not get that link to work but has now been changed to: ftp://ftp.hapmap.org/hapmap/recombination/2008-03_rel22_B36/rates/

ADD REPLYlink modified 10 weeks ago by RamRS19k • written 8.1 years ago by Zach Stednick650
5
gravatar for Pierre Lindenbaum
8.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum114k wrote:

I hope I understand your question: The UCSC genome database contains the position in both base index and cM for the STS. e.g.:

~> mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -D hg18 -e 'select name,chrom,chromStart,chromEnd,genethonPos from stsMap where genethonPos!=0 limit 2\G'
*************************** 1. row ***************************
       name: AFM280WE5
      chrom: chr1
 chromStart: 3574721
   chromEnd: 3575045
genethonPos: 6.2
*************************** 2. row ***************************
       name: AFM344WE9
      chrom: chr1
 chromStart: 4358261
   chromEnd: 4358654
genethonPos: 11.1

So you can use this STS map as a 'reference' to map your collection of SNP from bp to cM.

ADD COMMENTlink modified 10 weeks ago by RamRS19k • written 8.7 years ago by Pierre Lindenbaum114k
4
gravatar for Steve Buyske
8.1 years ago by
Steve Buyske40
Steve Buyske40 wrote:

The Rutgers Map has an interpolate feature that can be used for any SNP not too near the telomeres or centromere.

http://compgen.rutgers.edu/RutgersMap/InterpolateMap.aspx

ADD COMMENTlink written 8.1 years ago by Steve Buyske40
2
gravatar for Pierre
8.7 years ago by
Pierre480
Spain
Pierre480 wrote:

Hey Pierre,

thanks for your answer! I already got the genetic map provided by USCS database with STS markers. But my issue remains the same: the density of markers of the available genetic maps (the USCS one or HapMap II one) is much lower than the density of my SNP set. For example I have 4054 SNPs between the 2 STS markers above (in your comment).

Therefore I don't know how to attribute to each SNPs of my set a position in cM without doing a very important approximation.

Anyway, what do you suggest? To use a rule of three approach to attribute a cM position to my collection's SNPs?

Thanks for your advice! Pierre

ADD COMMENTlink written 8.7 years ago by Pierre480

I doubt you can get better resolution since centimorgan is experimental unit for recombination frequency. It does not correlate well with base pairs, which depends on genomic position. In human it approximately 1 Mbp in average.

ADD REPLYlink written 8.7 years ago by Yuri1.5k

ok Thanks. I finally adopted this solution which provides satisfying results!

ADD REPLYlink written 8.6 years ago by Pierre480

What solution did you eventually adopt?

ADD REPLYlink written 8.1 years ago by Zach Stednick650

Hi Pierre,

How did you solve it finally? I just have the same problem with you.

ADD REPLYlink modified 10 weeks ago by RamRS19k • written 3.2 years ago by ice4prince10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1166 users visited in the last hour