Evolutionary profiles of β-galactosidase (A), XyG-active β-glucosidase (B) and α-xylosidase (C) genes in green plants. The topologies were obtained by Neighbor-Joining method with genetic distances calculated by p-distance. The bootstrap values higher than 50% are shown for 1000 replicates along with Maximum-likelihood (ML) aLRT values higher than 0.50 (in the figure [NJ/ML]). Accession numbers for all genes in the trees are described in Additional File 5 for β-galactosidase, Additional File 7 for β-glucosidase and Additional File 9 for α-xylosidase. The detailed trees are shown in Additional File 4 for β-galactosidase, Additional File 6 for β-glucosidase and Additional File 8 for α-xylosidase. The tree containing the 100 first blastp hits of AtXYL1 is shown in Additional File 10.