Toxin⁻antitoxin (TA) systems are diverse genetic modules with demonstrated roles in plasmid stability, stress management, biofilm formation and antibiotic persistence. However, relatively little is known about their functional significance in plant pathogens. In this study we characterize type II and IV TA systems in the economically important plant pathogen Erwinia amylovora. Hidden Markov Model (HMM) and BLAST-based programs were used to predict the identity and distribution of putative TA systems among sequenced genomes of E. amylovora and other plant-associated Erwinia spp. Of six conserved TA systems tested for function from E. amylovora, three (CbtA/CbeA, ParE/RHH and Doc/PhD) were validated as functional. CbtA was toxic to E. amylovora, but not to Escherichia coli. While the E. coli homolog of CbtA elicits the formation of lemon-shaped cells upon overexpression and targets cytoskeletal proteins FtsZ and MreB, E. amylovora CbtA led to cell elongation and did not interact with these cytoskeletal proteins. Phylogenetic analysis revealed that E. amylovora CbtA belongs to a distinct clade from the CbtA of pathogenic E. coli. This study expands the repertoire of experimentally validated TA systems in plant pathogenic bacteria, and suggests that the E. amylovora homolog of CbtA is functionally distinct from that of E. coli.
Keywords: CbtA; Doc; Erwinia amylovora; ParE; YeeV; cell elongation; toxin–antitoxin system.