Question: Recombination rate per base per generation
0
gravatar for Alice
20 months ago by
Alice280
USA
Alice280 wrote:

Hello biostars,

I am trying to understand what rho means and how it is calculated in coalescent simulations. We know for Hundon's 'ms' program for example, that cross-over rate parameter is measured by 4N0 * r (r is probability of cross-over per generation between the ends of the locus). Now I want to simulate a locus in ms using -r parameter and I know average recombination rate for my species. How can I calculate per base per generation crossover rate?

Average recombination rate is in cM\Mb, so if, let's say, human recombination rate is 2 cM\Mb, genome size is 3 * 10^9 , then 2 / 1 * 10 ^6 = 2 * 10^-6 recombinations per base pair. This is way higher than a number used by default in different papers (such as 2*10^-8). And if I divide it by 3 * 10^9, the number is so small that it does not make sense at all. I know my calculations are incorrect, but I just don't see another way to derive this parameter.

coalescence ms recombination psmc • 1.6k views
ADD COMMENTlink modified 19 months ago by rjwang30 • written 20 months ago by Alice280
3
gravatar for rjwang
19 months ago by
rjwang30
Indiana University
rjwang30 wrote:

Hey,

I don't know if you've already figured this out, but aside from the complexities of estimating recombination rate - the big problem in your estimate is conflating cM/Mb with crossovers/Mb. The genetic distance of 1 centimorgan represents only a 1% chance of crossover in the region. If you take 2cM/Mb to be (2 x 10^-2 crossovers) / (1 x 10^-6 bp), you get the 2 x 10^-8 crossovers/bp.

ADD COMMENTlink modified 19 months ago • written 19 months ago by rjwang30
1
gravatar for Alice
20 months ago by
Alice280
USA
Alice280 wrote:

Answering myself.

TL;DR - it is hard to get this number.

Here is the paper about recombination rates and pop gen models: https://www.nature.com/articles/nrg1227

ADD COMMENTlink written 20 months ago by Alice280
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: 1687 users visited in the last hour