Fig. 2

Phylogenetic analysis of 16S rRNA gene sequences of potential P-endosymbionts of the five species of Bhutanese psyllids (emboldened). The evolutionary history was inferred using maximum likelihood based on the Hasegawa-Kishino-Yano model [52]. The tree with the highest log likelihood (− 1375.27) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 0.2691)). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The tree is rooted with ‘Ca. Portiera aleyrodidarum’. (NB OTU numbers in bold classified according to Greengenes database, see Supplementary Table 3; OTU 25406 described in [53])