Diversity and biogeography of plant phyllosphere bacteria are governed by latitude-dependent mechanisms.

Predicting and managing the structure and function of plant microbiomes requires quantitative understanding of community assembly and predictive models of spatial distributions at broad geographic scales. Here, we quantified the relative contribution of abiotic and biotic factors to the assembly of phyllosphere bacterial communities, and developed spatial distribution models for keystone bacterial taxa along a latitudinal gradient, by analyzing 16S rRNA gene sequences from 1453 leaf samples taken from 329 plant species in China. We demonstrated a latitudinal gradient in phyllosphere bacterial diversity and community composition, which was mostly explained by climate and host plant factors. We found that host-related factors were increasingly important in explaining bacterial assembly at higher latitudes while nonhost factors including abiotic environments, spatial proximity and plant neighbors were more important at lower latitudes. We further showed that local plant-bacteria associations were interconnected by hub bacteria taxa to form metacommunity-level networks, and the spatial distribution of these hub taxa was controlled by hosts and spatial factors with varying importance across latitudes. For the first time, we documented a latitude-dependent importance in the driving factors of phyllosphere bacteria assembly and distribution, serving as a baseline for predicting future changes in plant phyllosphere microbiomes under global change and human activities.