Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Low diversity alignments and GARD (Read 5465 times)
Liz622
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 4
Low diversity alignments and GARD
Sep 16th, 2009 at 9:53am
 
Hello,

My sequences have 0.1%-2.5% diversity, is this too low to use SBP and GARD?  Other recombination tools say you need a minimum diversity (~5%).  I've noticed others publishing using GARD with similar data to mine though, and in the paper about these methods there were some low diversity alignments included, so I decided to give it a try!

I ran sequences from the person with the highest diversity and SBP resulted in Yes for recombination (cAIC 4.25 improvement, 89.7% support), then the GARD said no recomb (used the model selection tool and 4 rate categories).  I'm very new at this, is it normal for the 2 to disagree?  My plan was to use SBP to determine if there is recombination, then follow-up with GARD to find the breakpoints.  Are the differences in my alignments too few to expect findings of recombination?  I did notice in table 1 of  the paper where delta was <100 there were no breakpoints.

The sequences are ~1000 bp and ~15 per alignment.

Thank you!
Liz  
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Low diversity alignments and GARD
Reply #1 - Sep 17th, 2009 at 2:51pm
 
Hi Liz,

A c-AIC difference of 4.25 is borderline (89.7%) significant, and I would say in this case insufficient to support recombination. It is possible to find strong recombination signal in a low diversity sample (assuming the fragments are long enough).

If I recall correctly, the breakpoints in the GARD paper were reported in they were confirmed by another (Kishino-Hasegawa) test for phylogenetic incongruence given the breakpoints.

Cheers,
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
 
Liz622
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 4
Re: Low diversity alignments and GARD
Reply #2 - Sep 21st, 2009 at 8:27am
 
Hi Sergei,

Thank you!  I have one additional question (for now!!).  I read some of the other posts, and am trying to use HyPhy GARDProcessor.bf

I uploaded my alignment, and the [CSV] output from GARD.  It's been running on my (old) computer for 3 days, and it's still fitting the single tree (currently tree_300).  Did I do something wrong, is is caught in a strange loop, or is it just my horribly slow computer? 

I know this is slightly off-topic, but if I need to seek out a faster computer I figure it's best to find out sooner, I've got a couple of these to do, and a deadline  Shocked.
Thank you,
Liz
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Low diversity alignments and GARD
Reply #3 - Sep 21st, 2009 at 10:35am
 
Hi Liz,

You no longer need to run GARDProcessor on your machine. KH tests and other metrics are now reported on the datamonkey.org results page.

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
 
Liz622
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 4
Re: Low diversity alignments and GARD
Reply #4 - Sep 21st, 2009 at 2:21pm
 
Thank you!  I was pretty confused as to why there was a KH report at the bottom of GARD if we had to run it separately, but I'm still at the following directions stage, trying to move into actual understanding!!

Liz
Back to top
 
 
IP Logged