# Least common ancestor

**URL:** https://nextstrain.discourse.group/t/least-common-ancestor/1043
**Category:** General
**Created:** [March 17, 2022, 7:14pm UTC](https://nextstrain.discourse.group/t/least-common-ancestor/1043 "2022-03-17T19:14:16Z")
**Posts on this page:** 4
**Page:** 1

<div class="post-metadata">

### Author: ![Patrick](https://avatars.discourse-cdn.com/v4/letter/p/b2d939/32.png) [@Patrick](https://nextstrain.discourse.group/u/Patrick)
#### Post date: [March 17, 2022, 7:14pm UTC](https://nextstrain.discourse.group/t/least-common-ancestor/1043/1 "2022-03-17T19:14:16Z")

</div>

I am creating a tree for my sequences of an avian reovirus. The data is from 1973-2019. The tree keeps creating a common ancestor in the 1700’s which is beyond my data. I am just making sure my tree is correct. If there is a problem how do I fix it.

 ![image](https://canada1.discourse-cdn.com/flex031/uploads/nextstrain/original/1X/4d279cb5a71652df7b85192ad12b7563fc19bef2.png)

Thank you for your help

---

<div class="post-metadata">

### Author: ![corneliusroemer](https://yyz1.discourse-cdn.com/flex031/user_avatar/nextstrain.discourse.group/corneliusroemer/32/175_2.png) [@corneliusroemer](https://nextstrain.discourse.group/u/corneliusroemer)
#### Post date: [March 19, 2022, 1:07pm UTC](https://nextstrain.discourse.group/t/least-common-ancestor/1043/2 "2022-03-19T13:07:22Z")

</div>

What makes you think that this tree is wrong? May the common ancestor not be from a few hundred years ago?

It can indeed be that the inferred common ancestor is wrong - but one would need more information. Is there a time tree in the literature for this virus? So you know what you expect?

It could be that you use a hard coded clock rate that is set too low for this virus. What is your `augur refine` command? In particular, what are the `timetree` parameters?

---

<div class="post-metadata">

### Author: ![Patrick](https://avatars.discourse-cdn.com/v4/letter/p/b2d939/32.png) [@Patrick](https://nextstrain.discourse.group/u/Patrick)
#### Post date: [March 22, 2022, 4:29pm UTC](https://nextstrain.discourse.group/t/least-common-ancestor/1043/3 "2022-03-22T16:29:52Z")

</div>

Thank you for responding and thanks for the help

There doesn’t seem to be a time tree in the literature. The virus is endemic in poultry. It was discovered in the 1950’s. I expect there to be a common ancestor just not that far back.

My refine command in my snakemake file is the following.

rule refine:  
message:  
“”"  
Refining tree  
- estimate timetree  
- use {params.coalescent} coalescent timescale  
- estimate {params.date\_inference} node dates  
- filter tips more than {params.clock\_filter\_iqd} IQDs from clock expectation  
“”"  
input:  
tree = rules.tree.output.tree,  
alignment = rules.align.output,  
metadata = input\_metadata  
output:  
tree = “results/tree.nwk”,  
node\_data = “results/branch\_lengths.json”  
params:  
coalescent = “opt”,  
date\_inference = “marginal”,  
clock\_filter\_iqd = 4  
shell:  
“”"  
augur refine   
–tree {input.tree}   
–alignment {input.alignment}   
–metadata {input.metadata}   
–output-tree {output.tree}   
–output-node-data {output.node\_data}   
–timetree   
–coalescent {params.coalescent}   
–date-confidence   
–date-inference {params.date\_inference}   
–clock-filter-iqd {params.clock\_filter\_iqd}  
“”"

---

<div class="post-metadata">

### Author: ![corneliusroemer](https://yyz1.discourse-cdn.com/flex031/user_avatar/nextstrain.discourse.group/corneliusroemer/32/175_2.png) [@corneliusroemer](https://nextstrain.discourse.group/u/corneliusroemer)
#### Post date: [March 23, 2022, 2:05am UTC](https://nextstrain.discourse.group/t/least-common-ancestor/1043/4 "2022-03-23T02:05:06Z")

</div>

Can’t see any obvious problem in the command.

What does the tree look like in divergence view?

Is there a non-time-tree in the literature?

It could be a tree building mistake but it’s hard to check remotely.

What makes you expect the common ancestor to be recent? Do you have good reasons for that prior?

It could also be that there’s reassortment/recombination happening which would throw off tree building I would imagine.
