Gene family evolution
We inferred gene family evolution (expansion and contraction) based on the dating tree using CAFÉ v4.2.1 (Han et al., 2013), and the mode of the single birth-death parameter ‘lambda’ was employed with a significance level of 0.01. PCGs from significantly expanded and species-specific gene families were further enriched for GO and KEGG terms using clusterProfiler with the default parameters (p-value of 0.01, q-value of 0.05).