For all BEAST analyses, data partitioning and nucleotide substitution models were determined using PartitionFinder 2.1.1 [34, 35]. The Markov chain Monte Carlo chains were run for 100 million generations, sampling every 10,000 generations. Tracer v1.6  was used to assess appropriate burn-in and the adequate effective sample size values (> 200). A burn-in of 25% was applied, and the maximum clade credibility (MCC) tree with the mean ages and 95% HPD intervals on nodes were conducted in TreeAnnotator v1.8.2 (part of the BEAST package) and edited in FigTree v.1.4.2 (http://beast.bio. ed.ac.uk/FigTree).