Fig. 4

Unrooted phylogenetic trees of Wolbachia based on 2 and 6 markers by Maximum Likelihood. a Analysis based on concatenation of ftsZ and coxA; the total length of datasets is 1217 bp. The topology was inferred using Maximum Likelihood (ML) inference using IQ-TREE [39]. The Best-fit model calculated using ModelFinder according to BIC index was TIM3 + I + G4. b Analysis based on concatenation of 16S rDNA, dnaA, ftsZ, coxA, fbpA and gatB; the total length of datasets is 3481 bp. The Best-fit model calculated using ModelFinder according to BIC index was GTR + I + G4. Nodes are associated with Bootstrap values based on 1000 replicates, only bootstrap value superior to 70 are indicated. The Wolbachia supergroups (A–S) are indicated. The bolded names indicate data produced in this study and in the case of wApol, the asterisk (*) indicates sequences produced by PCR amplification