search

gaynorr / AlphaSimR

R package for breeding program simulations
https://gaynorr.github.io/AlphaSimR/
Other
43 stars 20 forks source link

General Question About QTL Effect #202

Closed zjwinn closed 1 month ago

zjwinn commented 1 month ago

Hi again,

This is more of a question rather than a bug or comment on a new feature. I am trying to pull the "true" effect of QTL out of the simulation parameters. I do this by going into the SP object and this allows me to pull something out of the slot called "addEff" which I am assuming are the additive effects of the QTL assigned for the trait. I am trying to compare this to GWAS estimate results. I am wondering, how is the effect calculated and what allele is considered the reference allele? I am assuming that the effect is calculated in reference to the major allele (e.g., 0) and the effect represents the additive effect of inheriting a copy of the alternate allele.

So, If I am comparing the estimated effect of an MTA identified by GWAS and the "true" additive effect of the QTL, then I need to make sure the effect is always calculated in reference to the reference allele (0) and expressed as inheriting the alternate allele (2).

Is this assumption correct?

Thanks, Zach

gregorgorjanc commented 1 month ago

Allele 0 in AlphaSimR is assumed to the the ancestral allele if you have used runMacs(). If you have used quickHaplo() then the 0 and 1 have no real meaning. At any rate, which allele is 0 and 1 does not matter in terms of the additive effect - if you swap the 0 and 1 you just have to swap the sign of the additive effect;)

If you are simulating just an additive trait, then "addEff" is the thing you are after for comparing with GWAS, but note that "addEff" is a simulated effect (pulled from the assumed trait model), hence it's assumed known, and therefore does not need to be calculated. genParam() does return theoretical regression coefficients (alpha, allele substitution effects), which are akin to what you get from GWAS, but it really depends on what GWAS model you are fitting. Also, these true effects are calculated one at a time without incurring linkage&LD, so it's a theoretical quantity.

As to comparing between AlphaSimR and GWAS, yes you should be mindful how GWAS software you use will treat 0 and 1 (also be mindful that some software might consider 1 as major if it's more frequent). Do you really need to focus on the sign though? You could just look at the magnitude of the absolute effects too!?

@gaynorr another one to be moved to discussion forum;)

zjwinn commented 1 month ago

Oh yes, now that I see the discussion tab, this would definitely fit in that category. So, if I am pulling the "addEff" slot from "SP" for the specific trait, then that is the additive effect in either direction. I did start with runMacs() so I will go forth with the assumption that 0 is the ancestral allele and 1 is the alternate allele. It is true that I could look at the absolute value of the effect and compare the magnitudes. I think I will go forward with this assumption.