Hi HyPhy community,

I am trying to decide between running only PARRIS with 2 synonymous rate classes or also BranchSiteREL and MEME on a set of alignments with 6 taxa.

In general, there is a low level of sequence divergence in these alignments, which is why PARRIS was originally chosen for the positive selection analysis. The average dS is 0.01 among the ~11200 alignments in the dataset.

Given the additional information that can be gleaned from a "branch-site" approach and the reportedly increased power of the new methods BranchSiteREL and MEME, should I also run these new analyses on my genes or at least say on the top 20% most variable subset of them? Or, am I just not going to have enough information in my alignments to estimate the parameters of the more complex models in BranchSiteREL and MEME?

I realize BranchSiteREL and MEME are somewhat complementary to PARRIS, as they're designed to detect episodic as well as pervasive positive selection. However, the issue of lack of information relative to more complex models still applies.

So, what analyses make the most sense to run? Only PARRIS, or
PARRIS plus BranchSiteREL and MEME?

Thanks in advance!
Dan.

Hey all!

Sorry for hijacking the thread, but my question is related to the initial question. I want to find which lineages are subject to episodic diversifying selection using BranchSiteREL, but I have a recombination breakpoint. Is there any workaround to be able to use it with recombinant sequences? Or do you have any recommendation on how proceed?

Thank you very much!

Thanks, Sergei!