The structure and distribution patterns of Himalayan vegetation are poorly explored, and research on species composition along the elevation gradient in these mountain ranges is still deficient. The current study was undertaken to analyze the variation and pattern of plant species composition along a vertical gradient in northwestern Himalaya, India. A total of 18 sites were selected along an elevation gradient ranging from 2200 to 3900 m asl positioned at an interval of 100 m. The Renyi diversity profile, non-metric multidimensional scaling based on the Bray–Curtis dissimilarity metric and beta diversity components among the elevation belts were calculated. Furthermore, to study the influence of altitude on species richness and diversity, a generalized additive model was created. Two hundred and ten plant species representing 66 families and 147 genera were recorded. The Renyi diversity profiles show that the lower and mid-altitudes had rich species diversity. The results of the non-metric multidimensional scaling analysis show a considerable variation in the total plant species composition among the studied elevation belts. The observed multiple-site Sorensen dissimilarity index across the studied elevation belts was very high. The contribution of species replacement or the turnover component to the observed dissimilarity was much higher than the nestedness component. Furthermore, the herbaceous and tree richness showed a significant decrease with increase in elevation; however, the richness of shrubs showed a bimodal pattern. The present study increases our understanding of the trends and patterns of species richness along the vertical gradient in the Himalayan region. © 2022 by the authors.