Genomics of microbial plasmids: classification and identification based on replication and transfer systems and host taxonomy.

Plasmids are key vectors of horizontal gene transfer and essential genetic engineering tools. They code for genes involved in many aspects of microbial biology, including detoxication, virulence, ecological interactions, and antibiotic resistance. While many studies have decorticated the mechanisms of mobility in model plasmids, the identification and characterization of plasmid mobility from genome data are unexplored. By reviewing the available data and literature, we established a computational protocol to identify and classify conjugation and mobilization genetic modules in 1,730 plasmids. This allowed the accurate classification of proteobacterial conjugative or mobilizable systems in a combination of four mating pair formation and six relaxase families. The available evidence suggests that half of the plasmids are nonmobilizable and that half of the remaining plasmids are conjugative. Some conjugative systems are much more abundant than others and preferably associated with some clades or plasmid sizes. Most very large plasmids are nonmobilizable, with evidence of ongoing domestication into secondary chromosomes. The evolution of conjugation elements shows ancient divergence between mobility systems, with relaxases and type IV coupling proteins (T4CPs) often following separate paths from type IV secretion systems. Phylogenetic patterns of mobility proteins are consistent with the phylogeny of the host prokaryotes, suggesting that plasmid mobility is in general circumscribed within large clades. Our survey suggests the existence of unsuspected new relaxases in archaea and new conjugation systems in cyanobacteria and actinobacteria. Few genes, e.g., T4CPs, relaxases, and VirB4, are at the core of plasmid conjugation, and together with accessory genes, they have evolved into specific systems adapted to specific physiological and ecological contexts.

The genome of the bacterium Borrelia burgdorferi B31, the aetiologic agent of Lyme disease, contains a linear chromosome of 910,725 base pairs and at least 17 linear and circular plasmids with a combined size of more than 533,000 base pairs. The chromosome contains 853 genes encoding a basic set of proteins for DNA replication, transcription, translation, solute transport and energy metabolism, but, like Mycoplasma genitalium, it contains no genes for cellular biosynthetic reactions. Because B. burgdorferi and M. genitalium are distantly related eubacteria, we suggest that their limited metabolic capacities reflect convergent evolution by gene loss from more metabolically competent progenitors. Of 430 genes on 11 plasmids, most have no known biological function; 39% of plasmid genes are paralogues that form 47 gene families. The biological significance of the multiple plasmid-encoded genes is not clear, although they may be involved in antigenic variation or immune evasion.