Genes with evolutionary rates that are significantly associated with a change in carnivory score
| Gene | Rho | Parametric P-value | Parametric Q-value | Permulation P-value | Permulation Q-value |
|---|---|---|---|---|---|
| SLC14A2 | −0.446 | 1.54 × 10−8 | 2.14 × 10−4 | <1.00 × 10−5 a | <0.035a |
| SLC13A2 | −0.358 | 6.94 × 10−6 | 0.037 | <1.00 × 10−5 a | <0.035a |
| CLDN16 | −0.400 | 7.93 × 10−6 | 0.037 | 1.00 × 10−5 | 0.035 |
| CPB1 | −0.325 | 8.02 × 10−5 | 0.150 | 1.00 × 10−5 | 0.035 |
| PNLIP | −0.350 | 1.85 × 10−5 | 0.065 | 2.00 × 10−5 | 0.046 |
| ACADSB | −0.306 | 3.36 × 10−4 | 0.260 | 2.00 × 10−5 | 0.046 |
[i] Permulation P-values represent the proportion of 100,000 permulations that produced a stronger correlation with the change in carnivory score than the observed value for each gene. Multiple hypothesis testing corrections were performed by generating Q-values using Storey's correction method (Storey et al. 2020) (FDR = 0.05). A negative correlation (Rho) signifies the following pattern: the greater the decrease in carnivory, the higher the rate of evolution of the gene.
[ii] aAfter generating 100,000 null statistics, none produced a stronger correlation with diet than the observed values for SLC13A2 and SLC14A2. However, if the P-values are adjusted to the smallest observed nonzero P-value (1.00 × 10−5), they would produce a significant empirical Q-value (FDR = 0.05).