Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis

Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis, whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis. We found 590 amplifications in M. tuberculosis, and like B. pertussis, these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis. This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis, highlighting the need for a more holistic understanding of bacterial genetics.

Fundamental differences in physiology of Bordetella pertussis dependent on the two-component system Bvg revealed by gene essentiality studies

The identification of genes essential for a bacterium's growth reveals much about its basic physiology under different conditions. Bordetella pertussis, the causative agent of whooping cough, adopts both virulent and avirulent states through the activity of the two-component system, Bvg. The genes essential for B. pertussis growth in vitro were defined using transposon sequencing, for different Bvg-determined growth states. In addition, comparison of the insertion indices of each gene between Bvg phases identified those genes whose mutation exerted a significantly different fitness cost between phases. As expected, many of the genes identified as essential for growth in other bacteria were also essential for B. pertussis. However, the essentiality of some genes was dependent on Bvg. In particular, a number of key cell wall biosynthesis genes, including the entire mre/mrd locus, were essential for growth of the avirulent (Bvg minus) phase but not the virulent (Bvg plus) phase. In addition, cell wall biosynthesis was identified as a fundamental process that when disrupted produced greater fitness costs for the Bvg minus phase compared to the Bvg plus phase. Bvg minus phase growth was more susceptible than Bvg plus phase growth to the cell wall-disrupting antibiotic ampicillin, demonstrating the increased susceptibility of the Bvg minus phase to disruption of cell wall synthesis. This Bvg-dependent conditional essentiality was not due to Bvg-regulation of expression of cell wall biosynthesis genes; suggesting that this fundamental process differs between the Bvg phases in B. pertussis and is more susceptible to disruption in the Bvg minus phase. The ability of a bacterium to modify its cell wall synthesis is important when considering the action of antibiotics, particularly if developing novel drugs targeting cell wall synthesis.

A qPCR assay for Bordetella pertussis cells that enumerates both live and dead bacteria

Bordetella pertussis is the causative agent of whooping cough, commonly referred to as pertussis. Although the incidence of pertussis was reduced through vaccination, during the last thirty years it has returned to high levels in a number of countries. This resurgence has been linked to the switch from the use of whole-cell to acellular vaccines. Protection afforded by acellular vaccines appears to be short-lived compared to that afforded by whole cell vaccines. In order to inform future vaccine improvement by identifying immune correlates of protection, a human challenge model of B. pertussis colonisation has been developed. Accurate measurement of colonisation status in this model has required development of a qPCR-based assay to enumerate B. pertussis in samples that distinguishes between viable and dead bacteria. Here we report the development of this assay and its performance in the quantification of B. pertussis from human challenge model samples. This assay has future utility in diagnostic labs and in research where a quantitative measure of both B. pertussis number and viability is required.

The evolution of Bordetella pertussis has selected for mutations of acr that lead to sensitivity to hydrophobic molecules and fatty acids

Whooping cough, or pertussis, is resurgent in numerous countries worldwide. This has renewed interest in Bordetella pertussis biology and vaccinology. The in vitro growth of B. pertussis has been a source of difficulty, both for the study of the organism and the production of pertussis vaccines. It is inhibited by fatty acids and other hydrophobic molecules. The AcrAB efflux system is present in many different bacteria and in combination with an outer membrane factor exports acriflavine and other small hydrophobic molecules from the cell. Here, we identify that the speciation of B. pertussis has selected for an Acr system that is naturally mutated and displays reduced activity compared to B. bronchiseptica, in which the system appears intact. Replacement of the B. pertussis locus with that of B. bronchiseptica conferred higher levels of resistance to growth inhibition by acriflavine and fatty acids. In addition, we identified that the transcription of the locus is repressed by a LysR-type transcriptional regulator. Palmitate de-represses the expression of the acr locus, dependent on the LysR regulator, strongly suggesting that it is a transcriptional repressor that is regulated by palmitate. It is intriguing that the speciation of B. pertussis has selected for a reduction in activity of the Acr efflux system that typically is regarded as protective to bacteria.

Comparative Genomics of Rhodococcus equi Virulence Plasmids Indicates Host-Driven Evolution of the vap Pathogenicity Island

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.

Pangenome and Phylogenomic Analysis of the Pathogenic Actinobacterium Rhodococcus equi

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.

Novel transferable erm(46) determinant responsible for emerging macrolide resistance in Rhodococcus equi

OBJECTIVES

The objective of this study was to identify the molecular mechanism of macrolide resistance in the actinomycete Rhodococcus equi, a major equine pathogen and zoonotic agent causing opportunistic infections in people.

METHODS

Macrolide-resistant (n = 62) and macrolide-susceptible (n = 62) clinical isolates of R. equi from foals in the USA were studied. WGS of 18 macrolide-resistant and 6 macrolide-susceptible R. equi was performed. Representative sequences of all known macrolide resistance genes identified to date were used to search the genome assemblies for putative homologues. PCR was used to screen for the presence of the identified resistance determinant in the rest of the isolates. Mating experiments were performed to verify mobility of the gene.

RESULTS

A novel erm gene, erm(46), was identified in all sequenced resistant isolates, but not in susceptible isolates. There was complete association between macrolide resistance and the presence of erm(46) as detected by PCR screening of all 124 clinical isolates of R. equi. Expression of erm(46) in a macrolide-susceptible strain of R. equi induced high-level resistance to macrolides, lincosamides and streptogramins B, but not to other classes of antimicrobial agents. Transfer of erm(46) to macrolide-susceptible R. equi was confirmed. The transfer frequency ranged from 3 × 10(-3) to 1 × 10(-2).

CONCLUSIONS

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.

Mouse lung infection model to assess Rhodococcus equi virulence and vaccine protection

The pathogenic actinomycete Rhodococcus equi causes severe purulent lung infections in foals and immunocompromised people. Although relatively unsusceptible to R. equi, mice are widely used for in vivo studies with this pathogen. The most commonly employed mouse model is based on systemic (intravenous) infection and determination of R. equi burdens in spleen and liver. Here, we investigated the murine lung for experimental infection studies with R. equi. Using a 10(7)CFU intranasal challenge in BALB/c mice, virulent R. equi consistently survived in quantifiable numbers up to 10 days in the lungs whereas virulence-deficient R. equi bacteria were rapidly cleared. An internally controlled virulence assay was developed in which the test R. equi strains are co-inoculated and monitored in the same mouse. Isogenic R. equi bacteria lacking either the plasmid vapA gene or the entire virulence plasmid were compared using this competitive assay. Both strains showed no significant differences in in vivo fitness in the lung, indicating that the single loss of the virulence factor VapA was sufficient to account for the full attenuation seen in the absence of the virulence plasmid. To test the adequacy of the lung infection model for monitoring R. equi vaccine efficacy, BALB/c mice were immunized with live R. equi and challenged intranasally. Vaccination conferred protection against acute pulmonary challenge with virulent R. equi. Our data indicate that the murine lung infection model provides a useful tool for both R. equi virulence and vaccine studies.

Rhodococcus equi: the many facets of a pathogenic actinomycete

Rhodococcus equi is a soil-dwelling pathogenic actinomycete that causes pulmonary and extrapulmonary pyogranulomatous infections in a variety of animal species and people. Young foals are particularly susceptible and develop a life-threatening pneumonic disease that is endemic at many horse-breeding farms worldwide. R. equi is a facultative intracellular parasite of macrophages that replicates within a modified phagocytic vacuole. Its pathogenicity depends on a virulence plasmid that promotes intracellular survival by preventing phagosome-lysosome fusion. Species-specific tropism of R. equi for horses, pigs and cattle appears to be determined by host-adapted virulence plasmid types. Molecular epidemiological studies of these plasmids suggest that human R. equi infection is zoonotic. Analysis of the recently determined R. equi genome sequence has identified additional virulence determinants on the bacterial chromosome. This review summarizes our current understanding of the clinical aspects, biology, pathogenesis and immunity of this fascinating microbe with plasmid-governed infectivity.

The genome of a pathogenic rhodococcus: cooptive virulence underpinned by key gene acquisitions

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.

Rhodococcus is a prototypic genus within the Actinobacteria, one of the largest microbial groups on Earth. Many of the ubiquitous rhodococcal species are biotechnologically useful due to their metabolic versatility and biodegradative properties. We have deciphered the genome of a facultatively parasitic Rhodococcus, the animal and human pathogen R. equi. Comparative genomic analyses of related species provide a unique opportunity to increase our understanding of niche-adaptive genome evolution and specialization. The environmental rhodococci have much larger genomes, richer in metabolic and degradative pathways, due to gene duplication and acquisition, not genome contraction in R. equi. This probably reflects that the host-associated R. equi habitat is more stable and favorable than the chemically diverse but nutrient-poor environmental niches of nonpathogenic rhodococci, necessitating metabolically more complex, expanded genomes. Our work also highlights that the recruitment or cooption of core microbial traits, following the horizontal acquistion of a few critical genes that provide access to the host niche, is an important mechanism in actinobacterial virulence evolution. Gene cooption is a key evolutionary mechanism allowing rapid adaptive change and novel trait acquisition. Recognizing the contribution of cooption to virulence provides a rational framework for understanding and interpreting the emergence and evolution of microbial pathogenicity.

IEIIS Meeting minireview: Bordetella evolution: lipid A and Toll-like receptor 4

The evolution of Bordetella pertussis and Bordetella parapertussis from Bordetella bronchiseptica involved changes in host range and pathogenicity. Recent data suggest that the human-adapted Bordetella modified their interaction with host immune systems to effect these changes and that decreased stimulation of Toll-like receptor 4 (TLR4) by lipid A is central to this. We discuss Bordetella lipid A structure and genetics within the context of evolution and host immunity.