Summary: 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.

Keystone species are more important than others for community dynamics and stability. Keystone species can be identified and evaluated by their centrality (i.e., a relative ranking of the topological positional importance of a species) in ecological networks. Studies of node centrality of plant–fungus bipartite networks, for example, have identified the keystone species that are important for maintaining network structure and stability. However, the underlying drivers of the importance of species in a network have rarely been examined. We assessed the centrality (degree, closeness, and betweenness) of plant and fungal species in a plant–ectomycorrhizal fungus network in a subtropical forest in southern China. Based on the phylogenies of plants and fungi and plant traits, we explored ecological factors that led to a species taking a central position or not. We found one plant species (Ternstroemia gymnanthera) and four species of ectomycorrhizal fungi (Russula citrina, Scleroderma sp., and two Cenococcum sp.) were characterized by the highest centrality of degree, closeness, and betweenness among the bipartite network nodes and thus played key roles in maintaining network structure. Centrality for fungi (not for plants) was phylogenetically constrained. Plant traits and abundance together explained 46.36%, 46.0%, and 43.7% of variation in the centrality of degree, closeness, and betweenness of plant species in the bipartite network, respectively. When plant or fungal species were sequentially removed on the order of higher to lower centrality, network was less stable than randomly removed. We suggest that abundance and traits determine the positional importance of plant species in a network. This work helps understand how plant–fungus association networks will respond to species extinction and changes in species abundance and functional traits due to habitat fragmentation and human activities.

Host specialization plays a critical role in the ecology and evolution of plant–microbe symbiosis. Theory predicts that host specialization is associated with microbial genome streamlining and is influenced by the abundance of host species, both of which can vary across latitudes, leading to a latitudinal gradient in host specificity. Here, we quantified the host specificity and composition of plant–bacteria symbioses on leaves across 329 tree species spanning a latitudinal gradient. Our analysis revealed a predominance of host‐specialized leaf bacteria. The degree of host specificity was negatively correlated with bacterial genome size and the local abundance of host plants. Additionally, we found an increased host specificity at lower latitudes, aligning with the high prevalence of small bacterial genomes and rare host species in the tropics. These findings underscore the importance of genome streamlining and host abundance in the evolution of host specificity in plant‐associated bacteria along the latitudinal gradient.

Pulmonary dysbiosis may predispose people living with HIV (PLWH) to chronic lung disease. Herein, we assessed whether intrapulmonary HIV reservoir size and immune disruption are associated with reduced bacterial lung diversity in PLWH. Bacterial DNA was extracted and PCR-amplified from cell-free bronchoalveolar lavage (BAL) fluid from 28 PLWH and 9 HIV-negative controls. Amplicon sequence variant (ASV) relative abundances and taxonomic identities were analyzed using joint species distribution modeling. HIV-DNA was quantified from blood and pulmonary CD4+ T-cells using ultra-sensitive qPCR. Immunophenotyping of BAL T-cells was performed using flow cytometry. Lung microbiome diversity was lower in smokers than non-smokers and microbiome composition was more variable in PLWH than HIV-negative individuals. Frequencies of effector memory BAL CD4+ and CD8+ T-cells positively correlated with abundance of several bacterial families while frequencies of BAL activated CD4+ T-cells negatively correlated with abundance of most lung bacterial families. Higher HIV-DNA levels in blood, but not in BAL, as well as frequencies of senescent CD4+ T-cells were associated with reduced bacterial diversity. These findings suggest that HIV infection may weaken the relationship between the lung microbiome and smoking status. Viral reservoir and immune activation levels may impact the lung microbiome, predisposing PLWH to pulmonary comorbidities.