Greenhouse gas (GHG) emissions from constructed wetlands (CWs) have raised environmental concern and thus offset their environmental and ecological benefits. This study evaluated the influence of plant species, i.e., Canna indica (C. indica), Cyperus alternifolius (C. alternifolius), Phragmites australis (P. australis) and unplanted control, on GHG emissions, pollutant removal and associated microbial abundance in subsurface flow constructed wetland (SSFCW) mesocosms. C. indica outperformed the other tested plant species in pollutant removal, and the presence of plants irrespective of species enhanced the removal efficiencies of nitrogen, phosphorus and organics in SSFCW mesocosms compared to unplanted control. The greatest carbon dioxide (CO2) flux (582.01 +/- 89.25 mg/m(2)/h), methane (CH4) flux (21.88 +/- 2.51 mu g/m(2)/h) and nitrous oxide (N2O) flux (37.27 +/- 15.82 mu g/m(2)/h) were observed in mesocosms planted with C. indica, P. australis and Calternifolius, respectively. Unexpectedly, the mcrA and pmoA genes were not detected in any mesocosms. For denitrifiers, the N2O fluxes showed a significantly (p < 0.05) positive correlation with nirS and nirK genes abundance. The abundance of nosZ gene (ranged from 0.18 x 10(4) to 0.75 x 10(4) copies/mg gravel) and nosZ/(nirS + nirK) (ranged from 1.29 x 10(-4) to 2.12 x 10(-4) copies/mg gravel) in this study was lower than that in most reported studies. Regarding the global warming potential (GWP), the lowest value was observed in mesocosms planted with C. indica. In conclusion, C. indica is selected as the optimal plant species in this study due to its lower GWP and excellent pollutant removal performance. (C) 2019 Elsevier Ltd. All rights reserved.