Table 3

Results of PAML analyses for the primate dataset
Data set 2ΔlnL df P-value Proportion of sites with ω > 1 (average ω for these sites) Positively selected codonsa
All NAT coding sequences (n= 43)
 M1a vs. M2a 0 2 P = 1.0 none
 M7 vs. M8 3.49 2 P = 0.17 none
 M0 vs. free ratios 114.89 80 P= 0.006
Orthologous sequences to human NAT1 (n= 19)
 M1a vs. M2a 0 2 P = 1.0 none
 M7 vs. M8 0 2 P = 1.0 none
 M0 vs. free ratios 49.09 34 P= 0.04
Orthologous sequences to human NAT2 (n= 19)
 M1a vs. M2a 8.48 2 P= 0.01 0.038 (ω = 4.8)
 M7 vs. M8 11.24 2 P= 0.004 0.041 (ω = 4.6) 191*, 173*
 M0 vs. free ratios 27.52 34 P = 0.77

n, number of sequences; 2ΔlnL, twice the log-likelihood difference of the models compared; df, degrees of freedom; P-value, level of significance (P-values < 0.05 are shown in bold).

a Sites pinpointed to be under positive selection under the selection model with the highest likelihood (M8 for the set of NAT2 orthologous sequences) by Bayes Empirical Bayes (BEB) analysis. Only sites with posterior probability greater than 0.95 (*) or 0.99 (**) are shown. Sites are numbered according to the full human coding sequence.

Sabbagh et al.

Sabbagh et al. BMC Evolutionary Biology 2013 13:62   doi:10.1186/1471-2148-13-62

Open Data