Bacterial communities in the activated sludge (AS) determine the wastewater treatment performance in the municipal wastewater treatment plants (WWTPs). Aiming at identifying the affecting factors and the variation patterns of the bacterial assemblages in AS, a 2-year time-series AS samples were collected from two separated WWTPs and metagenomic sequencing was conducted. Obvious seasonal shift and succession of the bacterial community were observed in both WWTPs on the genus and species levels, especially for the persistent taxa, implying that temperature was a decisive factor for maintaining bacterial assemblage patterns in long-term period. Taxa abundance distribution (TAD) concerning occurrence frequency and average abundance were found fitting for exponential formulations, and the approximately equal total abundance of persistent taxa suggested that stable and high abundance (~ 90%) of core functional bacterial groups would help to maintain wastewater treatment performance. Drastic changes of environmental factors were found causing temporally significant bacterial structure variation, while the innate correlations between bacterial species could recover the community gradually and maintain relative stability of the AS system. Delayed correlations between environmental factors and abundant (persistent or intermittent) bacterial species were observed widely, while synchronous biotic interactions were identified more frequently. Besides, bacterial species with similar functions were prone to cluster together and shared the same seasonal variation pattern, implicating that the cooperation of functional correlated taxa played the most dominant role in shaping the bacterial assemblages. Furthermore, rare bacterial groups were to be explored for removing emerging pollutants with lower concentrations. The results of this study would assist dealing with operational defect and optimize the treatment system in WWTPs.