Dynamics of bacterial community abundance and structure of a newly established horizontal subsurface flow (HSSF) pilot-scale wetland were studied using high-throughput sequencing and quantitative polymerase chain reaction (PCR) methods. Bacterial community abundance increased rapidly within one month and stabilised thereafter in three replicate HSSF constructed wetland (CW) mesocosms. The most dominant phylum was Proteobacteria, followed by Bacteroidetes in wetland media biofilms and Firmicutes in influent wastewater. CW bacterial community diversity increased over time and was positively related to the wastewater treatment efficiency. Increase in the abundance of total bacteria in the community was accompanied with the abundance of denitrifying bacteria that promoted nitrate and nitrite removal from the wastewater. During the 150-day study period, similar patterns of bacterial community successions were observed in replicate HSSF CW mesocosms. The data indicate that successions in the bacterial community in HSSF CW are shaped by biotic interactions, with a significant contribution made by external abiotic factors such as influent chemical parameters. Network analysis of the bacterial community revealed that organic matter and nitrogen removal in HSSF CW could be, in large part, allocated to a small subset of tightly interconnected bacterial species. The diversity of bacterial community and abundance of denitrifiers were good predictors of the removal efficiency of ammonia, nitrate and total organic C in HSSF CW mesocosms, while the removal of the seven-day biochemical oxygen demand (BOD 7 ) was best predicted by the abundance of a small set of bacterial phylotypes. The results suggest that nitrogen removal in HSSF CW consist of two main pathways. The first is heterotrophic nitrification, which is coupled with aerobic denitrification and mediated by mixotrophic nitrite-oxidizers. The second pathway is anaerobic denitrification, which leads to gaseous intermediates and loss of nitrogen as N 2 .