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).