Abstract
Phytoremediation using willows is thought to be a sustainable alternative to traditional remediation techniques involving excavation, transport, and landfilling. However, the complexity of the interaction between the willow and its associated highly diverse microbial communities makes the optimization of phytoremediation very difficult. Here, we have sequenced the rhizosphere metatranscriptome of four willow species and the plant root metatranscriptome for two willow species growing in petroleum hydrocarbon-contaminated and non-contaminated soils on a former petroleum refinery site. Significant differences in the abundance of transcripts related to different bacterial and fungal taxa were observed between willow species, mostly in contaminated soils. When comparing transcript abundance in contaminated vs. non-contaminated soil for each willow species individually, transcripts for many microbial taxa and functions were significantly more abundant in contaminated rhizosphere soil for Salix eriocephala, S. miyabeana and S. purpurea, in contrast to what was observed in the rhizosphere of S. caprea. This agrees with the previously reported sensitivity of S. caprea to contamination, and the superior tolerance of S. miyabeana and S. purpurea to soil contamination at that site. The root metatranscriptomes of two species were compared and revealed that plants transcripts are mainly influenced by willow species, while microbial transcripts mainly responded to contamination. A comparison of the rhizosphere and root metatranscriptomes in the S. purpurea species revealed a complete reorganization of the linkages between root and rhizosphere pathways when comparing willows growing in contaminated and non-contaminated soils, mainly because of large shifts in the rhizosphere metatranscriptome.