Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

A small predatory core genome in the divergent marine Bacteriovorax marinus SJ and the terrestrial Bdellovibrio bacteriovorus

Bacteriovorax marinus SJ is a predatory delta-proteobacterium isolated from a marine environment. The genome sequence of this strain provides an interesting contrast to that of the terrestrial predatory bacterium Bdellovibrio bacteriovorus HD100. Based on their predatory lifestyle, Bacteriovorax were originally designated as members of the genus Bdellovibrio but subsequently were re-assigned to a new genus and family based on genetic and phenotypic differences. B. marinus attaches to gram-negative bacteria, penetrates through the cell wall to form a bdelloplast, in which it replicates, as shown using microscopy. Bacteriovorax is distinct, as it shares only 30% of its gene products with its closest sequenced relatives. Remarkably, 34% of predicted genes over 500 nt in length were completely unique with no significant matches in the databases. As expected, Bacteriovorax shares several characteristic loci with the other delta-proteobacteria. A geneset shared between Bacteriovorax and Bdellovibrio that is not conserved among other delta-proteobacteria such as Myxobacteria (which destroy prey bacteria externally via lysis), or the non-predatory Desulfo-bacteria and Geobacter species was identified. These 291 gene orthologues common to both Bacteriovorax and Bdellovibrio may be the key indicators of host-interaction predatory-specific processes required for prey entry. The locus from Bdellovibrio bacteriovorus is implicated in the switch from predatory to prey/host-independent growth. Although the locus is conserved in B. marinus, the sequence has only limited similarity. The results of this study advance understanding of both the similarities and differences between Bdellovibrio and Bacteriovorax and confirm the distant relationship between the two and their separation into different families.

Characterizing the flagellar filament and the role of motility in bacterial prey-penetration by Bdellovibrio bacteriovorus

The predatory bacterium Bdellovibrio bacteriovorus swims rapidly by rotation of a single, polar flagellum comprised of a helical filament of flagellin monomers, contained within a membrane sheath and powered by a basal motor complex. Bdellovibrio collides with, enters and replicates within bacterial prey, a process previously suggested to firstly require flagellar motility and then flagellar shedding upon prey entry. Here we show that flagella are not always shed upon prey entry and we study the six fliC flagellin genes of B. bacteriovorus, finding them all conserved and expressed in genome strain HD100 and the widely studied lab strain 109J. Individual inactivation of five of the fliC genes gave mutant Bdellovibrio that still made flagella, and which were motile and predatory. Inactivation of the sixth fliC gene abolished normal flagellar synthesis and motility, but a disordered flagellar sheath was still seen. We find that this non-motile mutant was still able to predate when directly applied to lawns of YFP-labelled prey bacteria, showing that flagellar motility is not essential for prey entry but important for efficient encounters with prey in liquid environments.

Within-host evolution of Burkholderia pseudomallei over a twelve-year chronic carriage infection

UNLABELLED: Burkholderia pseudomallei causes the potentially fatal disease melioidosis. It is generally accepted that B. pseudomallei is a noncommensal bacterium and that any culture-positive clinical specimen denotes disease requiring treatment. Over a 23-year study of melioidosis cases in Darwin, Australia, just one patient from 707 survivors has developed persistent asymptomatic B. pseudomallei carriage. To better understand the mechanisms behind this unique scenario, we performed whole-genome analysis of two strains isolated 139 months apart. During this period, B. pseudomallei underwent several adaptive changes. Of 23 point mutations, 78% were nonsynonymous and 43% were predicted to be deleterious to gene function, demonstrating a strong propensity for positive selection. Notably, a nonsense mutation inactivated the universal stress response sigma factor RpoS, with pleiotropic implications. The genome underwent substantial reduction, with four deletions in chromosome 2 resulting in the loss of 221 genes. The deleted loci included genes involved in secondary metabolism, environmental survival, and pathogenesis. Of 14 indels, 11 occurred in coding regions and 9 resulted in frameshift mutations that dramatically affected predicted gene products. Disproportionately, four indels affected lipopolysaccharide biosynthesis and modification. Finally, we identified a frameshift mutation in both P314 isolates within wcbR, an important component of the capsular polysaccharide I locus, suggesting virulence attenuation early in infection. Our study illustrates a unique clinical case that contrasts a high-consequence infectious agent with a long-term commensal infection and provides further insights into bacterial evolution within the human host.

IMPORTANCE: Some bacterial pathogens establish long-term infections that are difficult or impossible to eradicate with current treatments. Rapid advances in genome sequencing technologies provide a powerful tool for understanding bacterial persistence within the human host. Burkholderia pseudomallei is considered a highly pathogenic bacterium because infection is commonly fatal. Here, we document within-host evolution of B. pseudomallei in a unique case of human infection with ongoing chronic carriage. Genomic comparison of isolates obtained 139 months (11.5 years) apart showed a strong signal of adaptation within the human host, including inactivation of virulence and immunogenic factors, and deletion of pathways involved in environmental survival. Two global regulatory genes were mutated in the 139-month isolate, indicating extensive regulatory changes favoring bacterial persistence. Our study provides insights into B. pseudomallei pathogenesis and, more broadly, identifies parallel evolutionary mechanisms that underlie chronic persistence of all bacterial pathogens.

Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor

BACKGROUND: Klebsiella pneumoniae strains are pathogenic to animals and humans, in which they are both a frequent cause of nosocomial infections and a re-emerging cause of severe community-acquired infections. K. pneumoniae isolates of the capsular serotype K2 are among the most virulent. In order to identify novel putative virulence factors that may account for the severity of K2 infections, the genome sequence of the K2 reference strain Kp52.145 was determined and compared to two K1 and K2 strains of low virulence and to the reference strains MGH 78578 and NTUH-K2044.

RESULTS: In addition to diverse functions related to host colonization and virulence encoded in genomic regions common to the four strains, four genomic islands specific for Kp52.145 were identified. These regions encoded genes for the synthesis of colibactin toxin, a putative cytotoxin outer membrane protein, secretion systems, nucleases and eukaryotic-like proteins. In addition, an insertion within a type VI secretion system locus included sel1 domain containing proteins and a phospholipase D family protein (PLD1). The pld1 mutant was avirulent in a pneumonia model in mouse. The pld1 mRNA was expressed in vivo and the pld1 gene was associated with K. pneumoniae isolates from severe infections. Analysis of lipid composition of a defective E. coli strain complemented with pld1 suggests an involvement of PLD1 in cardiolipin metabolism.

CONCLUSIONS: Determination of the complete genome of the K2 reference strain identified several genomic islands comprising putative elements of pathogenicity. The role of PLD1 in pathogenesis was demonstrated for the first time and suggests that lipid metabolism is a novel virulence mechanism of K. pneumoniae.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

New ways to break an old bond: the bacterial carbon-phosphorus hydrolases and their role in biogeochemical phosphorus cycling

Phosphonates are organophosphorus molecules that contain the highly stable C-P bond, rather than the more common, and more labile, C-O-P phosphate ester bond. They have ancient origins but their biosynthesis is widespread among more primitive organisms and their importance in the contemporary biosphere is increasingly recognized; for example phosphonate-P is believed to play a particularly significant role in the productivity of the oceans. The microbial degradation of phosphonates was originally thought to occur only under conditions of phosphate limitation, mediated exclusively by the poorly characterized C-P lyase multienzyme system, under Pho regulon control. However, more recent studies have demonstrated the Pho-independent mineralization by environmental bacteria of three of the most widely distributed biogenic phosphonates: 2-aminoethylphosphonic acid (ciliatine), phosphonoacetic acid, and 2-amino-3-phosphonopropionic acid (phosphonoalanine). The three phosphonohydrolases responsible have unique specificities and are members of separate enzyme superfamilies; their expression is regulated by distinct members of the LysR family of bacterial transcriptional regulators, for each of which the phosphonate substrate of the respective degradative operon serves as coinducer. Previously no organophosphorus compound was known to induce the enzymes required for its own degradation. Whole-genome and metagenome sequence analysis indicates that the genes encoding these newly described C-P hydrolases are distributed widely among prokaryotes. As they are able to function under conditions in which C-P lyases are inactive, the three enzymes may play a hitherto-unrecognized role in phosphonate breakdown in the environment and hence make a significant contribution to global biogeochemical P-cycling.

Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita

Plant-parasitic nematodes are major agricultural pests worldwide and novel approaches to control them are sorely needed. We report the draft genome sequence of the root-knot nematode Meloidogyne incognita, a biotrophic parasite of many crops, including tomato, cotton and coffee. Most of the assembled sequence of this asexually reproducing nematode, totaling 86 Mb, exists in pairs of homologous but divergent segments. This suggests that ancient allelic regions in M. incognita are evolving toward effective haploidy, permitting new mechanisms of adaptation. The number and diversity of plant cell wall-degrading enzymes in M. incognita is unprecedented in any animal for which a genome sequence is available, and may derive from multiple horizontal gene transfers from bacterial sources. Our results provide insights into the adaptations required by metazoans to successfully parasitize immunocompetent plants, and open the way for discovering new antiparasitic strategies.

A Burkholderia cenocepacia MurJ (MviN) homolog is essential for cell wall peptidoglycan synthesis and bacterial viability

The cell wall peptidoglycan (PG) of Burkholderia cenocepacia, an opportunistic pathogen, has not yet been characterized. However, the B. cenocepacia genome contains homologs of genes encoding PG biosynthetic functions in other bacteria. PG biosynthesis involves the formation of the undecaprenyl-pyrophosphate-linked N-acetyl glucosamine-N-acetyl muramic acid-pentapeptide, known as lipid II, which is built on the cytosolic face of the cell membrane. Lipid II is then translocated across the membrane and its glycopeptide moiety becomes incorporated into the growing cell wall mesh; this translocation step is critical to PG synthesis. We have investigated candidate flippase homologs of the MurJ family in B. cenocepacia. Our results show that BCAL2764, herein referred to as murJBc, is indispensable for viability. Viable B. cenocepacia could only be obtained through a conditional mutagenesis strategy by placing murJBc under the control of a rhamnose-inducible promoter. Under rhamnose depletion, the conditional strain stopped growing and individual cells displayed morphological abnormalities consistent with a defect in PG synthesis. Bacterial cells unable to express MurJBc underwent cell lysis, while partial MurJBc depletion sensitized the mutant to the action of β-lactam antibiotics. Depletion of MurJBc caused accumulation of PG precursors consistent with the notion that this protein plays a role in lipid II flipping to the periplasmic compartment. Reciprocal complementation experiments of conditional murJ mutants in B. cenocepacia and Escherichia coli with plasmids expressing MurJ from each strain indicated that MurJBc and MurJEc are functional homologs. Together, our results are consistent with the notion that MurJBc is a PG lipid II flippase in B. cenocepacia.

Discrete cyclic di-GMP-dependent control of bacterial predation versus axenic growth in Bdellovibrio bacteriovorus

Bdellovibrio bacteriovorus is a Delta-proteobacterium that oscillates between free-living growth and predation on Gram-negative bacteria including important pathogens of man, animals and plants. After entering the prey periplasm, killing the prey and replicating inside the prey bdelloplast, several motile B. bacteriovorus progeny cells emerge. The B. bacteriovorus HD100 genome encodes numerous proteins predicted to be involved in signalling via the secondary messenger cyclic di-GMP (c-di-GMP), which is known to affect bacterial lifestyle choices. We investigated the role of c-di-GMP signalling in B. bacteriovorus, focussing on the five GGDEF domain proteins that are predicted to function as diguanylyl cyclases initiating c-di-GMP signalling cascades. Inactivation of individual GGDEF domain genes resulted in remarkably distinct phenotypes. Deletion of dgcB (Bd0742) resulted in a predation impaired, obligately axenic mutant, while deletion of dgcC (Bd1434) resulted in the opposite, obligately predatory mutant. Deletion of dgcA (Bd0367) abolished gliding motility, producing bacteria capable of predatory invasion but unable to leave the exhausted prey. Complementation was achieved with wild type dgc genes, but not with GGAAF versions. Deletion of cdgA (Bd3125) substantially slowed predation; this was restored by wild type complementation. Deletion of dgcD (Bd3766) had no observable phenotype. In vitro assays showed that DgcA, DgcB, and DgcC were diguanylyl cyclases. CdgA lacks enzymatic activity but functions as a c-di-GMP receptor apparently in the DgcB pathway. Activity of DgcD was not detected. Deletion of DgcA strongly decreased the extractable c-di-GMP content of axenic Bdellovibrio cells. We show that c-di-GMP signalling pathways are essential for both the free-living and predatory lifestyles of B. bacteriovorus and that obligately predatory dgcC- can be made lacking a propensity to survive without predation of bacterial pathogens and thus possibly useful in anti-pathogen applications. In contrast to many studies in other bacteria, Bdellovibrio shows specificity and lack of overlap in c-di-GMP signalling pathways.

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers. © 2012 Biek et al.

Genome Expression Profiling-Based Identification and Administration Efficacy of Host-Directed Antimicrobial Drugs against Respiratory Infection by Nontypeable Haemophilus influenzae

Therapies that are safe, effective, and not vulnerable to developing resistance are highly desirable to counteract bacterial infections. Host-directed therapeutics is an antimicrobial approach alternative to conventional antibiotics based on perturbing host pathways subverted by pathogens during their life cycle by using host-directed drugs. In this study, we identified and evaluated the efficacy of a panel of host-directed drugs against respiratory infection by nontypeable Haemophilus influenzae (NTHi). NTHi is an opportunistic pathogen that is an important cause of exacerbation of chronic obstructive pulmonary disease (COPD). We screened for host genes differentially expressed upon infection by the clinical isolate NTHi375 by analyzing cell whole-genome expression profiling and identified a repertoire of host target candidates that were pharmacologically modulated. Based on the proposed relationship between NTHi intracellular location and persistence, we hypothesized that drugs perturbing host pathways used by NTHi to enter epithelial cells could have antimicrobial potential against NTHi infection. Interfering drugs were tested for their effects on bacterial and cellular viability, on NTHi-epithelial cell interplay, and on mouse pulmonary infection. Glucocorticoids and statins lacked in vitro and/or in vivo efficacy. Conversely, the sirtuin-1 activator resveratrol showed a bactericidal effect against NTHi, and the PDE4 inhibitor rolipram showed therapeutic efficacy by lowering NTHi375 counts intracellularly and in the lungs of infected mice. PDE4 inhibition is currently prescribed in COPD, and resveratrol is an attractive geroprotector for COPD treatment. Together, these results expand our knowledge of NTHi-triggered host subversion and frame the antimicrobial potential of rolipram and resveratrol against NTHi respiratory infection.

A genome-wide linkage scan for genes controlling variation in urinary albumin excretion in type ll diabetes

The main hallmark of diabetic nephropathy is elevation in urinary albumin excretion. We performed a genome-wide linkage scan in 63 extended families with multiple members with type II diabetes. Urinary albumin excretion, measured as the albumin-to-creatinine ratio (ACR), was determined in 426 diabetic and 431 nondiabetic relatives who were genotyped for 383 markers. The data were analyzed using variance components linkage analysis. Heritability (h2) of ACR was significant in diabetic (h2=0.23, P=0.0007), and nondiabetic (h2=0.39, P=0.0001) relatives. There was no significant difference in genetic variance of ACR between diabetic and nondiabetic relatives (P=0.16), and the genetic correlation (rG=0.64) for ACR between these two groups was not different from 1 (P=0.12). These results suggested that similar genes contribute to variation in ACR in diabetic and nondiabetic relatives. This hypothesis was supported further by the linkage results.