To understand the relationship between elevation and bacterial communities in wastewater treatment plants (WWTPs), bacterial communities in 21 municipal WWTPs across China, located 9 to 3,660 m above sea level (masl), were investigated by 454 pyrosequencing. A threshold for the association of elevation with bacterial community richness and evenness was observed at approximately 1,200 masl. At lower elevations, both richness and evenness were not significantly associated with elevation. At higher elevations, significant declines with increased elevations were observed for community richness and evenness. The declining evenness trend at the phylum level was reflected by distinct trends in relative abundance for individual bacterial phyla. Betaproteobacteria, Bacteroidetes, and Firmicutes displayed significant increases, while most other phyla showed declines. Spearman correlation analysis indicated that the community richness and evenness at high elevations were more correlated with elevation than with any other single environmental variable. Redundancy analysis indicated that the contribution of elevation to community composition variances increased from 3% at lower elevations to 11% at higher elevations whereas the community composition variance at higher elevations remained much more explained by operational variables (39.2%) than by elevation. The influent total phosphorus concentration, food/microorganism ratio, and treatment process were the three shared dominant contributors to the community composition variance across the whole elevation gradient, followed by effluent ammonia nitrogen and temperature at higher elevations. W astewater treatment plants (WWTPs) are broadly applied for the efficient removal of contaminants from wastewater before its discharge into the natural environment. A complex microbial community is the most dominant contributor to the removal of contaminants in WWTPs. A better understanding of the microbial community structure and its relationship to influential variables could provide important guidance for the improvement of wastewater treatment efficiency (1). Multiple studies have investigated the significant associations between WWTP microbial community structures and operational parameters or wastewater characteristics, such as dissolved oxygen concentration (2), temperature (3), treatment process (4), influent ammonia nitrogen concentration (5), biological oxygen demand (6), and total phosphorus (TP) concentration (7). However, regardless of the weakly neutral assembly of bacterial communities (8), there is still a large part of the bacterial community variance that cannot be explained by these well-known variables. Thus, it is presumed that there are still other environmental variables that potentially influence the variance of WWTP microbial communities.In recent years, a few studies were performed to explore the relationship between microbial communities and geographical locations of WWTPs. Zhang et al. (9) found significant geographical differences in bacterial communities amo...