Comparative genomics unravels mechanisms of genetic adaptation for the catabolism of the phenylurea herbicide linuron inVariovorax

Biodegradation of the phenylurea herbicide linuron appears a specialization within a specific clade of the Variovorax genus. The linuron catabolic ability is likely acquired by horizontal gene transfer but the mechanisms involved are not known. The full genome sequences of six linuron degrading Variovorax strains isolated from geographically distant locations were analyzed to acquire insight in the mechanisms of genetic adaptation towards linuron metabolism in Variovorax . Whole genome sequence analysis confirmed the phylogenetic position of the linuron degraders in a separate clade within Variovorax and indicated their unlikely origin from a common ancestral linuron degrader. The linuron degraders differentiated from non-degraders by the presence of multiple plasmids of 20 to 839 kb, including plasmids of unknown plasmid groups. The linuron catabolic gene clusters showed (i) high conservation and synteny and (ii) strain-dependent distribution among the different plasmids. All were bordered by IS 1071 elements forming composite transposon structures appointing IS 1071 as key for catabolic gene recruitment. Most of the strain carried at least one broad host range plasmid that might have been a second instrument for catabolic gene acquisition. We conclude that clade 1 Variovorax strains, despite their different geographical origin, made use of a limited genetic repertoire to acquire linuron biodegradation. Importance The genus Variovorax and especially a clade of strains that phylogenetically separates from the majority of Variovorax species, appears to be a specialist in the biodegradation of the phenyl urea herbicide linuron. Horizontal gene transfer (HGT) likely played an essential role in the genetic adaptation of those strain to acquire the linuron catabolic genotype. However, we do not know the genetic repertoire involved in this adaptation both regarding catabolic gene functions as well as gene functions that promote HGT neither do we know how this varies between the different strains. These questions are addressed in this paper by analyzing the full genome sequences of six linuron degrading Variovorax strains. This knowledge is important for understanding the mechanisms that steer world-wide genetic adaptation in a particular species and this for a particular phenotypic trait as linuron biodegradation.