Species data of 249 National Nature Reserves in China was used to identify potential underlying drivers of latitudinal gradients in plant diversity. We used generalized linear models (GLMs) to assess the correlations between predictor and response variables. We also used SAM (Spatial Analysis in Macroecology) to eliminate autocorrelation along each of the 249 studied locations. We used the Akaike information criterion (AICc; Montoya et al. 2007) to select the independent variables were those included in the best models from different combinations of climate, habitat and animal variables. Variance partitioning was used to decompose the variation in plant richness across different taxonomic levels among the three groups of predictors. We found that: Total plant species, gymnosperms, angiosperms and ferns showed significant latitudinal trends in richness (p < 0.001). Water-energy and habitat variables generally explained more variation in richness across different plant groups than did animal richness. Annual precipitation was selected as the best water-energy variable across different taxonomic plants groups, soil PH and elevation range were selected as the best habitat variables across different taxonomic plant groups. The independent effects of habitatvariables were higher than that of water-energy and animal variables across different taxonomic plants groups. Water-energy, habitat heterogeneity, and animal variables explain 48.8% of the variation in total species richness, 28.2% in gymnosperm richness, 44.2% in angiosperm richness, and 38.9% in fern richness.Plants showed significant latitudinal trends in richness (p < 0.001). Water-energy and habitat variables generally explained more variation in richness across different taxonomic plants groups than did animal variables. The independent effects of habitat variables were higher than those of water-energy and animal variables across different taxonomic plants groups.