Skip to content

Latest commit

 

History

History

tmrca

Folders and files

NameName
Last commit message
Last commit date

parent directory

..
GISAID-HA-genes.txt
GISAID-HA-genes.txt
 
 
README.md
README.md
 
 
TMRCA-estimates-allRuns.txt
TMRCA-estimates-allRuns.txt
 
 
cattle-taxa-set.txt
cattle-taxa-set.txt
 
 
mcc-gtrg-ucld-gmrf-colored.pruned.tre
mcc-gtrg-ucld-gmrf-colored.pruned.tre
 
 
mcc-gtrg-ucld-gmrf-colored.tre
mcc-gtrg-ucld-gmrf-colored.tre
 
 
r1-MP-mcc-gtrg-strict-gmrf.tre
r1-MP-mcc-gtrg-strict-gmrf.tre
 
 
r1-NA-mcc-gtrg-strict-gmrf.tre
r1-NA-mcc-gtrg-strict-gmrf.tre
 
 
r1-NP-mcc-gtrg-strict-gmrf.tre
r1-NP-mcc-gtrg-strict-gmrf.tre
 
 
r1-NS-mcc-gtrg-strict-gmrf.tre
r1-NS-mcc-gtrg-strict-gmrf.tre
 
 
r1-PA-mcc-gtrg-strict-gmrf.tre
r1-PA-mcc-gtrg-strict-gmrf.tre
 
 
r1-PB1-mcc-gtrg-strict-gmrf.tre
r1-PB1-mcc-gtrg-strict-gmrf.tre
 
 
r1-PB2-mcc-gtrg-strict-gmrf.tre
r1-PB2-mcc-gtrg-strict-gmrf.tre
 
 
r1-mcc-gtrg-ucld-gmrf.tre
r1-mcc-gtrg-ucld-gmrf.tre
 
 
r8-100-gtrg-ucld-gmrf-noEPI.xml
r8-100-gtrg-ucld-gmrf-noEPI.xml
 
 

README.md

BEAST tMRCA analysis on H5N1 2.3.4.4b strains from the US, including the strains from the March-April 2024 outbreak in dairy cattle

This analysis uses BEAST v1.10.4: https://beast.community.

File structure

  • mcc-gtrg-ucld-gmrf-colored.tre: a time-scaled maximum clade credibility HA gene tree generated from this analysis.
  • mcc-gtrg-ucld-gmrf-colored.pruned.tre: a paraphyletically subsampled HA phylogenetic tree used in Fig-S2.
  • GISAID-HA-genes.txt: lists the GISAID EPI identifiers for the seqeuences used in this analysis.
  • r8-100-gtrg-ucld-gmrf-noEPI.xml: a non-functional xml file generated by BEAUTI v1.10.4. This file has the GISAID sequence data removed and will not work.

Step-by-step:

  1. Combine the GISAID sequences with the IRMA sequences from this analysis - genbank accessions pending.

  2. input into BEAUTI, parse the dates, and configure the analysis (GTR + gamma with estimated base frequencies; an uncorrelated relaxed lognormal clock; GMRF Bayesian Skyride with time aware smoothlng; uninformative priors; 100 million MCMC with sampling every 10,000)

  3. Analyses were run 8 times, and after inspection the logs and trees were combined (LogCombiner v1.10.4), and a single MCC was generated using TreeAnnotator v1.10.4.