We report the genome of the facultative intracellular parasite Rhodococcus equi, the only animal pathogen within the biotechnologically important actinobacterial genus Rhodococcus. The 5.0-Mb R. equi 103S genome is significantly smaller than those of environmental rhodococci. This is due to genome expansion in nonpathogenic species, via a linear gain of paralogous genes and an accelerated genetic flux, rather than reductive evolution in R. equi. The 103S genome lacks the extensive catabolic and secondary metabolic complement of environmental rhodococci, and it displays unique adaptations for host colonization and competition in the short-chain fatty acid–rich intestine and manure of herbivores—two main R. equi reservoirs. Except for a few horizontally acquired (HGT) pathogenicity loci, including a cytoadhesive pilus determinant (rpl) and the virulence plasmid vap pathogenicity island (PAI) required for intramacrophage survival, most of the potential virulence-associated genes identified in R. equi are conserved in environmental rhodococci or have homologs in nonpathogenic Actinobacteria. This suggests a mechanism of virulence evolution based on the cooption of existing core actinobacterial traits, triggered by key host niche–adaptive HGT events. We tested this hypothesis by investigating R. equi virulence plasmid-chromosome crosstalk, by global transcription profiling and expression network analysis. Two chromosomal genes conserved in environmental rhodococci, encoding putative chorismate mutase and anthranilate synthase enzymes involved in aromatic amino acid biosynthesis, were strongly coregulated with vap PAI virulence genes and required for optimal proliferation in macrophages. The regulatory integration of chromosomal metabolic genes under the control of the HGT–acquired plasmid PAI is thus an important element in the cooptive virulence of R. equi.

dWe report a novel host-associated virulence plasmid in Rhodococcus equi, pVAPN, carried by bovine isolates of this facultative intracellular pathogenic actinomycete. Surprisingly, pVAPN is a 120-kb invertron-like linear replicon unrelated to the circular virulence plasmids associated with equine (pVAPA) and porcine (pVAPB variant) R. equi isolates. pVAPN is similar to the linear plasmid pNSL1 from Rhodococcus sp. NS1 and harbors six new vap multigene family members (vapN to vapS) in a vap pathogenicity locus presumably acquired via en bloc mobilization from a direct predecessor of equine pVAPA. Loss of pVAPN rendered R. equi avirulent in macrophages and mice. Mating experiments using an in vivo transconjugant selection strategy demonstrated that pVAPN transfer is sufficient to confer virulence to a plasmid-cured R. equi recipient. Phylogenetic analyses assigned the vap multigene family complement from pVAPN, pVAPA, and pVAPB to seven monophyletic clades, each containing plasmid type-specific allelic variants of a precursor vap gene carried by the nearest vap island ancestor. Deletion of vapN, the predicted "bovine-type" allelic counterpart of vapA, essential for virulence in pVAPA, abrogated pVAPN-mediated intramacrophage proliferation and virulence in mice. Our findings support a model in which R. equi virulence is conferred by host-adapted plasmids. Their central role is mediating intracellular proliferation in macrophages, promoted by a key vap determinant present in the common ancestor of the plasmid-specific vap islands, with host tropism as a secondary trait selected during coevolution with specific animal species.

No abstract

We report a comparative study of 29 representative genomes of the animal pathogen Rhodococcus equi. The analyses showed that R. equi is genetically homogeneous and clonal, with a large core genome accounting for ≈80% of an isolates’ gene content. An open pangenome, even distribution of accessory genes among the isolates, and absence of significant core–genome recombination, indicated that gene gain/loss is a main driver of R. equi genome evolution. Traits previously predicted to be important in R. equi physiology, virulence and niche adaptation were part of the core genome. This included the lack of a phosphoenolpyruvate:carbohydrate transport system (PTS), unique among the rhodococci except for the closely related Rhodococcus defluvii, reflecting selective PTS gene loss in the R. equi–R. defluvii sublineage. Thought to be asaccharolytic, rbsCB and glcP non-PTS sugar permease homologues were identified in the core genome and, albeit inefficiently, R. equi utilized their putative substrates, ribose and (irregularly) glucose. There was no correlation between R. equi whole-genome phylogeny and host or geographical source, with evidence of global spread of genomovars. The distribution of host-associated virulence plasmid types was consistent with the exchange of the plasmids (and corresponding host shifts) across the R. equi population, and human infection being zoonotically acquired. Phylogenomic analyses demonstrated that R. equi occupies a central position in the Rhodococcus phylogeny, not supporting the recently proposed transfer of the species to a new genus.

The conjugative virulence plasmid is a key component of the Rhodococcus equi accessory genome essential for pathogenesis. Three host-associated virulence plasmid types have been identified: the equine pVAPA and porcine pVAPB circular variants, and the linear pVAPN found in bovine (ruminant) isolates. We recently characterized the R. equi pangenome (Anastasi E, et al. 2016. Pangenome and phylogenomic analysis of the pathogenic actinobacterium Rhodococcus equi. Genome Biol Evol. 8:3140–3148.) and we report here the comparative analysis of the virulence plasmid genomes. Plasmids within each host-associated type were highly similar despite their diverse origins. Variation was accounted for by scattered single nucleotide polymorphisms and short nucleotide indels, while larger indels—mostly in the plasticity region near the vap pathogencity island (PAI)—defined plasmid genomic subtypes. Only one of the plasmids analyzed, of pVAPN type, was exceptionally divergent due to accumulation of indels in the housekeeping backbone. Each host-associated plasmid type carried a unique PAI differing in vap gene complement, suggesting animal host-specific evolution of the vap multigene family. Complete conservation of the vap PAI was observed within each host-associated plasmid type. Both diversity of host-associated plasmid types and clonality of specific chromosomal-plasmid genomic type combinations were observed within the same R. equi phylogenomic subclade. Our data indicate that the overall strong conservation of the R. equi host-associated virulence plasmids is the combined result of host-driven selection, lateral transfer between strains, and geographical spread due to international livestock exchanges.

This is the first molecular characterization of resistance to macrolides, lincosamides and streptogramins B in R. equi. Resistance was due to the presence of a novel erm(46) gene mobilizable likely by conjugation, which has spread among equine isolates of R. equi in the USA.

PagL and LpxO are enzymes that modify lipid A. PagL is a 3-O deacylase that removes the primary acyl chain from the 3 position, and LpxO is an oxygenase that 2-hydroxylates specific acyl chains in the lipid A. pagL and lpxO homologues have been identified in the genome of Bordetella bronchiseptica, but in the current structure for B. bronchiseptica lipid A the 3 position is acylated and 2-OH acylation is not reported. We have investigated the role of B. bronchiseptica pagL and lpxO in lipid A biosynthesis. We report a different structure for wild-type (WT) B. bronchiseptica lipid A, including the presence of 2-OH-myristate, the presence of which is dependent on lpxO. We also demonstrate that the 3 position is not acylated in the major WT lipid A structures but that mutation of pagL results in the presence of 3-OH-decanoic acid at this position, suggesting that lipid A containing this acylation is synthesized but that PagL removes most of it from the mature lipid A. These data refine the structure of B. bronchiseptica lipid A and demonstrate that pagL and lpxO are involved in its biosynthesis.

No abstract