Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Interpreting MEME results (Read 10558 times)
Vinny Lynch
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 14
Interpreting MEME results
Oct 27th, 2013 at 8:02am
 
Hi Sergei,

I'm curious how one should interpret MEME results when no sites are identified with dN>dS (i.e., in the table where  synonymous and non-synonymous substitution rates over sites are shown), but when a particular site (which we a priori know is a function altering aa change) has post. prob. of 0.995 and empirical Bayes factor of selection >1000 along the lineage where we know the function altering aa change occurred?

In other words does one interpret this situation as no evidence of diversifying selection at this site across all species (from the FEL model fit), but some evidence for an episodic selection event at that site in a particular lineage (from the MEME model fit)?

Thanks
Vinny
Back to top
 
 
IP Logged
 
Vinny Lynch
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 14
Re: Interpreting MEME results
Reply #1 - Oct 28th, 2013 at 8:21am
 
(With the obvious caution against placing too much importance on the finding of selection at any particular branch-site combination.)
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Interpreting MEME results
Reply #2 - Oct 31st, 2013 at 3:13pm
 
Hi Vinny,

Am I correct to assume that the p-value from MEME is significant for site X, but the p-value from FEL is NOT significant for the same site?

What MEME tells you is that there is evidence for selection at a proportion of branches, inferred as a part of the model, along that site (i.e. there is episodic but NOT pervasive positive diversifying selection).

The actual branches under selection at that site are really only good for exploratory data analysis; there is no strong statistical significance you could or should derive from this result.

Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged
 
Vinny Lynch
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 14
Re: Interpreting MEME results
Reply #3 - May 15th, 2014 at 9:52am
 
Hi Sergei,

Thanks for the reply. Actually I'm curious what it means when the p-value is not significant (in either MEME or FEL), but sites nevertheless have high PP and BF of selection?

We would never use the MEME results to argue that any particular aa change was selected in any particular lineage. But are using the results in exploratory analyses, we have functional data indicating a couple sites are functionally divergent in different species and were curious what MEME might infer was going on there.
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Interpreting MEME results
Reply #4 - May 15th, 2014 at 10:12am
 
Hi Vinny,

The PP and BF for the exploratory MEME analysis are not terribly reliable and that's because Empirical Bayes (used to find them) does not account for variability in point estimates.
For example, if a site has a point estimate of beta- = 0.1 and beta+ = 2, but the variance of beta+ is high (e.g. we can confidently bound it away from 1), EB will take the point estimate of 2 at face value, even though it could have also been anything between 0.7 and 4 (hypothetically speaking).

Exploratory analysis is another matter -- here you probably want to err on the side of sensitivity versus specificity, so you suggested use is probably OK.

Sergei

Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged