Understanding the spatial distribution of tree species in subtropical evergreen-deciduous broadleaf karst forest is fundamental to studying species coexistence and karst species diversity. Here, complete spatial randomness and heterogeneous Poisson process models were used to analyze the spatial distribution patterns of 146 species with at least one individual per ha in a 25-ha plot in southwest China. We used canonical correspondence analysis (CCA) and the torus-translation test (TTT) to explain the distributions of observed species. Our results show that an aggregated distribution was the dominant pattern in Mulun karst forests; the percentage and intensity of aggregated decreased with increasing spatial scale, abundance, mean diameter at breast height (DBH), and maximum DBH. Rare species were more aggregated than intermediately abundant and abundant species. However, functional traits (e.g., growth form and phenological guild) had no significant effects on the distributions of species. The CCA revealed that the four analyzed topographic variables (elevation, slope, aspect, and convexity) had significant influences on species distributions. The TTT showed that not all species have habitat preferences and that 68.5% (100 out of 146 species) show a strongly positive or negative association with at least one habitat. Most species were inclined to grow on slopes and hilltops.