Higher-order interactions enhance the latitudinal tree diversity gradient

The global decrease in species diversity from low to high latitudes is among the most robust biogeographic patterns^1,2. There is continuing debate on the contribution of conspecific negative density dependence (CNDD) to the latitudinal diversity gradient evident for trees^3,4. Theory suggests that CNDD based on pairwise interactions alone is not sufficient to explain the intricacies of diverse communities, because higher-order interactions (HOIs) may greatly modify these interactions^5,6. However, there has been a lack of empirical studies investigating how HOIs intertwine with pairwise interactions and how they may contribute to the latitudinal tree diversity gradient. Here we examined both pairwise interactions and HOIs across 32 large permanent forest plots, most in the northern hemisphere. We detected evidence of HOIs in 40% of the 1,543 species-plot combinations for tree growth, and 23% of the 1,340 such combinations for tree survival, with the strength of these interactions declining with latitude. HOIs were found to benefit rare species but disadvantage common species, suggesting a potential mechanism promoting species diversity. This stabilizing effect weakened towards higher latitudes, consistent with the latitudinal tree diversity gradient. Our findings reveal an important interplay between pairwise interactions and HOIs in promoting the latitudinal tree diversity gradient and help to clarify the contribution of CNDD to this biogeographic pattern.