Host haplotype identification

The mitochondrial variable open reading frame (ORF) was amplified with FATP6.1 (5′-TTTGGGSATTCGTTTAGCAG-3′) and RORF (5′-SCCAATATGTTAAACASCATGTCA-3′) primers [38] and submitted to sanger sequencing (GenBank submission ID 2033564). Protein-coding sequences were analyzed using MEGA version 6 [39]. Sequence alignment was performed using MUSCLE. The best model (Kimura-2 parameters with uniform substitution rates) was selected for the lowest BIC (Bayesian Information Criterion). Maximum-likelihood tree was computed with the best model, and the robustness of the tree was tested with 1000 bootstrap replicates.