Bordetella species are distinguished by patterns of substantial gene loss and host adaptation

Genomic evidence and virulence properties decipher the extra-host origin of Bordetella bronchiseptica

Until recently, members of the classical Bordetella species comprised only pathogenic bacteria that were thought to live exclusively in warm-blooded animals. The close phylogenetic relationship of Bordetella with Achromobacter and Alcaligenes, which include primarily environmental bacteria, suggests that the ancestral Bordetellae were probably free-living. Eventually, the Bordetella species evolved to infect and live within warm-blooded animals. The modern history of pathogens related to the genus Bordetella started towards the end of the 19th century when it was discovered in the infected respiratory epithelium of mammals, including humans. The first identified member was Bordetella pertussis, which causes whooping cough, a fatal disease in young children. In due course, B. bronchiseptica was recovered from the trachea and bronchi of dogs with distemper. Later, a second closely related human pathogen, B. parapertussis, was described as causing milder whooping cough. The classical Bordetellae are strictly host-associated pathogens transmitted via the host-to-host aerosol route. Recently, the B. bronchiseptica strain HT200 has been reported from a thermal spring exhibiting unique genomic features that were not previously observed in clinical strains. Therefore, it advocates that members of classical Bordetella species have evolved from environmental sources. This organism can be transmitted via environmental reservoirs as it can survive nutrient-limiting conditions and possesses a motile flagellum. This study aims to review the molecular basis of origin and virulence properties of obligate host-restricted and environmental strains of classical Bordetella.

Comparative genomics of Bordetella pertussis and prediction of new vaccines and drug targets

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis, a Gram-negative bacterium described over a century ago. Despite broad vaccine coverage and treatment options, the disease is remerging as a public health problem especially in infants and older children. Recent data indicate re-emergence of the disease is related to bacterial resistance to immune defences and decreased vaccine effectiveness, which obviously suggests the need of new effective vaccines and drugs. In an attempt to contribute with solutions to this great challenge, bioinformatics tools were used to genetically comprehend the species of these bacteria and predict new vaccines and drug targets. In fact, approaches were used to analysis genomic plasticity, gene synteny and species similarities between the 20 genomes of Bordetella pertussis already available. Furthermore, it was conducted reverse vaccinology and docking analysis to identify proteins with potential to become vaccine and drug targets, respectively. The analyses showed the 20 genomes belongs to a homogeneous group that has preserved most of the genes over time. Besides that, were found genomics islands and good proteins to be candidates for vaccine and drugs. Taken together, these results suggests new possibilities that may be useful to develop new vaccines and drugs that will help the prevention and treatment strategies of pertussis disease caused by these Bordetella strains. Communicated by Ramaswamy H. Sarma.

Population Genomics Reveals Molecular Determinants of Specialization to Tomato in the Polyphagous Fungal Pathogen Botrytis cinerea in France

Many fungal plant pathogens encompass multiple populations specialized on different plant species. Understanding the factors underlying pathogen adaptation to their hosts is a major challenge of evolutionary microbiology, and it should help to prevent the emergence of new specialized pathogens on novel hosts. Previous studies have shown that French populations of the gray mold pathogen Botrytis cinerea parasitizing tomato and grapevine are differentiated from each other, and have higher aggressiveness on their host of origin than on other hosts, indicating some degree of host specialization in this polyphagous pathogen. Here, we aimed at identifying the genomic features underlying the specialization of B. cinerea populations to tomato and grapevine. Based on whole genome sequences of 32 isolates, we confirmed the subdivision of B. cinerea pathogens into two genetic clusters on grapevine and another, single cluster on tomato. Levels of genetic variation in the different clusters were similar, suggesting that the tomato-specific cluster has not recently emerged following a bottleneck. Using genome scans for selective sweeps and divergent selection, tests of positive selection based on polymorphism and divergence at synonymous and nonsynonymous sites, and analyses of presence and absence variation, we identified several candidate genes that represent possible determinants of host specialization in the tomato-associated population. This work deepens our understanding of the genomic changes underlying the specialization of fungal pathogen populations.

A Unique Reverse Adaptation Mechanism Assists Bordetella pertussis in Resistance to Both Scarcity and Toxicity of Manganese

The ability of bacterial pathogens to acquire essential micronutrients is critical for their survival in the host environment. Manganese plays a complex role in the virulence of a variety of pathogens due to its function as an antioxidant and enzymatic cofactor. Therefore, host cells deprive pathogens of manganese to prevent or attenuate infection. Here, we show that evolution of the human-restricted pathogen Bordetella pertussis has selected for an inhibitory duplication within a manganese exporter of the calcium:cation antiporter superfamily. Intriguingly, upon exposure to toxic levels of manganese, the nonfunctional exporter becomes operative in resister cells due to a unique reverse adaptation mechanism. However, compared with wild-type (wt) cells, the resisters carrying a functional copy of the exporter displayed strongly reduced intracellular levels of manganese and impaired growth under oxidative stress. Apparently, inactivation of the manganese exporter and the resulting accumulation of manganese in the cytosol benefited the pathogen by improving its survival under stress conditions. The inhibitory duplication within the exporter gene is highly conserved among B. pertussis strains, absent from all other Bordetella species and from a vast majority of organisms across all kingdoms of life. Therefore, we conclude that inactivation of the exporter gene represents an exceptional example of a flexible genome decay strategy employed by a human pathogen to adapt to its exclusive host.

Genetic Diversity of Clinical Bordetella Pertussis ST2 Strains in comparison with Vaccine Reference Strains of India

Objectives: Pertussis is a highly contagious disease of the respiratory tract caused by Bordetella pertussis, a bacterium that lives in the mouth, nose, and throat. Current study reports the highly accurate complete genomes of two clinical B. pertussis strains from India for the first time.

Methods: Complete genome sequencing was performed for two B. pertussis strains using Ion Torrent PGM and Oxford nanopore sequencing method. Data was assembled de novo and the sequence annotation was performed through PATRIC and NCBI server. Downstream analyses of the isolates were performed using CGE server databases for antimicrobial resistance genes, plasmids, and sequence types. The phylogenetic analysis was performed using Roary.

Results: The analysis revealed insertional elements flanked by IS481, which has been previously regarded as the important component for bacterial evolution. The two B. pertussis clinical strains exhibited diversity through genome degradation when compared to whole-cell vaccine reference strains of India. These isolates harboured multiple genetic virulence traits and toxin subunits, which belonged to sequence type ST2.

Conclusion: The genome information of Indian clinical B. pertussis strains will serve as a baseline data to decipher more information on the genome evolution, virulence factors and their role in pathogenesis for effective vaccine strategies.

Comparative Omics Analysis of Historic and Recent Isolates of Bordetella pertussis and Effects of Genome Rearrangements on Evolution

Despite high vaccination coverage, pertussis is increasing in many industrialized countries, including the Czech Republic. To better understand Bordetella pertussis resurgence, we analyzed historic strains and recent clinical isolates by using a comparative omics approach. Whole-genome sequencing showed that historic and recent isolates of B. pertussis have substantial variation in genome organization and form separate phylogenetic clusters. Subsequent RNA sequence analysis and liquid chromatography with mass tandem spectrometry analyses showed that these variations translated into discretely separated transcriptomic and proteomic profiles. When compared with historic strains, recent isolates showed increased expression of flagellar genes and genes involved in lipopolysaccharide biosynthesis and decreased expression of polysaccharide capsule genes. Compared with reference strain Tohama I, all strains had increased expression and production of the type III secretion system apparatus. We detected the potential link between observed effects and insertion sequence element–induced changes in gene context only for a few genes.

Genomic plasticity and antibody response of Bordetella bronchiseptica strain HT200, a natural variant from a thermal spring

Classical Bordetella species are primarily isolated from animals and humans causing asymptomatic infection to lethal pneumonia. However, isolation of these bacteria from any extra-host environmental niche has not been reported so far. Here, we have characterized the genomic plasticity and antibody response of Bordetella bronchiseptica strain HT200, isolated from a thermal spring. Genomic ANI value and SNPs-based phylogenetic tree suggest a divergent evolution of strain HT200 from a human-adapted lineage of B. bronchiseptica. Growth and survivability assay showed strain HT200 retained viability for more than 5 weeks in the filter-sterilized spring water. In addition, genes or loci encoding the Bordetella virulence factors such as DNT, ACT and LPS O-antigen were absent in strain HT200, while genes encoding other virulence factors were highly divergent. Phenotypically, strain HT200 was non-hemolytic and showed weak hemagglutination activity, but was able to colonize in the respiratory organs of mice. Further, both infection and vaccination with strain HT200 induced protective antibody response in mouse against challenge infection with virulent B. bronchiseptica strain RB50. In addition, genome of strain HT200 (DSM 26023) showed presence of accessory genes and operons encoding predicted metabolic functions pertinent to the ecological conditions of the thermal spring.

Comparative Phosphoproteomics of Classical Bordetellae Elucidates the Potential Role of Serine, Threonine and Tyrosine Phosphorylation in Bordetella Biology and Virulence

The Bordetella genus is divided into two groups: classical and non-classical. Bordetella pertussis, Bordetella bronchiseptica and Bordetella parapertussis are known as classical bordetellae, a group of important human pathogens causing whooping cough or whooping cough-like disease and hypothesized to have evolved from environmental non-classical bordetellae. Bordetella infections have increased globally driving the need to better understand these pathogens for the development of new treatments and vaccines. One unexplored component in Bordetella is the role of serine, threonine and tyrosine phosphorylation. Therefore, this study characterized the phosphoproteome of classical bordetellae and examined its potential role in Bordetella biology and virulence. Applying strict identification of localization criteria, this study identified 70 unique phosphorylated proteins in the classical bordetellae group with a high degree of conservation. Phosphorylation was a key regulator of Bordetella metabolism with proteins involved in gluconeogenesis, TCA cycle, amino acid and nucleotide synthesis significantly enriched. Three key virulence pathways were also phosphorylated including type III secretion system, alcaligin synthesis and the BvgAS master transcriptional regulatory system for virulence genes in Bordetella. Seven new phosphosites were identified in BvgA with 6 located in the DNA binding domain. Of the 7, 4 were not present in non-classical bordetellae. This suggests that serine/threonine phosphorylation may play an important role in stabilizing/destabilizing BvgA binding to DNA for fine-tuning of virulence gene expression and that BvgA phosphorylation may be an important factor separating classical from non-classical bordetellae. This study provides the first insight into the phosphoproteome of classical Bordetella species and the role that Ser/Thr/Tyr phosphorylation may play in Bordetella biology and virulence.

Bordetella Type III Secretion Injectosome and Effector Proteins

Pertussis, also known as whooping cough, is a resurging acute respiratory disease of humans primarily caused by the Gram-negative coccobacilli Bordetella pertussis, and less commonly by the human-adapted lineage of B. parapertussis _HU. The ovine-adapted lineage of B. parapertussis _OV infects only sheep, while B. bronchiseptica causes chronic and often asymptomatic respiratory infections in a broad range of mammals but rarely in humans. A largely overlapping set of virulence factors inflicts the pathogenicity of these bordetellae. Their genomes also harbor a pathogenicity island, named bsc locus, that encodes components of the type III secretion injectosome, and adjacent btr locus with the type III regulatory proteins. The Bsc injectosome of bordetellae translocates the cytotoxic BteA effector protein, also referred to as BopC, into the cells of the mammalian hosts. While the role of type III secretion activity in the persistent colonization of the lower respiratory tract by B. bronchiseptica is well recognized, the functionality of the type III secretion injectosome in B. pertussis was overlooked for many years due to the adaptation of laboratory-passaged B. pertussis strains. This review highlights the current knowledge of the type III secretion system in the so-called classical Bordetella species, comprising B. pertussis, B. parapertussis, and B. bronchiseptica, and discusses its functional divergence. Comparison with other well-studied bacterial injectosomes, regulation of the type III secretion on the transcriptional and post-transcriptional level, and activities of BteA effector protein and BopN protein, homologous to the type III secretion gatekeepers, are addressed.

Comparative genomics of whole-cell pertussis vaccine strains from India

BACKGROUND

Despite high vaccination coverage using acellular (ACV) and whole-cell pertussis (WCV) vaccines, the resurgence of pertussis is observed globally. Genetic divergence in circulating strains of Bordetella pertussis has been reported as one of the contributing factors for the resurgence of the disease. Our current knowledge of B. pertussis genetic evolution in circulating strains is mostly based on studies conducted in countries using ACVs targeting only a few antigens used in the production of ACVs. To better understand the adaptation to vaccine-induced selection pressure, it will be essential to study B. pertussis populations in developing countries which are using WCVs. India is a significant user and global supplier of WCVs. We report here comparative genome analyses of vaccine and clinical isolates reported from India. Whole-genome sequences obtained from vaccine strains: WCV (J445, J446, J447 and J448), ACV (BP165) were compared with Tohama-I reference strain and recently reported clinical isolates from India (BPD1, BPD2). Core genome-based phylogenetic analysis was also performed using 166 isolates reported from countries using ACV.

RESULTS

Whole-genome analysis of vaccine and clinical isolates reported from India revealed high genetic similarity and conserved genome among strains. Phylogenetic analysis showed that clinical and vaccine strains share genetic closeness with reference strain Tohama-I. The allelic profile of vaccine strains (J445:ptxP1/ptxA2/prn1/fim2-1/fim3-1; J446: ptxP2/ptxA4/prn7/fim2-2/fim3-1; J447 and J448: ptxP1/ptxA1/ prn1/fim2-1/fim3-1), which matched entirely with clinical isolates (BPD1:ptxP1/ptxA1/prn1/fim2-1 and BPD2: ptxP1/ptxA1/prn1/fim2-1) reported from India. Multi-locus sequence typing (MLST) demonstrated the presence of dominant sequence types ST2 and primitive ST1 in vaccine strains which will allow better coverage against circulating strains of B. pertussis.

CONCLUSIONS

The study provides a detailed characterization of vaccine and clinical strains reported from India, which will further facilitate epidemiological studies on genetic shifts in countries which are using WCVs in their immunization programs.

Genomic epidemiology of Iranian Bordetella pertussis: 50 years after the implementation of whole cell vaccine

Pertussis caused by Bordetella pertussis, remains a public health problem worldwide, despite high vaccine coverage in infants and children in many countries. Iran has been using whole cell vaccine for the last 50 years with more than 95% vaccination rate since 1988 and has experienced pertussis resurgence in recent years. Here, we sequenced 55 B. pertussis isolates mostly collected from three provinces with the highest number of pertussis cases in Iran, including Tehran, Mazandaran, and Eastern-Azarbayjan from the period of 2008-2016. Most isolates carried ptxP3/prn2 alleles (42/55, 76%), the same genotype as isolates circulating in acellular vaccine-administrating countries. The second most frequent genotype was ptxP3/prn9 (8/55, 14%). Only three isolates (5%) were ptxP1. Phylogenetic analysis showed that Iranian ptxP3 isolates can be divided into eight clades (Clades 1-8) with no temporal association. Most of the isolates from Tehran grouped together as one distinctive clade (Clade 8) with six unique single nucleotide polymorphisms (SNPs). In addition, the prn9 isolates were grouped together as Clade 5 with 12 clade-supporting SNPs. No pertactin deficient isolates were found among the 55 Iranian isolates. Our findings suggest that there is an ongoing adaptation and evolution of B. pertussis regardless of the types of vaccine used.

The BvgASR virulence regulon of Bordetella pertussis

The BvgAS two-component system of Bordetella pertussis directly activates the expression of a large number of virulence genes in an environmentally responsive manner. The Bvg⁺ mode also promotes the expression of the phosphodiesterase BvgR, which turns off the expression of another set of genes, the vrgs, by reducing levels of c-di-GMP. Increased levels of c-di-GMP in the Bvg^- mode are required, together with the phosphorylated response regulator protein RisA∼P, to activate vrg expression. Phosphorylation of RisA requires RisK, a non-co-operonic sensor kinase, but not its co-operonic sensor kinase RisS which is truncated in B. pertussis but intact in the ancestral B. bronchiseptica. The loss of RisS during evolution of B. pertussis led to the ability to express the vrgs, potentially enhancing aerosol transmission of B. pertussis.

Comparative genomics of Czech vaccine strains of Bordetella pertussis

Bordetella pertussis is a strictly human pathogen causing the respiratory infectious disease called whooping cough or pertussis. B. pertussis adaptation to acellular pertussis vaccine pressure has been repeatedly highlighted, but recent data indicate that adaptation of circulating strains started already in the era of the whole cell pertussis vaccine (wP) use. We sequenced the genomes of five B. pertussis wP vaccine strains isolated in the former Czechoslovakia in the pre-wP (1954-1957) and early wP (1958-1965) eras, when only limited population travel into and out of the country was possible. Four isolates exhibit a similar genome organization and form a distinct phylogenetic cluster with a geographic signature. The fifth strain is rather distinct, both in genome organization and SNP-based phylogeny. Surprisingly, despite isolation of this strain before 1966, its closest sequenced relative appears to be a recent isolate from the US. On the genome content level, the five vaccine strains contained both new and already described regions of difference. One of the new regions contains duplicated genes potentially associated with transport across the membrane. The prevalence of this region in recent isolates indicates that its spread might be associated with selective advantage leading to increased strain fitness.

A Novel Bvg-Repressed Promoter Causes vrg-Like Transcription of fim3 but Does Not Result in the Production of Serotype 3 Fimbriae in Bvg- Mode Bordetella pertussis

In Bordetella pertussis, two serologically distinct fimbriae, FIM2 and FIM3, undergo on/off phase variation independently of each other via variation in the lengths of C stretches in the promoters for their major subunit genes, fim2 and fim3 These two promoters are also part of the BvgAS virulence regulon and therefore, if in an on configuration, are activated by phosporylated BvgA (BvgA~P) under normal growth conditions (Bvg⁺ mode) but not in the Bvg^- mode, inducible by growth in medium containing MgSO₄ or other compounds, termed modulators. In the B. pertussis Tohama I strain (FIM2⁺ FIM3^-), the fim3 promoter is in the off state. However, a high level of transcription of the fim3 gene is observed in the Bvg^- mode. In this study, we provide an explanation for this anomalous behavior by defining a Bvg-repressed promoter (BRP), located approximately 400 bp upstream of the Pfim3 transcriptional start. Although transcription of the fim3 gene in the Bvg^- mode resulted in Fim3 translation, as measured by LacZ translational fusions, no accumulation of Fim3 protein was detectable. We propose that Fim3 protein resulting from translation of mRNA driven by BRP in the Bvg^- mode is unstable due to a lack of the fimbrial assembly apparatus encoded by the fimBC genes, located within the fha operon, and therefore is not expressed in the Bvg^- mode.IMPORTANCE In Bordetella pertussis, the promoter Pfim3-15C for the major fimbrial subunit gene fim3 is activated by the two-component system BvgAS in the Bvg⁺ mode but not in the Bvg^- mode. However, many transcriptional profiling studies have shown that fim3 is transcribed in the Bvg^- mode even when Pfim3 is in a nonpermissive state (Pfim3-13C), suggesting the presence of a reciprocally regulated element upstream of Pfim3 Here, we provide evidence that BRP is the cause of this anomalous behavior of fim3 Although BRP effects vrg-like transcription of fim3 in the Bvg^- mode, it does not lead to stable production of FIM3 fimbriae, because expression of the chaperone and usher proteins FimB and FimC occurs only in the Bvg⁺ mode.

The Transcriptional Regulator BpsR Controls the Growth of Bordetella bronchiseptica by Repressing Genes Involved in Nicotinic Acid Degradation

Many of the pathogenic species of the genus Bordetella have an absolute requirement for nicotinic acid (NA) for laboratory growth. These Gram-negative bacteria also harbor a gene cluster homologous to the nic cluster of Pseudomonas putida which is involved in the aerobic degradation of NA and its transcriptional control. We report here that BpsR, a negative regulator of biofilm formation and Bps polysaccharide production, controls the growth of Bordetella bronchiseptica by repressing the expression of nic genes. The severe growth defect of the ΔbpsR strain in Stainer-Scholte medium was restored by supplementation with NA, which also functioned as an inducer of nic genes at low micromolar concentrations that are usually present in animals and humans. Purified BpsR protein bound to the nic promoter region, and its DNA binding activity was inhibited by 6-hydroxynicotinic acid (6-HNA), the first metabolite of the NA degradative pathway. Reporter assays with the isogenic mutant derivative of the wild-type (WT) strain harboring deletion in nicA, which encodes a putative nicotinic acid hydroxylase responsible for conversion of NA to 6-HNA, showed that 6-HNA is the actual inducer of the nic genes in the bacterial cell. Gene expression profiling further showed that BpsR dually activated and repressed the expression of genes associated with pathogenesis, transcriptional regulation, metabolism, and other cellular processes. We discuss the implications of these findings with respect to the selection of pyridines such as NA and quinolinic acid for optimum bacterial growth depending on the ecological niche.IMPORTANCE BpsR, the previously described regulator of biofilm formation and Bps polysaccharide production, controls Bordetella bronchiseptica growth by regulating the expression of genes involved in the degradation of nicotinic acid (NA). 6-Hydroxynicotinic acid (6-HNA), the first metabolite of the NA degradation pathway prevented BpsR from binding to DNA and was the actual in vivo inducer. We hypothesize that BpsR enables Bordetella bacteria to efficiently and selectively utilize NA for their survival depending on the environment in which they reside. The results reported herein lay the foundation for future investigations of how BpsR and the alteration of its activity by NA orchestrate the control of Bordetella growth, metabolism, biofilm formation, and pathogenesis.

Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis

Bordetella pertussis is the causative agent of whooping cough, a respiratory disease still considered as a major public health threat and for which recent re-emergence has been observed. Constant reshuffling of Bordetella pertussis genome organization was observed during evolution. These rearrangements are essentially mediated by Insertion Sequences (IS), a mobile genetic elements present in more than 230 copies in the genome, which are supposed to be one of the driving forces enabling the pathogen to escape from vaccine-induced immunity. Here we use high-throughput sequencing approaches (RNA-seq and differential RNA-seq), to decipher Bordetella pertussis transcriptome characteristics and to evaluate the impact of IS elements on transcriptome architecture. Transcriptional organization was determined by identification of transcription start sites and revealed also a large variety of non-coding RNAs including sRNAs, leaderless mRNAs or long 3' and 5'UTR including seven riboswitches. Unusual topological organizations, such as overlapping 5'- or 3'-extremities between oppositely orientated mRNA were also unveiled. The pivotal role of IS elements in the transcriptome architecture and their effect on the transcription of neighboring genes was examined. This effect is mediated by the introduction of IS harbored promoters or by emergence of hybrid promoters. This study revealed that in addition to their impact on genome rearrangements, most of the IS also impact on the expression of their flanking genes. Furthermore, the transcripts produced by IS are strain-specific due to the strain to strain variation in IS copy number and genomic context.

Ectopic Expression of O Antigen in Bordetella pertussis by a Novel Genomic Integration System

Some bacterial phenotypes emerge through the cooperative functions of a number of genes residing within a large genetic locus. To transfer the phenotype of one bacterium to another, a means to introduce the large genetic locus into the recipient bacterium is needed. Therefore, we developed a novel system by combining the advantages of a bacterial artificial chromosome vector and phage-derived gene integration machinery. In this study, we succeeded for the first time in introducing a gene locus involved in O antigen biosynthesis of Bordetella bronchiseptica into the chromosome of B. pertussis, which intrinsically lacks O antigen, and using this system we analyzed phenotypic alterations in the resultant mutant strain of B. pertussis. The present results demonstrate that this system successfully accomplished the above-described purpose. We consider this system to be applicable to a number of bacteria other than Bordetella.

We describe a novel genome integration system that enables the introduction of DNA fragments as large as 50 kbp into the chromosomes of recipient bacteria. This system, named BPI, comprises a bacterial artificial chromosome vector and phage-derived gene integration machinery. We introduced the wbm locus of Bordetella bronchiseptica, which is required for O antigen biosynthesis, into the chromosome of B. pertussis, which intrinsically lacks O antigen, using the BPI system. After the introduction of the wbm locus, B. pertussis presented an additional substance in the lipooligosaccharide fraction that was specifically recognized by the anti-B. bronchiseptica antibody but not the anti-B. pertussis antibody, indicating that B. pertussis expressed O antigen corresponding to that of B. bronchiseptica. O antigen-expressing B. pertussis was less sensitive to the bactericidal effects of serum and polymyxin B than the isogenic parental strain. In addition, an in vivo competitive infection assay showed that O antigen-expressing B. pertussis dominantly colonized the mouse respiratory tract over the parental strain. These results indicate that the BPI system provides a means to alter the phenotypes of bacteria by introducing large exogenous DNA fragments.

IMPORTANCE Some bacterial phenotypes emerge through the cooperative functions of a number of genes residing within a large genetic locus. To transfer the phenotype of one bacterium to another, a means to introduce the large genetic locus into the recipient bacterium is needed. Therefore, we developed a novel system by combining the advantages of a bacterial artificial chromosome vector and phage-derived gene integration machinery. In this study, we succeeded for the first time in introducing a gene locus involved in O antigen biosynthesis of Bordetella bronchiseptica into the chromosome of B. pertussis, which intrinsically lacks O antigen, and using this system we analyzed phenotypic alterations in the resultant mutant strain of B. pertussis. The present results demonstrate that this system successfully accomplished the above-described purpose. We consider this system to be applicable to a number of bacteria other than Bordetella.

Characterization of a Bvg-regulated fatty acid methyl-transferase in Bordetella pertussis

The whooping cough agent Bordetella pertussis controls the expression of its large virulence regulon in a coordinated manner through the two-component signal transduction system BvgAS. In addition to the genes coding for bona fide virulence factors, the Bvg regulon comprises genes of unknown function. In this work, we characterized a new Bvg-activated gene called BP2936. Homologs of BP2936 are found in other pathogenic Bordetellae and in several other species, including plant pathogens and environmental bacteria. We showed that the gene product of BP2936 is a membrane-associated methyl-transferase of free fatty acids. We thus propose to name it FmtB, for fatty acid methyl-transferase of Bordetella. The role of this protein was tested in cellular and animal models of infection, but the loss of BP2936 did not appear to affect host-pathogen interactions in those assays. The high level of conservation of BP2936 among B. pertussis isolates nevertheless argues that it probably plays a role in the life cycle of this pathogen.

The extent of the temperature-induced membrane remodeling in two closely related Bordetella species reflects their adaptation to diverse environmental niches

Changes in environmental temperature represent one of the major stresses faced by microorganisms as they affect the function of the cytoplasmic membrane. In this study, we have analyzed the thermal adaptation in two closely related respiratory pathogens Bordetella pertussis and Bordetella bronchiseptica Although B. pertussis represents a pathogen strictly adapted to the human body temperature, B. bronchiseptica causes infection in a broad range of animals and survives also outside of the host. We applied GC-MS to determine the fatty acids of both Bordetella species grown at different temperatures and analyzed the membrane fluidity by fluorescence anisotropy measurement. In parallel, we also monitored the effect of growth temperature changes on the expression and production of several virulence factors. In response to low temperatures, B. pertussis adapted its fatty acid composition and membrane fluidity to a considerably lesser extent when compared with B. bronchiseptica Remarkably, B. pertussis maintained the production of virulence factors at 24 °C, whereas B. bronchiseptica cells resumed the production only upon temperature upshift to 37 °C. This growth temperature-associated differential modulation of virulence factor production was linked to the phosphorylation state of transcriptional regulator BvgA. The observed differences in low-temperature adaptation between B. pertussis and B. bronchiseptica may result from selective adaptation of B. pertussis to the human host. We propose that the reduced plasticity of the B. pertussis membranes ensures sustained production of virulence factors at suboptimal temperatures and may play an important role in the transmission of the disease.

Etiologic Framework for the Study of Neurodegenerative Disorders as Well as Vascular and Metabolic Comorbidities on the Grounds of Shared Epidemiologic and Biologic Features

Background: During the last two decades, protein aggregation at all organismal levels, from viruses to humans, has emerged from a neglected area of protein science to become a central issue in biology and biomedicine. This article constitutes a risk-based review aimed at supporting an etiologic scenario of selected, sporadic, protein-associated, i.e., conformational, neurodegenerative disorders (NDDs), and their vascular- and metabolic-associated ailments.

Methods: A rationale is adopted, to incorporate selected clinical data and results from animal-model research, complementing epidemiologic evidences reported in two prior articles.

Findings: Theory is formulated assuming an underlying conformational transmission mechanism, mediated either by horizontal transfer of mammalian genes coding for specific aggregation-prone proteins, or by xeno-templating between bacterial and host proteins. We build a few population-based and experimentally-testable hypotheses focusing on: (1) non-disposable surgical instruments for sporadic Creutzfeldt-Jakob disease (sCJD) and other rapid progressive neurodegenerative dementia (sRPNDd), multiple system atrophy (MSA), and motor neuron disease (MND); and (2) specific bacterial infections such as B. pertussis and E. coli for all forms, but particularly for late-life sporadic conformational, NDDs, type 2 diabetes mellitus (T2DM), and atherosclerosis where natural protein fibrils present in such organisms as a result of adaptation to the human host induce prion-like mechanisms.

Conclusion: Implications for cohort alignment and experimental animal research are discussed and research lines proposed.

Genomic perspectives on the evolution and spread of bacterial pathogens

Since the first complete sequencing of a free-living organism, Haemophilus influenzae, genomics has been used to probe both the biology of bacterial pathogens and their evolution. Single-genome approaches provided information on the repertoire of virulence determinants and host-interaction factors, and, along with comparative analyses, allowed the proposal of hypotheses to explain the evolution of many of these traits. These analyses suggested many bacterial pathogens to be of relatively recent origin and identified genome degradation as a key aspect of host adaptation. The advent of very-high-throughput sequencing has allowed for detailed phylogenetic analysis of many important pathogens, revealing patterns of global and local spread, and recent evolution in response to pressure from therapeutics and the human immune system. Such analyses have shown that bacteria can evolve and transmit very rapidly, with emerging clones showing adaptation and global spread over years or decades. The resolution achieved with whole-genome sequencing has shown considerable benefits in clinical microbiology, enabling accurate outbreak tracking within hospitals and across continents. Continued large-scale sequencing promises many further insights into genetic determinants of drug resistance, virulence and transmission in bacterial pathogens.

The Bordetella Secreted Regulator BspR Is Translocated into the Nucleus of Host Cells via Its N-Terminal Moiety: Evaluation of Bacterial Effector Translocation by the Escherichia coli Type III Secretion System

Bordetella bronchiseptica is genetically related to B. pertussis and B. parapertussis, which cause respiratory tract infections in humans. These pathogens possess a large number of virulence factors, including the type III secretion system (T3SS), which is required for the delivery of effectors into the host cells. In a previous study, we identified a transcriptional regulator, BspR, that is involved in the regulation of the T3SS-related genes in response to iron-starved conditions. A unique feature of BspR is that this regulator is secreted into the extracellular milieu via the T3SS. To further characterize the role of BspR in extracellular localization, we constructed various truncated derivatives of BspR and investigated their translocation into the host cells using conventional translocation assays. In this study, the effector translocation was evaluated by the T3SS of enteropathogenic E. coli (EPEC), since the exogenous expression of BspR triggers severe repression of the Bordetella T3SS expression. The results of the translocation assays using the EPEC T3SS showed that the N-terminal 150 amino acid (aa) residues of BspR are sufficient for translocation into the host cells in a T3SS-dependent manner. In addition, exogenous expression of BspR in HeLa cells demonstrated that the N-terminal 100 aa residues are involved in the nuclear localization. In contrast, the N-terminal 54 aa residues are sufficient for the extracellular secretion into the bacterial culture supernatant via the EPEC T3SS. Thus, BspR is not only a transcriptional regulator in bacteria cytosol, but also functions as an effector that translocates into the nuclei of infected host cells.

Genomic characterization of Sinorhizobium meliloti AK21, a wild isolate from the Aral Sea Region

The symbiotic, nitrogen-fixing bacterium Sinorhizobium meliloti has been widely studied due to its ability to improve crop yields through direct interactions with leguminous plants. S. meliloti AK21 is a wild type strain that forms nodules on Medicago plants in saline and drought conditions in the Aral Sea Region. The aim of this work was to establish the genetic similarities and differences between S. meliloti AK21 and the reference strain S. meliloti 1021. Comparative genome hybridization with the model reference strain S. meliloti 1021 yielded 365 variable genes, grouped into 11 regions in the three main replicons in S. meliloti AK21. The most extensive regions of variability were found in the symbiotic plasmid pSymA, which also contained the largest number of orthologous and polymorphic sequences identified by suppression subtractive hybridization. This procedure identified a large number of divergent sequences and others without homology in the databases, the further investigation of which could provide new insight into the alternative metabolic pathways present in S. meliloti AK21. We identified a plasmid replication module from the repABC replicon family, together with plasmid mobilization-related genes (traG and a VirB9-like protein), which suggest that this indigenous isolate harbors an accessory plasmid. Furthermore, the transcriptomic profiles reflected differences in gene content and regulation between S. meliloti AK21 and S. meliloti 1021 (ExpR and PhoB regulons), but provided evidence for an as yet unknown, alternative mechanism involving activation of the cbb3 terminal oxidase. Finally, phenotypic microarrays characterization revealed a greater versatility of substrate use and chemical degradation than for S. meliloti 1021.

Polymorphisms influencing expression of dermonecrotic toxin in Bordetella bronchiseptica

Bordetella bronchiseptica is a pathogenic bacterium causing respiratory infections in a broad range of mammals. Recently, we determined the whole genome sequence of B. bronchiseptica S798 strain isolated from a pig infected with atrophic rhinitis and found four single-nucleotide polymorphisms (SNPs) at positions -129, -72, +22, and +38 in the region upstream of dnt encoding dermonecrotic toxin (DNT), when compared with a rabbit isolate, RB50. DNT is known to be involved in turbinate atrophy observed in atrophic rhinitis. Immunoblotting, quantitative real-time PCR, and β-galactosidase reporter assay revealed that these SNPs resulted in the increased promoter activity of dnt and conferred the increased ability to produce DNT on the bacteria. Similar or identical SNPs were also found in other pig isolates kept in our laboratory, all of which produce a larger amount of DNT than RB50. Our analysis revealed that substitution of at least two of the four bases, at positions -72 and +22, influenced the promoter activity for dnt. These results imply that these SNPs are involved in the pathogenicity of bordetellae specific to pig diseases.

Investigating genome reduction of Bordetella pertussis using a multiplex PCR-based reverse line blot assay (mPCR/RLB)

Background

The genetic composition of the bacterium causing whooping cough, Bordetella pertussis, has been investigated using microarray studies in order to examine potential genetic contributors to the disease re-emergence in the past decade. Regions of difference (RDs) have been previously identified as clusters of genes flanked by insertion sequences which are variably present in different sets of isolates, and have also been shown to be potential markers of B. pertussis evolution.

This study used microarray data to identify and select a panel of RDs; primers and probes for these RDs were then designed to test for the presence or absence of these regions in a novel and less expensive multiplex PCR-based reverse line blot (mPCR/RLB) assay. By comparing the presence or absence of RDs, we aimed to determine the genomic variability of a diverse collection of B. pertussis strains and how they have changed over time.

Results

A B. pertussis specific mPCR/RLB using 43 genes representing 30 RDs, was developed and used to characterise a set of 42 B. pertussis isolates. When mapped against the previously identified evolutionary relationships of the strains, the losses of two RDs - BP0910A - BP00930 and BP1948-BP1962 - were found to be associated with significant events in B. pertussis history: the loss of BP0910A - BP00930 coincided with introduction of whole cell vaccines in the 1950s while that of BP1948-BP1962 occurred after the introduction of acellular vaccines. The loss of BP1948-BP1962 also coincided with expansion of the most recent B. pertussis strains.

Conclusions

The mPCR/RLB assay offers an inexpensive and fast method of determining the gene content of B. pertussis strains and also confirms that gene losses are an ongoing feature of B. pertussis evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-727) contains supplementary material, which is available to authorized users.

Complete Genome Sequence of Bordetella bronchiseptica S798, an Isolate from a Pig with Atrophic Rhinitis

Bordetella bronchiseptica colonizes the respiratory tracts of a wide variety of mammals and causes a range of diseases, from lethal pneumonia to asymptomatic chronic infection. We report the complete genome sequence of Bordetella bronchiseptica strain S798, isolated from a pig with atrophic rhinitis in Japan.

Genome implosion elicits host-confinement in Alcaligenaceae: evidence from the comparative genomics of Tetrathiobacter kashmirensis, a pathogen in the making

This study elucidates the genomic basis of the evolution of pathogens alongside free-living organisms within the family Alcaligenaceae of Betaproteobacteria. Towards that end, the complete genome sequence of the sulfur-chemolithoautotroph Tetrathiobacter kashmirensis WT001^T was determined and compared with the soil isolate Achromobacter xylosoxidans A8 and the two pathogens Bordetella bronchiseptica RB50 and Taylorella equigenitalis MCE9. All analyses comprehensively indicated that the RB50 and MCE9 genomes were almost the subsets of A8 and WT001^T, respectively. In the immediate evolutionary past Achromobacter and Bordetella shared a common ancestor, which was distinct from the other contemporary stock that gave rise to Tetrathiobacter and Taylorella. The Achromobacter-Bordetella precursor, after diverging from the family ancestor, evolved through extensive genome inflation, subsequent to which the two genera separated via differential gene losses and acquisitions. Tetrathiobacter, meanwhile, retained the core characteristics of the family ancestor, and Taylorella underwent massive genome degeneration to reach an evolutionary dead-end. Interestingly, the WT001^T genome, despite its conserved architecture, had only 85% coding density, besides which 578 out of its 4452 protein-coding sequences were found to be pseudogenized. Translational impairment of several DNA repair-recombination genes in the first place seemed to have ushered the rampant and indiscriminate frame-shift mutations across the WT001^T genome. Presumably, this strain has just come out of a recent evolutionary bottleneck, representing a unique transition state where genome self-degeneration has started comprehensively but selective host-confinement has not yet set in. In the light of this evolutionary link, host-adaptation of Taylorella clearly appears to be the aftereffect of genome implosion in another member of the same bottleneck. Remarkably again, potent virulence factors were found widespread in Alcaligenaceae, corroborating which hemolytic and mammalian cell-adhering abilities were discovered in WT001^T. So, while WT001^T relatives/derivatives in nature could be going the Taylorella way, the lineage as such was well-prepared for imminent host-confinement.

Predicting the proteins of Angomonas deanei, Strigomonas culicis and their respective endosymbionts reveals new aspects of the trypanosomatidae family

Endosymbiont-bearing trypanosomatids have been considered excellent models for the study of cell evolution because the host protozoan co-evolves with an intracellular bacterium in a mutualistic relationship. Such protozoa inhabit a single invertebrate host during their entire life cycle and exhibit special characteristics that group them in a particular phylogenetic cluster of the Trypanosomatidae family, thus classified as monoxenics. In an effort to better understand such symbiotic association, we used DNA pyrosequencing and a reference-guided assembly to generate reads that predicted 16,960 and 12,162 open reading frames (ORFs) in two symbiont-bearing trypanosomatids, Angomonas deanei (previously named as Crithidia deanei) and Strigomonas culicis (first known as Blastocrithidia culicis), respectively. Identification of each ORF was based primarily on TriTrypDB using tblastn, and each ORF was confirmed by employing getorf from EMBOSS and Newbler 2.6 when necessary. The monoxenic organisms revealed conserved housekeeping functions when compared to other trypanosomatids, especially compared with Leishmania major. However, major differences were found in ORFs corresponding to the cytoskeleton, the kinetoplast, and the paraflagellar structure. The monoxenic organisms also contain a large number of genes for cytosolic calpain-like and surface gp63 metalloproteases and a reduced number of compartmentalized cysteine proteases in comparison to other TriTryp organisms, reflecting adaptations to the presence of the symbiont. The assembled bacterial endosymbiont sequences exhibit a high A+T content with a total of 787 and 769 ORFs for the Angomonas deanei and Strigomonas culicis endosymbionts, respectively, and indicate that these organisms hold a common ancestor related to the Alcaligenaceae family. Importantly, both symbionts contain enzymes that complement essential host cell biosynthetic pathways, such as those for amino acid, lipid and purine/pyrimidine metabolism. These findings increase our understanding of the intricate symbiotic relationship between the bacterium and the trypanosomatid host and provide clues to better understand eukaryotic cell evolution.

Comparative genomics of the classical Bordetella subspecies: the evolution and exchange of virulence-associated diversity amongst closely related pathogens

Background

The classical Bordetella subspecies are phylogenetically closely related, yet differ in some of the most interesting and important characteristics of pathogens, such as host range, virulence and persistence. The compelling picture from previous comparisons of the three sequenced genomes was of genome degradation, with substantial loss of genome content (up to 24%) associated with adaptation to humans.

Results

For a more comprehensive picture of lineage evolution, we employed comparative genomic and phylogenomic analyses using seven additional diverse, newly sequenced Bordetella isolates. Genome-wide single nucleotide polymorphism (SNP) analysis supports a reevaluation of the phylogenetic relationships between the classical Bordetella subspecies, and suggests a closer link between ovine and human B. parapertussis lineages than has been previously proposed. Comparative analyses of genome content revealed that only 50% of the pan-genome is conserved in all strains, reflecting substantial diversity of genome content in these closely related pathogens that may relate to their different host ranges, virulence and persistence characteristics. Strikingly, these analyses suggest possible horizontal gene transfer (HGT) events in multiple loci encoding virulence factors, including O-antigen and pertussis toxin (Ptx). Segments of the pertussis toxin locus (ptx) and its secretion system locus (ptl) appear to have been acquired by the classical Bordetella subspecies and are divergent in different lineages, suggesting functional divergence in the classical Bordetellae.

Conclusions

Together, these observations, especially in key virulence factors, reveal that multiple mechanisms, such as point mutations, gain or loss of genes, as well as HGTs, contribute to the substantial phenotypic diversity of these versatile subspecies in various hosts.

Small mutations in Bordetella pertussis are associated with selective sweeps

Bordetella pertussis is the causative agent of pertussis, a highly contagious disease of the human respiratory tract. Despite high vaccination coverage, pertussis has resurged and has become one of the most prevalent vaccine-preventable diseases in developed countries. We have proposed that both waning immunity and pathogen adaptation have contributed to the persistence and resurgence of pertussis. Allelic variation has been found in virulence-associated genes coding for the pertussis toxin A subunit (ptxA), pertactin (prn), serotype 2 fimbriae (fim2), serotype 3 fimbriae (fim3) and the promoter for pertussis toxin (ptxP). In this study, we investigated how more than 60 years of vaccination has affected the Dutch B. pertussis population by combining data from phylogeny, genomics and temporal trends in strain frequencies. Our main focus was on the ptxA, prn, fim3 and ptxP genes. However, we also compared the genomes of 11 Dutch strains belonging to successful lineages. Our results showed that, between 1949 and 2010, the Dutch B. pertussis population has undergone as least four selective sweeps that were associated with small mutations in ptxA, prn, fim3 and ptxP. Phylogenetic analysis revealed a stepwise adaptation in which mutations accumulated clonally. Genomic analysis revealed a number of additional mutations which may have a contributed to the selective sweeps. Five large deletions were identified which were fixed in the pathogen population. However, only one was linked to a selective sweep. No evidence was found for a role of gene acquisition in pathogen adaptation. Our results suggest that the B. pertussis gene repertoire is already well adapted to its current niche and required only fine tuning to persist in the face of vaccination. Further, this work shows that small mutations, even single SNPs, can drive large changes in the populations of bacterial pathogens within a time span of six to 19 years.

Phenotypic and genomic analysis of hypervirulent human-associated Bordetella bronchiseptica

Background

B. bronchiseptica infections are usually associated with wild or domesticated animals, but infrequently with humans. A recent phylogenetic analysis distinguished two distinct B. bronchiseptica subpopulations, designated complexes I and IV. Complex IV isolates appear to have a bias for infecting humans; however, little is known regarding their epidemiology, virulence properties, or comparative genomics.

Results

Here we report a characterization of the virulence of human-associated complex IV B. bronchiseptica strains. In in vitro cytotoxicity assays, complex IV strains showed increased cytotoxicity in comparison to a panel of complex I strains. Some complex IV isolates were remarkably cytotoxic, resulting in LDH release levels in A549 cells that were 10- to 20-fold greater than complex I strains. In vivo, a subset of complex IV strains was found to be hypervirulent, with an increased ability to cause lethal pulmonary infections in mice. Hypercytotoxicity in vitro and hypervirulence in vivo were both dependent on the activity of the bsc T3SS and the BteA effector. To clarify differences between lineages, representative complex IV isolates were sequenced and their genomes were compared to complex I isolates. Although our analysis showed there were no genomic sequences that can be considered unique to complex IV strains, there were several loci that were predominantly found in complex IV isolates.

Conclusion

Our observations reveal a T3SS-dependent hypervirulence phenotype in human-associated complex IV isolates, highlighting the need for further studies on the epidemiology and evolutionary dynamics of this B. bronchiseptica lineage.

The influence of rickettsiologists on post-modern microbiology

Many of the definitions in microbiology are currently false. We have reviewed the great denominations of microbiology and attempted to free microorganisms from the theories of the twentieth century. The presence of compartmentation and a nucleoid in Planctomycetes clearly calls into question the accuracy of the definitions of eukaryotes and prokaryotes. Archaea are viewed as prokaryotes resembling bacteria. However, the name archaea, suggesting an archaic origin of lifestyle, is inconsistent with the lifestyle of this family. Viruses are defined as small, filterable infectious agents, but giant viruses challenge the size criteria used for the definition of a virus. Pathogenicity does not require the acquisition of virulence factors (except for toxins), and in many cases, gene loss is significantly inked to the emergence of virulence. Species classification based on 16S rRNA is useless for taxonomic purposes of human pathogens, as a 2% divergence would classify all Rickettsiae within the same species and would not identify bacteria specialized for mammal infection. The use of metagenomics helps us to understand evolution and physiology by elucidating the structure, function, and interactions of the major microbial communities, but it neglects the minority populations. Finally, Darwin’s descent with modification theory, as represented by the tree of life, no longer matches our current genomic knowledge because genomics has revealed the occurrence of de novo-created genes and the mosaic structure of genomes, the Rhizome of life is therefore more appropriate.

Array CGH phylogeny: how accurate are comparative genomic hybridization-based trees?

Background

Array-based Comparative Genomic Hybridization (CGH) data have been used to infer phylogenetic relationships. However, the reliability of array CGH analysis to determine evolutionary relationships has not been well established. In most CGH work, all species and strains are compared to a single reference species, whose genome was used to design the array. In the accompanying work, we critically evaluated CGH-based phylogeny using simulated competitive hybridization data. This work showed that a limited number of conditions, principally the tree topology and placement of the reference taxon in the tree, had a strong effect on the ability to recover the correct tree topology. Here, we add to our simulation study by testing the use of CGH as a phylogenetic tool with experimental CGH data from competitive hybridizations between N. crassa and other Neurospora species. In the discussion, we add to our empirical study of Neurospora by reanalyzing of data from a previous CGH phylogenetic analysis of the yeast sensu stricto complex.

Results

Array ratio data for Neurospora and related species were normalized with loess, robust spline, and linear ratio based methods, and then used to construct Neighbor-Joining and parsimony trees. These trees were compared to published phylogenetic analyses for Neurospora based on multilocus sequence analysis (MLSA). For the Neurospora dataset, the best combination of methods resulted in recovery of the MLSA tree topology less than half the time. Our reanalysis of a yeast dataset found that trees identical to established phylogeny were recovered only by pruning taxa - including the reference taxon - from the analysis.

Conclusion

Our results indicate that CGH data can be problematic for phylogenetic analysis. Success fluctuates based on the methods utilized to construct the tree and the taxa included. Selective pruning of the taxa improves the results - an impractical approach for normal phylogenetic analysis. From the more successful methods we make suggestions on the normalization and post-normalization methods that work best in estimating genetic distance between taxa.

Array Comparative Genomic Hybridizations: assessing the ability to recapture evolutionary relationships using an in silico approach

Background

Comparative Genomic Hybridization (CGH) with DNA microarrays has many biological applications including surveys of copy number changes in tumorogenesis, species detection and identification, and functional genomics studies among related organisms. Array CGH has also been used to infer phylogenetic relatedness among species or strains. Although the use of the entire genome can be seen as a considerable advantage for use in phylogenetic analysis, few such studies have questioned the reliability of array CGH to correctly determine evolutionary relationships. A potential flaw in this application lies in the fact that all comparisons are made to a single reference species. This situation differs from traditional DNA sequence, distance-based phylogenetic analyses where all possible pairwise comparisons are made for the isolates in question. By simulating array data based on the Neurospora crassa genome, we address this potential flaw and other questions regarding array CGH phylogeny.

Results

Our simulation data indicates that having a single reference can, in some cases, be a serious limitation when using this technique. Additionally, the tree building process with a single reference is sensitive to many factors including tree topology, choice of tree reconstruction method, and the distance metric used.

Conclusions

Without prior knowledge of the topology and placement of the reference taxon in the topology, the outcome is likely to be wrong and the error undetected. Given these limitations, using CGH to reveal phylogeny based on sequence divergence does not offer a robust alternative to traditional phylogenetic analysis.

SNP-based typing: a useful tool to study Bordetella pertussis populations

To monitor changes in Bordetella pertussis populations, mainly two typing methods are used; Pulsed-Field Gel Electrophoresis (PFGE) and Multiple-Locus Variable-Number Tandem Repeat Analysis (MLVA). In this study, a single nucleotide polymorphism (SNP) typing method, based on 87 SNPs, was developed and compared with PFGE and MLVA. The discriminatory indices of SNP typing, PFGE and MLVA were found to be 0.85, 0.95 and 0.83, respectively. Phylogenetic analysis, using SNP typing as Gold Standard, revealed false homoplasies in the PFGE and MLVA trees. Further, in contrast to the SNP-based tree, the PFGE- and MLVA-based trees did not reveal a positive correlation between root-to-tip distance and the isolation year of strains. Thus PFGE and MLVA do not allow an estimation of the relative age of the selected strains. In conclusion, SNP typing was found to be phylogenetically more informative than PFGE and more discriminative than MLVA. Further, in contrast to PFGE, it is readily standardized allowing interlaboratory comparisons. We applied SNP typing to study strains with a novel allele for the pertussis toxin promoter, ptxP3, which have a worldwide distribution and which have replaced the resident ptxP1 strains in the last 20 years. Previously, we showed that ptxP3 strains showed increased pertussis toxin expression and that their emergence was associated with increased notification in the Netherlands. SNP typing showed that the ptxP3 strains isolated in the Americas, Asia, Australia and Europe formed a monophyletic branch which recently diverged from ptxP1 strains. Two predominant ptxP3 SNP types were identified which spread worldwide. The widespread use of SNP typing will enhance our understanding of the evolution and global epidemiology of B. pertussis.

A novel sensor kinase is required for Bordetella bronchiseptica to colonize the lower respiratory tract

Bacterial virulence is influenced by the activity of two-component regulator systems (TCSs), which consist of membrane-bound sensor kinases that allow bacteria to sense the external environment and cytoplasmic, DNA-binding response regulator proteins that control appropriate gene expression. Respiratory pathogens of the Bordetella genus require the well-studied TCS BvgAS to control the expression of many genes required for colonization of the mammalian respiratory tract. Here we describe the identification of a novel gene in Bordetella bronchiseptica, plrS, the product of which shares sequence homology to several NtrY-family sensor kinases and is required for B. bronchiseptica to colonize and persist in the lower, but not upper, respiratory tract in rats and mice. The plrS gene is located immediately 5' to and presumably cotranscribed with a gene encoding a putative response regulator, supporting the idea that PlrS and the product of the downstream gene may compose a TCS. Consistent with this hypothesis, the PlrS-dependent colonization phenotype requires a conserved histidine that serves as the site of autophosphorylation in other sensor kinases, and in strains lacking plrS, the production and/or cellular localization of several immune-recognized proteins is altered in comparison to that in the wild-type strain. Because plrS is required for colonization and persistence only in the lower respiratory tract, a site where innate and adaptive immune mechanisms actively target infectious agents, we hypothesize that its role may be to allow Bordetella to resist the host immune response.

Detection of small RNAs in Bordetella pertussis and identification of a novel repeated genetic element

Background

Small bacterial RNAs (sRNAs) have been shown to participate in the regulation of gene expression and have been identified in numerous prokaryotic species. Some of them are involved in the regulation of virulence in pathogenic bacteria. So far, little is known about sRNAs in Bordetella, and only very few sRNAs have been identified in the genome of Bordetella pertussis, the causative agent of whooping cough.

Results

An in silico approach was used to predict sRNAs genes in intergenic regions of the B. pertussis genome. The genome sequences of B. pertussis, Bordetella parapertussis, Bordetella bronchiseptica and Bordetella avium were compared using a Blast, and significant hits were analyzed using RNAz. Twenty-three candidate regions were obtained, including regions encoding the already documented 6S RNA, and the GCVT and FMN riboswitches. The existence of sRNAs was verified by Northern blot analyses, and transcripts were detected for 13 out of the 20 additional candidates. These new sRNAs were named Bordetella pertussis RNAs, bpr. The expression of 4 of them differed between the early, exponential and late growth phases, and one of them, bprJ2, was found to be under the control of BvgA/BvgS two-component regulatory system of Bordetella virulence. A phylogenetic study of the bprJ sequence revealed a novel, so far undocumented repeat of ~90 bp, found in numerous copies in the Bordetella genomes and in that of other Betaproteobacteria. This repeat exhibits certain features of mobile elements.

Conclusion

We shown here that B. pertussis, like other pathogens, expresses sRNAs, and that the expression of one of them is controlled by the BvgA/BvgS system, similarly to most virulence genes, suggesting that it is involved in virulence of B. pertussis.

Comparative genomics of prevaccination and modern Bordetella pertussis strains

Background

Despite vaccination since the 1950s, pertussis has persisted and resurged. It remains a major cause of infant death worldwide and is the most prevalent vaccine-preventable disease in developed countries. The resurgence of pertussis has been associated with the expansion of Bordetella pertussis strains with a novel allele for the pertussis toxin (Ptx) promoter, ptxP3, which have replaced resident ptxP1 strains. Compared to ptxP1 strains, ptxP3 produce more Ptx resulting in increased virulence and immune suppression. To elucidate how B. pertussis has adapted to vaccination, we compared genome sequences of two ptxP3 strains with four strains isolated before and after the introduction vaccination.

Results

The distribution of SNPs in regions involved in transcription and translation suggested that changes in gene regulation play an important role in adaptation. No evidence was found for acquisition of novel genes. Modern strains differed significantly from prevaccination strains, both phylogenetically and with respect to particular alleles. The ptxP3 strains were found to have diverged recently from modern ptxP1 strains. Differences between ptxP3 and modern ptxP1 strains included SNPs in a number of pathogenicity-associated genes. Further, both gene inactivation and reactivation was observed in ptxP3 strains relative to modern ptxP1 strains.

Conclusions

Our work suggests that B. pertussis adapted by successive accumulation of SNPs and by gene (in)activation. In particular changes in gene regulation may have played a role in adaptation.

Antigenic Variation among Bordetella: Bordetella bronchiseptica strain MO149 expresses a novel o chain that is poorly immunogenic

The O chain polysaccharide (O PS) of Bordetella bronchiseptica and Bordetella parapertussis lipopolysaccharide is a homopolymer of 2,3-diacetamido-2,3-dideoxygalacturonic acid (GalNAc3NAcA) in which some of the sugars are present as uronamides. The terminal residue contains several unusual modifications. To date, two types of modification have been characterized, and a survey of numerous strains demonstrated that each contained one of these two modification types. Host antibody responses against the O PS are directed against the terminal residue modifications, and there is little cross-reactivity between the two types. This suggests that Bordetella O PS modifications represent a means of antigenic variation. Here we report the characterization of the O PS of B. bronchiseptica strain MO149. It consists of a novel two-sugar repeating unit and a novel terminal residue modification, with the structure Me-4-alpha-L-GalNAc3NAcA-(4-beta-D-GlcNAc3NAcA-4-alpha-L-GalNAc3NAcA-)(5-6)-, which we propose be defined as the B. bronchiseptica O3 PS. We show that the O3 PS is very poorly immunogenic and that the MO149 strain contains a novel wbm (O PS biosynthesis) locus. Thus, there is greater diversity among Bordetella O PSs than previously recognized, which is likely to be a result of selection pressure from host immunity. We also determine experimentally, for the first time, the absolute configuration of the diacetimido-uronic acid sugars in Bordetella O PS.

Evaluation of specificity of BP3385 for Bordetella pertussis detection

BP3385 has been proposed as a diagnostic PCR target for discriminating between Bordetella pertussis and other Bordetella species that also infect humans. Our results demonstrate that this gene is also present in some strains of Bordetella hinzii and Bordetella bronchiseptica.

Cross-species protection mediated by a Bordetella bronchiseptica strain lacking antigenic homologs present in acellular pertussis vaccines

The Bordetella species are Gram-negative bacterial pathogens that are characterized by long-term colonization of the mammalian respiratory tract and are causative agents of respiratory diseases in humans and animals. Despite widespread and efficient vaccination, there has been a world-wide resurgence of pertussis, which remains the leading cause of vaccine-preventable death in developed countries. It has been proposed that current acellular vaccines (Pa) composed of only a few bacterial proteins may be less efficacious because of vaccine-induced antigenic shifts and adaptations. To gain insight into the development of a newer generation of vaccines, we constructed a Bordetella bronchiseptica strain (LPaV) that does not express the antigenic homologs included in any of the Pa vaccines currently in use. This strain also lacks adenylate cyclase toxin, an essential virulence factor, and BipA, a surface protein. While LPaV colonized the mouse nose as efficiently as the wild-type strain, it was highly deficient in colonization of the lower respiratory tract and was attenuated in induction of inflammation and injury to the lungs. Strikingly, to our surprise, we found that in an intranasal murine challenge model, LPaV elicited cross-species protection against both B. bronchiseptica and Bordetella pertussis. Our data suggest the presence of immunogenic protective components other than those included in the pertussis vaccine. Combined with the whole-genome sequences of many Bordetella spp. that are available, the results of this study should serve as a platform for strategic development of the next generation of acellular pertussis vaccines.

Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents

Background

Rodents represent a high-risk reservoir for the emergence of new human pathogens. The recent completion of the 2.3 Mb genome of Bartonella grahamii, one of the most prevalent blood-borne bacteria in wild rodents, revealed a higher abundance of genes for host-cell interaction systems than in the genomes of closely related human pathogens. The sequence variability within the global B. grahamii population was recently investigated by multi locus sequence typing, but no study on the variability of putative host-cell interaction systems has been performed.

Results

To study the population dynamics of B. grahamii, we analyzed the genomic diversity on a whole-genome scale of 27 B. grahamii strains isolated from four different species of wild rodents in three geographic locations separated by less than 30 km. Even using highly variable spacer regions, only 3 sequence types were identified. This low sequence diversity contrasted with a high variability in genome content. Microarray comparative genome hybridizations identified genes for outer surface proteins, including a repeated region containing the fha gene for filamentous hemaggluttinin and a plasmid that encodes a type IV secretion system, as the most variable. The estimated generation times in liquid culture medium for a subset of strains ranged from 5 to 22 hours, but did not correlate with sequence type or presence/absence patterns of the fha gene or the plasmid.

Conclusion

Our study has revealed a geographic microstructure of B. grahamii in wild rodents. Despite near-identity in nucleotide sequence, major differences were observed in gene presence/absence patterns that did not segregate with host species. This suggests that genetically similar strains can infect a range of different hosts.

Genome sequence of the deep-rooted Yersinia pestis strain Angola reveals new insights into the evolution and pangenome of the plague bacterium

To gain insights into the origin and genome evolution of the plague bacterium Yersinia pestis, we have sequenced the deep-rooted strain Angola, a virulent Pestoides isolate. Its ancient nature makes this atypical isolate of particular importance in understanding the evolution of plague pathogenicity. Its chromosome features a unique genetic make-up intermediate between modern Y. pestis isolates and its evolutionary ancestor, Y. pseudotuberculosis. Our genotypic and phenotypic analyses led us to conclude that Angola belongs to one of the most ancient Y. pestis lineages thus far sequenced. The mobilome carries the first reported chimeric plasmid combining the two species-specific virulence plasmids. Genomic findings were validated in virulence assays demonstrating that its pathogenic potential is distinct from modern Y. pestis isolates. Human infection with this particular isolate would not be diagnosed by the standard clinical tests, as Angola lacks the plasmid-borne capsule, and a possible emergence of this genotype raises major public health concerns. To assess the genomic plasticity in Y. pestis, we investigated the global gene reservoir and estimated the pangenome at 4,844 unique protein-coding genes. As shown by the genomic analysis of this evolutionary key isolate, we found that the genomic plasticity within Y. pestis clearly was not as limited as previously thought, which is strengthened by the detection of the largest number of isolate-specific single-nucleotide polymorphisms (SNPs) currently reported in the species. This study identified numerous novel genetic signatures, some of which seem to be intimately associated with plague virulence. These markers are valuable in the development of a robust typing system critical for forensic, diagnostic, and epidemiological studies.

Changes in the genomic content of circulating Bordetella pertussis strains isolated from the Netherlands, Sweden, Japan and Australia: adaptive evolution or drift?

Background

Bordetella pertussis is the causative agent of human whooping cough (pertussis) and is particularly severe in infants. Despite worldwide vaccinations, whooping cough remains a public health problem. A significant increase in the incidence of whooping cough has been observed in many countries since the 1990s. Several reasons for the re-emergence of this highly contagious disease have been suggested. A particularly intriguing possibility is based on evidence indicating that pathogen adaptation may play a role in this process. In an attempt to gain insight into the genomic make-up of B. pertussis over the last 60 years, we used an oligonucleotide DNA microarray to compare the genomic contents of a collection of 171 strains of B. pertussis isolates from different countries.

Results

The CGH microarray analysis estimated the core genome of B. pertussis, to consist of 3,281 CDSs that are conserved among all B. pertussis strains, and represent 84.8% of all CDSs found in the 171 B. pertussis strains. A total of 64 regions of difference consisting of one or more contiguous CDSs were identified among the variable genes. CGH data also revealed that the genome size of B. pertussis strains is decreasing progressively over the past 60 years. Phylogenetic analysis of microarray data generated a minimum spanning tree that depicted the phylogenetic structure of the strains. B. pertussis strains with the same gene content were found in several different countries. However, geographic specificity of the B. pertussis strains was not observed. The gene content was determined to highly correlate with the ptxP-type of the strains.

Conclusions

An overview of genomic contents of a large collection of isolates from different countries allowed us to derive a core genome and a phylogenetic structure of B. pertussis. Our results show that B. pertussis is a dynamic organism that continues to evolve.

Bordetella pertussis and vaccination: the persistence of a genetically monomorphic pathogen

Before childhood vaccination was introduced in the 1950s, pertussis or whooping cough was a major cause of infant death worldwide. Widespread vaccination of children was successful in significantly reducing morbidity and mortality. However, despite vaccination, pertussis has persisted and, in the 1990s, resurged in a number of countries with highly vaccinated populations. Indeed, pertussis has become the most prevalent vaccine-preventable disease in developed countries with estimated infection frequencies of 1-6%. Recently vaccinated children are well protected against pertussis disease and its increase is mainly seen in adolescents and adults in which disease symptoms are often mild. The etiologic agent of pertussis, Bordetella pertussis, is extremely monomorphic and its ability to persist in the face of intensive vaccination is intriguing. Numerous studies have shown that B. pertussis populations changed after the introduction of vaccination suggesting adaptation. These adaptations did not involve the acquisition of novel genes but small genetic changes, mainly SNPs, and occurred in successive steps in a period of 40 years. The earliest adaptations resulted in antigenic divergence with vaccine strains. More recently, strains emerged with increased pertussis toxin (Ptx) production. Here I argue that the resurgence of pertussis is the compound effect of pathogen adaptation and waning immunity. I propose that the removal by vaccination of naïve infants as the major source for transmission was the crucial event which has driven the changes in B. pertussis populations. This has selected for strains which are more efficiently transmitted by primed hosts in which immunity has waned. The adaptation of B. pertussis to primed hosts involved delaying an effective immune response by antigenic divergence with vaccine strains and by increasing immune suppression through higher levels of Ptx production. Higher levels of Ptx may also benefit transmission by enhancing clinical symptoms. The study of B. pertussis populations has not only increased our understanding of pathogen evolution, but also suggests way to improve pertussis vaccines, underlining the public health significance of population-based studies of pathogens.