Supplementary Information

Table S1. Complete TiLV genome data used in this study.
Table S2. Summary of alignment regions used in this study.
Table S3. Root-to-tip regression analyses.
Table S4: tMRCA and evolutionary rate estimates of TiLV (median with 95% HPD interval) under the uncorrelated relaxed molecular clock assumption with lognormal rate distribution.
Table S5. tMRCA and evolutionary rate estimates of TiLV (median with 95% HPD interval) under the strict molecular clock and the Skygrid coalescent prior.
Table S6. Additional sequences from individual segments.
Table S7. tMRCA and evolutionary rate estimates of individual segments with additional sequences (median with 95% HPD interval) under the strict molecular clock.
Figure S1. RT-PCR amplified TiLV segments (from segments 1 to 10) from representative infected sample in 2019. M, DNA marker (New England Biolabs).
Figure S2. Time-calibrated phylogenies under the relaxed molecular clock assumption with uncorrelated lognormal rate distribution.
Figure S3. Demographic histories under a relaxed molecular clock with uncorrelated lognormal rate distribution.
Figure S4. Time-calibrated phylogenies under the Skygrid coalescent model and strict molecular clock. The time cutoff value was 14 years.
Figure S5. Demographic history under the Skygrid coalescent model and strict molecular clock. The time cutoff value was 14 years.
Figure S6. Time-calibrated phylogenies of segments 1, 3, 5, and 9 with additional sequence data under the strict molecular clock.
Figure S7. Demographic histories of segments 1, 3, 5, and 9 with additional sequence data under the strict molecular clock.
Figure S8. Time-calibrated phylogenies of segments 1, 3, 5, and 9 with additional sequence data under a relaxed molecular clock with uncorrelated lognormal rate distribution.
Figure S9. Demographic histories of segments 1, 3, 5, and 9 with additional sequence data under a relaxed molecular clock with uncorrelated lognormal rate distribution.