Skip to main content

Table 1 Positive selection in the glycolysis enzymes from Lanceocercata stick insects

From: Positive selection in glycolysis among Australasian stick insects

Gene pp n (m) bp ө k ω (M0) M0:M3 M7:M8 M8:M8a
Galm −1 28 (2) 681 0.11 1.86 0.16 210.93** 6.08* 0.68
Pgi1 2 20 (3) 1452 0.06 2.80 0.10 292.65** 21.92** 15.33**
Pgi2 2 28 (4) 675 0.08 4.69 0.10 216.97** 16.55** 4.59*
Tpi 5 28 (1) 651 0.09 2.37 0.14 53.06** 1.27 0.95
Gapdh 6 29 (2) 837 0.09 2.60 0.04 126.85** 10.78* 0.24
Pgk 7 29 (1) 915 0.13 2.70 0.19 283.07** 3.34 0.57
Pglym 8 27 (0) 486 0.10 3.67 0.08 153.05** 2.10 0.68
Eno 9 25 (1) 903 0.09 1.99 0.11 131.51** 1.89 0.01
Pyk 10 29 (2) 777 0.10 3.07 0.04 146.00** 0.00 0.00
Pdh +1 27 (2) 729 0.09 3.95 0.05 103.29** 0.46 0.00
  1. Selection inferred using likelihood ratio tests between nested site-specific models.
  2. Eno = Enolase; Galm = Aldose 1-epimerase; Gapdh = Glyceraldehyde 3-phosphate dehydrogenase; Pdh = Pyruvate dehydrogenase; Pgi1 = Phosphoglucose isomerase 1; Pgi2 = Phosphoglucose isomerase 2; Pgk = Phosphoglycerate kinase; Pglym = Phosphoglycerate mutase; Pyk = Pyruvate kinase; Tpi = Triosephosphate isomerase; PP = pathway position of the enzyme out of the ten key glycolysis reactions, N.B. Galm and Pdh are the steps immediately before and after the pathway. n = number of species with complete and partial (m) sequences; bp = sequence length; ө = nucleotide diversity; k = ratio of transitions to transversions; ω(M0) = ratio of synonymous to nonsynonymous substitutions under the M0 model; 2∆l is given for each model comparison (MO:M3; M7:M8 and M8:M8a), which is twice the difference between the log likelihood of the two nested site-specific models implemented in PAML. Models are judged to have a significantly better fit (* = P-value < 0.05; ** = P-value < 0.001) based on the χ2 distribution with degrees of freedom proportional to the difference in the number of parameters between models; M7/M8 = 2, M8/M8a = 50:50 mixture of point mass 0 and 1.