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

Acquisition and loss of virulence-associated factors during genome evolution and speciation in three clades of Bordetella species

Background

The genus Bordetella consists of nine species that include important respiratory pathogens such as the ‘classical’ species B. bronchiseptica, B. pertussis and B. parapertussis and six more distantly related and less extensively studied species. Here we analyze sequence diversity and gene content of 128 genome sequences from all nine species with focus on the evolution of virulence-associated factors.

Results

Both genome-wide sequence-based and gene content-based phylogenetic trees divide the genus into three species clades. The phylogenies are congruent between species suggesting genus-wide co-evolution of sequence diversity and gene content, but less correlated within species, mainly because of strain-specific presence of many different prophages. We compared the genomes with focus on virulence-associated genes and identified multiple clade-specific, species-specific and strain-specific events of gene acquisition and gene loss, including genes encoding O-antigens, protein secretion systems and bacterial toxins. Gene loss was more frequent than gene gain throughout the evolution, and loss of hundreds of genes was associated with the origin of several species, including the recently evolved human-restricted B. pertussis and B. holmesii, B. parapertussis and the avian pathogen B. avium.

Conclusions

Acquisition and loss of multiple genes drive the evolution and speciation in the genus Bordetella, including large scale gene loss associated with the origin of several species. Recent loss and functional inactivation of genes, including those encoding pertussis vaccine components and bacterial toxins, in individual strains emphasize ongoing evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3112-5) contains supplementary material, which is available to authorized users.

Pipeline for amplifying and analyzing amplicons of the V1-V3 region of the 16S rRNA gene

Background

Profiling of 16S rRNA gene sequences is an important tool for testing hypotheses in complex microbial communities, and analysis methods must be updated and validated as sequencing technologies advance. In host-associated bacterial communities, the V1–V3 region of the 16S rRNA gene is a valuable region to profile because it provides a useful level of taxonomic resolution; however, use of Illumina MiSeq data for experiments targeting this region needs validation.

Results

Using a MiSeq machine and the version 3 (300 × 2) chemistry, we sequenced the V1–V3 region of the 16S rRNA gene within a mock community. Nineteen bacteria and one archaeon comprised the mock community, and 12 replicate amplifications of the community were performed and sequenced. Sequencing the large fragment (490 bp) that encompasses V1–V3 yielded a higher error rate (3.6 %) than has been reported when using smaller fragment sizes. This higher error rate was due to a large number of sequences that occurred only one or two times among all mock community samples. Removing sequences that occurred one time among all samples (singletons) reduced the error rate to 1.4 %. Diversity estimates of the mock community containing all sequences were inflated, whereas estimates following singleton removal more closely reflected the actual mock community membership. A higher percentage of the sequences could be taxonomically assigned after singleton and doubleton sequences were removed, and the assignments reflected the membership of the input DNA.

Conclusions

Sequencing the V1–V3 region of the 16S rRNA gene on the MiSeq platform may require additional sequence curation in silico, and improved error rates and diversity estimates show that removing low-frequency sequences is reasonable. When datasets have a high number of singletons, these singletons can be removed from the analysis without losing statistical power while reducing error and improving microbiota assessment.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-2172-6) contains supplementary material, which is available to authorized users.

An Integrated Metabolomic and Genomic Mining Workflow To Uncover the Biosynthetic Potential of Bacteria

Microorganisms are a rich source of bioactives; however, chemical identification is a major bottleneck. Strategies that can prioritize the most prolific microbial strains and novel compounds are of great interest. Here, we present an integrated approach to evaluate the biosynthetic richness in bacteria and mine the associated chemical diversity. Thirteen strains closely related to Pseudoalteromonas luteoviolacea isolated from all over the Earth were analyzed using an untargeted metabolomics strategy, and metabolomic profiles were correlated with whole-genome sequences of the strains. We found considerable diversity: only 2% of the chemical features and 7% of the biosynthetic genes were common to all strains, while 30% of all features and 24% of the genes were unique to single strains. The list of chemical features was reduced to 50 discriminating features using a genetic algorithm and support vector machines. Features were dereplicated by tandem mass spectrometry (MS/MS) networking to identify molecular families of the same biosynthetic origin, and the associated pathways were probed using comparative genomics. Most of the discriminating features were related to antibacterial compounds, including the thiomarinols that were reported from P. luteoviolacea here for the first time. By comparative genomics, we identified the biosynthetic cluster responsible for the production of the antibiotic indolmycin, which could not be predicted with standard methods. In conclusion, we present an efficient, integrative strategy for elucidating the chemical richness of a given set of bacteria and link the chemistry to biosynthetic genes. IMPORTANCE We here combine chemical analysis and genomics to probe for new bioactive secondary metabolites based on their pattern of distribution within bacterial species. We demonstrate the usefulness of this combined approach in a group of marine Gram-negative bacteria closely related to Pseudoalteromonas luteoviolacea, which is a species known to produce a broad spectrum of chemicals. The approach allowed us to identify new antibiotics and their associated biosynthetic pathways. Combining chemical analysis and genetics is an efficient "mining" workflow for identifying diverse pharmaceutical candidates in a broad range of microorganisms and therefore of great use in bioprospecting.

An Integrated Metabolomic and Genomic Mining Workflow To Uncover the Biosynthetic Potential of Bacteria

Microorganisms are a rich source of bioactives; however, chemical identification is a major bottleneck. Strategies that can prioritize the most prolific microbial strains and novel compounds are of great interest. Here, we present an integrated approach to evaluate the biosynthetic richness in bacteria and mine the associated chemical diversity. Thirteen strains closely related to Pseudoalteromonas luteoviolacea isolated from all over the Earth were analyzed using an untargeted metabolomics strategy, and metabolomic profiles were correlated with whole-genome sequences of the strains. We found considerable diversity: only 2% of the chemical features and 7% of the biosynthetic genes were common to all strains, while 30% of all features and 24% of the genes were unique to single strains. The list of chemical features was reduced to 50 discriminating features using a genetic algorithm and support vector machines. Features were dereplicated by tandem mass spectrometry (MS/MS) networking to identify molecular families of the same biosynthetic origin, and the associated pathways were probed using comparative genomics. Most of the discriminating features were related to antibacterial compounds, including the thiomarinols that were reported from P. luteoviolacea here for the first time. By comparative genomics, we identified the biosynthetic cluster responsible for the production of the antibiotic indolmycin, which could not be predicted with standard methods. In conclusion, we present an efficient, integrative strategy for elucidating the chemical richness of a given set of bacteria and link the chemistry to biosynthetic genes. IMPORTANCE We here combine chemical analysis and genomics to probe for new bioactive secondary metabolites based on their pattern of distribution within bacterial species. We demonstrate the usefulness of this combined approach in a group of marine Gram-negative bacteria closely related to Pseudoalteromonas luteoviolacea, which is a species known to produce a broad spectrum of chemicals. The approach allowed us to identify new antibiotics and their associated biosynthetic pathways. Combining chemical analysis and genetics is an efficient "mining" workflow for identifying diverse pharmaceutical candidates in a broad range of microorganisms and therefore of great use in bioprospecting.

Whole genome sequencing revealed host adaptation-focused genomic plasticity of pathogenic Leptospira

Leptospirosis, caused by pathogenic Leptospira spp., has recently been recognized as an emerging infectious disease worldwide. Despite its severity and global importance, knowledge about the molecular pathogenesis and virulence evolution of Leptospira spp. remains limited. Here we sequenced and analyzed 102 isolates representing global sources. A high genomic variability were observed among different Leptospira species, which was attributed to massive gene gain and loss events allowing for adaptation to specific niche conditions and changing host environments. Horizontal gene transfer and gene duplication allowed the stepwise acquisition of virulence factors in pathogenic Leptospira evolved from a recent common ancestor. More importantly, the abundant expansion of specific virulence-related protein families, such as metalloproteases-associated paralogs, were exclusively identified in pathogenic species, reflecting the importance of these protein families in the pathogenesis of leptospirosis. Our observations also indicated that positive selection played a crucial role on this bacteria adaptation to hosts. These novel findings may lead to greater understanding of the global diversity and virulence evolution of Leptospira spp.

What Makes a Bacterial Species Pathogenic?:Comparative Genomic Analysis of the Genus Leptospira

Leptospirosis, caused by spirochetes of the genus Leptospira, is a globally widespread, neglected and emerging zoonotic disease. While whole genome analysis of individual pathogenic, intermediately pathogenic and saprophytic Leptospira species has been reported, comprehensive cross-species genomic comparison of all known species of infectious and non-infectious Leptospira, with the goal of identifying genes related to pathogenesis and mammalian host adaptation, remains a key gap in the field. Infectious Leptospira, comprised of pathogenic and intermediately pathogenic Leptospira, evolutionarily diverged from non-infectious, saprophytic Leptospira, as demonstrated by the following computational biology analyses: 1) the definitive taxonomy and evolutionary relatedness among all known Leptospira species; 2) genomically-predicted metabolic reconstructions that indicate novel adaptation of infectious Leptospira to mammals, including sialic acid biosynthesis, pathogen-specific porphyrin metabolism and the first-time demonstration of cobalamin (B12) autotrophy as a bacterial virulence factor; 3) CRISPR/Cas systems demonstrated only to be present in pathogenic Leptospira, suggesting a potential mechanism for this clade's refractoriness to gene targeting; 4) finding Leptospira pathogen-specific specialized protein secretion systems; 5) novel virulence-related genes/gene families such as the Virulence Modifying (VM) (PF07598 paralogs) proteins and pathogen-specific adhesins; 6) discovery of novel, pathogen-specific protein modification and secretion mechanisms including unique lipoprotein signal peptide motifs, Sec-independent twin arginine protein secretion motifs, and the absence of certain canonical signal recognition particle proteins from all Leptospira; and 7) and demonstration of infectious Leptospira-specific signal-responsive gene expression, motility and chemotaxis systems. By identifying large scale changes in infectious (pathogenic and intermediately pathogenic) vs. non-infectious Leptospira, this work provides new insights into the evolution of a genus of bacterial pathogens. This work will be a comprehensive roadmap for understanding leptospirosis pathogenesis. More generally, it provides new insights into mechanisms by which bacterial pathogens adapt to mammalian hosts.

Genomic dissection of Australian Bordetella pertussis isolates from the 2008-2012 epidemic

OBJECTIVES

Despite high pertussis vaccination coverage, Australia experienced a prolonged epidemic in 2008-2012. The predominant Bordetella pertussis genotype harboured pertussis toxin promoter allele, ptxP3, and pertactin gene allele, prn2. The emergence and expansion of prn non-expressing isolates (Prn negative), were also observed. We aimed to investigate the microevolution and genomic diversity of epidemic B. pertussis isolates.

METHODS

We sequenced 22 B. pertussis isolates collected in 2008-2012 from two states of Australia which are geographically widely separated. Ten of the 22 were Prn negative isolates with three different modes of silencing of prn (prn::IS481F, prn::IS481R and prn::IS1002). Five pre-epidemic isolates were also sequenced for comparison.

RESULTS

Five single nucleotide polymorphisms were common in the epidemic isolates and differentiated them from pre-epidemic isolates. The Australian epidemic isolates can be divided into five lineages (EL1-EL5) with EL1 containing only Prn negative isolates. Comparison with global isolates showed that three lineages remained geographically and temporally distinct whereas two lineages mixed with isolates from 2012 UK outbreak.

CONCLUSION

Our results suggest significant diversification and the microevolution of B. pertussis within the 2008-2012 Australian epidemic.

Transcriptional profiling of Bordetella pertussis reveals requirement of RNA chaperone Hfq for Type III secretion system functionality

Bordetella pertussis, the causative agent of human whooping cough (pertussis) produces a complex array of virulence factors in order to establish efficient infection in the host. The RNA chaperone Hfq and small regulatory RNAs are key players in posttranscriptional regulation in bacteria and have been shown to play an essential role in virulence of a broad spectrum of bacterial pathogens. This study represents the first attempt to characterize the Hfq regulon of the human pathogen B. pertussis under laboratory conditions as well as upon passage in the host and indicates that loss of Hfq has a profound effect on gene expression in B. pertussis. Comparative transcriptional profiling revealed that Hfq is required for expression of several virulence factors in B. pertussis cells including the Type III secretion system (T3SS). In striking contrast to the wt strain, T3SS did not become operational in the hfq mutant passaged either through mice or macrophages thereby proving that Hfq is required for the functionality of the B. pertussis T3SS. Likewise, expression of virulence factors vag8 and tcfA encoding autotransporter and tracheal colonization factor, respectively, was strongly reduced in the hfq mutant. Importantly, for the first time we demonstrate that B. pertussis T3SS can be activated upon contact with macrophage cells in vitro.

Complete Genome Sequences of 11 Bordetella pertussis Strains Representing the Pandemic ptxP3 Lineage

Pathogen adaptation has contributed to the resurgence of pertussis. To facilitate our understanding of this adaptation we report here 11 completely closed and annotated Bordetella pertussis genomes representing the pandemic ptxP3 lineage. Our analyses included six strains which do not produce the vaccine components pertactin and/or filamentous hemagglutinin.

Comparative genomic analysis of the swine pathogen Bordetella bronchisepticastrain KM22

The well-characterized Bordetella bronchiseptica strain KM22, originally isolated from a pig with atrophic rhinitis, has been used to develop a reproducible swine respiratory disease model. The goal of this study was to identify genetic features unique to KM22 by comparing the genome sequence of KM22 to the laboratory reference strain RB50. To gain a broader perspective of the genetic relationship of KM22 among other B. bronchiseptica strains, selected genes of KM22 were then compared to five other B. bronchiseptica strains isolated from different hosts. Overall, the KM22 genome sequence is more similar to the genome sequences of the strains isolated from animals than the strains isolated from humans. The majority of virulence gene expression in Bordetella is positively regulated by the two-component sensory transduction system BvgAS. bopN, bvgA, fimB, and fimC were the most highly conserved BvgAS-regulated genes present in all seven strains analyzed. In contrast, the BvgAS-regulated genes present in all seven strains with the highest sequence divergence werefimN, fim2, fhaL, andfhaS. A total of eight major fimbrial subunit genes were identified in KM22. Quantitative real-time PCR data demonstrated that seven of the eight fimbrial subunit genes identified in KM22 are expressed and regulated by BvgAS. The annotation of the KM22 genome sequence, coupled with the comparative genomic analyses reported in this study, can be used to facilitate the development of vaccines with improved efficacy towards B. bronchiseptica in swine to decrease the prevalence and disease burden caused by this pathogen.

Type Six Secretion System of Bordetella bronchiseptica and Adaptive Immune Components Limit Intracellular Survival During Infection

The Type Six Secretion System (T6SS) is required for Bordetella bronchiseptica cytotoxicity, cytokine modulation, infection, and persistence. However, one-third of recently sequenced Bordetella bronchiseptica strains of the predominantly human-associated Complex IV have lost their T6SS through gene deletion or degradation. Since most human B. bronchiseptica infections occur in immunocompromised patients, we determine here whether loss of Type Six Secretion is beneficial to B. bronchiseptica during infection of immunocompromised mice. Infection of mice lacking adaptive immunity (Rag1-/- mice) with a T6SS-deficient mutant results in a hypervirulent phenotype that is characterized by high numbers of intracellular bacteria in systemic organs. In contrast, wild-type B. bronchiseptica kill their eukaryotic cellular hosts via a T6SS-dependent mechanism that prevents survival in systemic organs. High numbers of intracellular bacteria recovered from immunodeficient mice but only low numbers from wild-type mice demonstrates that B. bronchiseptica survival in an intracellular niche is limited by B and T cell responses. Understanding the nature of intracellular survival during infection, and its effects on the generation and function of the host immune response, are important to contain and control the spread of Bordetella-caused disease.

Diversity of secretion systems associated with virulence characteristics of the classical bordetellae

Secretion systems are key virulence factors, modulating interactions between pathogens and the host's immune response. Six potential secretion systems (types 1-6; T1SS-T6SS) have been discussed in classical bordetellae, respiratory commensals/pathogens of mammals. The prototypical Bordetella bronchiseptica strain RB50 genome seems to contain all six systems, whilst two human-restricted subspecies, Bordetella parapertussis and Bordetella pertussis, have lost different subsets of these. This implicates secretion systems in the divergent evolutionary histories that have led to their success in different niches. Based on our previous work demonstrating that changes in secretion systems are associated with virulence characteristics, we hypothesized there would be substantial divergence of the loci encoding each amongst sequenced strains. Here, we describe extensive differences in secretion system loci; 10 of the 11 sequenced strains had lost subsets of genes or one entire secretion system locus. These loci contained genes homologous to those present in the respective loci in distantly related organisms, as well as genes unique to bordetellae, suggesting novel and/or auxiliary functions. The high degree of conservation of the T3SS locus, a complex machine with interdependent parts that must be conserved, stands in dramatic contrast to repeated loss of T5aSS 'autotransporters', which function as an autonomous unit. This comparative analysis provided insights into critical aspects of each pathogen's adaptation to its different niche, and the relative contributions of recombination, mutation and horizontal gene transfer. In addition, the relative conservation of various secretion systems is an important consideration in the ongoing search for more highly conserved protective antigens for the next generation of pertussis vaccines.

Contribution of pertussis toxin to the pathogenesis of pertussis disease

Pertussis toxin (PT) is a multisubunit protein toxin secreted by Bordetella pertussis, the bacterial agent of the disease pertussis or whooping cough. PT in detoxified form is a component of all licensed acellular pertussis vaccines, since it is considered to be an important virulence factor for this pathogen. PT inhibits G protein-coupled receptor signaling through Gi proteins in mammalian cells, an activity that has led to its widespread use as a cell biology tool. But how does this activity of PT contribute to pertussis, including the severe respiratory symptoms of this disease? In this minireview, the contribution of PT to the pathogenesis of pertussis disease will be considered based on evidence from both human infections and animal model studies. Although definitive proof of the role of PT in humans is lacking, substantial evidence supports the idea that PT is a major contributor to pertussis pathology, including the severe respiratory symptoms associated with this disease.

New Data on Vaccine Antigen Deficient Bordetella pertussis Isolates

Evolution of Bordetella pertussis is driven by natural and vaccine pressures. Isolates circulating in regions with high vaccination coverage present multiple allelic and antigenic variations as compared to isolates collected before introduction of vaccination. Furthermore, during the last epidemics reported in regions using pertussis acellular vaccines, isolates deficient for vaccine antigens, such as pertactin (PRN), were reported to reach high proportions of circulating isolates. More sporadic filamentous hemagglutinin (FHA) or pertussis toxin (PT) deficient isolates were also collected. The whole genome of some recent French isolates, deficient or non-deficient in vaccine antigens, were analyzed. Transcription profiles of the expression of the main virulence factors were also compared. The invasive phenotype in an in vitro human tracheal epithelial (HTE) cell model of infection was evaluated. Our genomic analysis focused on SNPs related to virulence genes known to be more likely to present allelic polymorphism. Transcriptomic data indicated that isolates circulating since the introduction of pertussis vaccines present lower transcription levels of the main virulence genes than the isolates of the pre-vaccine era. Furthermore, isolates not producing FHA present significantly higher expression levels of the entire set of genes tested. Finally, we observed that recent isolates are more invasive in HTE cells when compared to the reference strain, but no multiplication occurs within cells.

Bordetella pertussis transmission

Bordetella pertussis and B. bronchiseptica are Gram-negative bacterial respiratory pathogens. Bordetella pertussis is the causative agent of whooping cough and is considered a human-adapted variant of B. bronchiseptica. Bordetella pertussis and B. bronchiseptica share mechanisms of pathogenesis and are genetically closely related. However, despite the close genetic relatedness, these Bordetella species differ in several classic fundamental aspects of bacterial pathogens such as host range, pathologies and persistence. The development of the baboon model for the study of B. pertussis transmission, along with the development of the swine and mouse model for the study of B. bronchiseptica, has enabled the investigation of different aspects of transmission including the route, attack rate, role of bacterial and host factors, and the impact of vaccination on transmission. This review will focus on B. pertussis transmission and how animal models of B. pertussis transmission and transmission models using the closely related B. bronchiseptica have increased our understanding of B. pertussis transmission.

Whole-genome sequencing reveals the effect of vaccination on the evolution of Bordetella pertussis

Herd immunity can potentially induce a change of circulating viruses. However, it remains largely unknown that how bacterial pathogens adapt to vaccination. In this study, Bordetella pertussis, the causative agent of whooping cough, was selected as an example to explore possible effect of vaccination on the bacterial pathogen. We sequenced and analysed the complete genomes of 40 B. pertussis strains from Finland and China, as well as 11 previously sequenced strains from the Netherlands, where different vaccination strategies have been used over the past 50 years. The results showed that the molecular clock moved at different rates in these countries and in distinct periods, which suggested that evolution of the B. pertussis population was closely associated with the country vaccination coverage. Comparative whole-genome analyses indicated that evolution in this human-restricted pathogen was mainly characterised by ongoing genetic shift and gene loss. Furthermore, 116 SNPs were specifically detected in currently circulating ptxP3-containing strains. The finding might explain the successful emergence of this lineage and its spread worldwide. Collectively, our results suggest that the immune pressure of vaccination is one major driving force for the evolution of B. pertussis, which facilitates further exploration of the pathogenicity of B. pertussis.

A novel method of consensus pan-chromosome assembly and large-scale comparative analysis reveal the highly flexible pan-genome of Acinetobacter baumannii

Background

Infections by pan-drug resistant Acinetobacter baumannii plague military and civilian healthcare systems. Previous A. baumannii pan-genomic studies used modest sample sizes of low diversity and comparisons to a single reference genome, limiting our understanding of gene order and content. A consensus representation of multiple genomes will provide a better framework for comparison. A large-scale comparative study will identify genomic determinants associated with their diversity and adaptation as a successful pathogen.

Results

We determine draft-level genomic sequence of 50 diverse military isolates and conduct the largest bacterial pan-genome analysis of 249 genomes. The pan-genome of A. baumannii is open when the input genomes are normalized for diversity with 1867 core proteins and a paralog-collapsed pan-genome size of 11,694 proteins. We developed a novel graph-based algorithm and use it to assemble the first consensus pan-chromosome, identifying both the order and orientation of core genes and flexible genomic regions. Comparative genome analyses demonstrate the existence of novel resistance islands and isolates with increased numbers of resistance island insertions over time, from single insertions in the 1950s to triple insertions in 2011. Gene clusters responsible for carbon utilization, siderophore production, and pilus assembly demonstrate frequent gain or loss among isolates.

Conclusions

The highly variable and dynamic nature of the A. baumannii genome may be the result of its success in rapidly adapting to both abiotic and biotic environments through the gain and loss of gene clusters controlling fitness. Importantly, some archaic adaptation mechanisms appear to have reemerged among recent isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0701-6) contains supplementary material, which is available to authorized users.

Cooperative roles for fimbria and filamentous hemagglutinin in Bordetella adherence and immune modulation

Bordetella fimbriae (FIM) are generally considered to function as adhesins despite a lack of experimental evidence supporting this conclusion for Bordetella pertussis and evidence against a requirement for FIM in adherence of Bordetella bronchiseptica to mammalian cell lines. Using B. bronchiseptica and mice, we developed an in vivo adherence assay that revealed that FIM do function as critically important adhesins in the lower respiratory tract. In the first few days postinoculation, FIM-deficient B. bronchiseptica induced a more robust inflammatory response than wild-type bacteria did, suggesting that FIM, like filamentous hemagglutinin (FHA), allow B. bronchiseptica to suppress the innate immune response to infection. Localization analyses indicated that FIM are required for efficient attachment to airway epithelium, as bacteria lacking FIM localized to alveoli. FHA-deficient bacteria, in contrast, localized to airways. Bacteria unable to produce both FIM and FHA localized to alveoli and caused increased inflammation and histopathology identical to that caused by FIM-deficient bacteria, demonstrating that lack of FIM is epistatic to lack of FHA. Coinoculation experiments provided evidence that wild-type B. bronchiseptica suppresses inflammation locally within the respiratory tract and that both FHA and FIM are required for defense against clearance by the innate immune system. Altogether, our data suggest that FIM-mediated adherence to airway epithelium is a critical first step in Bordetella infection that allows FHA-dependent interactions to mediate tight adherence, suppression of inflammation, and resistance to inflammatory cell-mediated clearance. Our results suggest that mucosal antibodies capable of blocking FIM-mediated interactions could prevent bacterial colonization of the lower respiratory tract.

Although fimbriae (FIM) have been shown to be important mediators of adherence for many bacterial pathogens, there is surprisingly little experimental evidence supporting this role for Bordetella fimbria. Our results provide the first demonstration that Bordetella FIM function as adhesins in vivo, specifically to airway epithelium. Furthermore, our results suggest that FIM mediate initial interactions with airway epithelial cells that are followed by tight filamentous hemagglutinin (FHA)-mediated binding and that together, FIM and FHA allow Bordetella to suppress inflammation, leading to prolonged colonization. Given the shortcoming of the current acellular component pertussis (aP) vaccine in preventing colonization, these findings suggest that generation of antibodies capable of blocking FIM-mediated adherence could potentially prevent Bordetella colonization.

Draft Genome Sequences of 53 Genetically Distinct Isolates of Bordetella bronchiseptica Representing 11 Terrestrial and Aquatic Hosts

Bordetella bronchiseptica infects a variety of mammalian and avian hosts. Here, we report the genome sequences of 53 genetically distinct isolates acquired from a broad range of terrestrial and aquatic animals. These data will greatly facilitate ongoing efforts to better understand the evolution, host adaptation, and virulence mechanisms of B. bronchiseptica.

Novel, host-restricted genotypes of Bordetella bronchiseptica associated with phocine respiratory tract isolates

During a succession of phocine morbillivirus outbreaks spanning the past 25 years, Bordetella bronchiseptica was identified as a frequent secondary invader and cause of death. The goal of this study was to evaluate genetic diversity and the molecular basis for host specificity among seal isolates from these outbreaks. MLST and PvuII ribotyping of 54 isolates from Scottish, English or Danish coasts of the Atlantic or North Sea revealed a single, host-restricted genotype. A single, novel genotype, unique from that of the Atlantic and North Sea isolates, was found in isolates from an outbreak in the Caspian Sea. Phylogenetic analysis based either on MLST sequence, ribotype patterns or genome-wide SNPs consistently placed both seal-specific genotypes within the same major clade but indicates a distinct evolutionary history for each. An additional isolate from the intestinal tract of a seal on the south-west coast of England has a genotype otherwise found in rabbit, guinea pig and pig isolates. To investigate the molecular basis for host specificity, DNA and predicted protein sequences of virulence genes that mediate host interactions were used in comparisons between a North Sea isolate, a Caspian Sea isolate and each of their closest relatives as inferred from genome-wide SNP analysis. Despite their phylogenetic divergence, fewer nucleotide and amino acid substitutions were found in comparisons of the two seal isolates than in comparisons with closely related strains. These data indicate isolates of B. bronchiseptica associated with respiratory disease in seals comprise unique, host-adapted and highly clonal populations.

A brief review of software tools for pangenomics

Since the proposal for pangenomic study, there have been a dozen software tools actively in use for pangenomic analysis. By the end of 2014, Panseq and the pan-genomes analysis pipeline (PGAP) ranked as the top two most popular packages according to cumulative citations of peer-reviewed scientific publications. The functions of the software packages and tools, albeit variable among them, include categorizing orthologous genes, calculating pangenomic profiles, integrating gene annotations, and constructing phylogenies. As epigenomic elements are being gradually revealed in prokaryotes, it is expected that pangenomic databases and toolkits have to be extended to handle information of detailed functional annotations for genes and non-protein-coding sequences including non-coding RNAs, insertion elements, and conserved structural elements. To develop better bioinformatic tools, user feedback and integration of novel features are both of essence.

High correlation between genotypes and phenotypes of environmental bacteria Comamonas testosteroni strains

BACKGROUND

Members of Comamonas testosteroni are environmental microorganisms that are usually found in polluted environment samples. They utilize steroids and aromatic compounds but rarely sugars, and show resistance to multiple heavy metals and multiple drugs. However, comprehensive genomic analysis among the C. testosteroni strains is lacked.

RESULTS

To understand the genome bases of the features of C. testosteroni, we sequenced 10 strains of this species and analyzed them together with other related published genome sequences. The results revealed that: 1) the strains of C. testosteroni have genome sizes ranging from 5.1 to 6.0 Mb and G + C contents ranging from 61.1% to 61.8%. The pan-genome contained 10,165 gene families and the core genome contained 3,599 gene families. Heap's law analysis indicated that the pan-genome of C. testosteroni may be open (α = 0.639); 2) by analyzing 31 phenotypes of 11 available C. testosteroni strains, 99.4% of the genotypes (putative genes) were found to be correlated to the phenotypes, indicating a high correlation between phenotypes and genotypes; 3) gene clusters for nitrate reduction, steroids degradation and metal and multi-drug resistance were found and were highly conserved among all the genomes of this species; 4) the genome similarity of C. testosteroni may be related to the geographical distances.

CONCLUSIONS

This work provided an overview on the genomes of C. testosteroni and new genome resources that would accelerate the further investigations of this species. Importantly, this work focused on the analysis of potential genetic determinants for the typical characters and found high correlation between the phenotypes and their corresponding genotypes.

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.

Heterogeneity of Bordetella bronchiseptica adenylate cyclase (cyaA) RTX domain

Bordetella bronchiseptica is a widespread pathogen, with a broad host range, occasionally including humans. Diverse virulence factors (adhesins, toxins) allow its adaptation to its host, but this property of the adenylate cyclase (cyaA) toxin is not well understood. In this study, we analyzed the repeats-in-toxin domain of B. bronchiseptica cyaA with PCR, followed by restriction fragment length analysis. Of ninety-two B. bronchiseptica strains collected from different hosts and geographic regions, 72 (78.3 %) carried cyaA and four RFLP types (A-D) were established using NarI and SalI. However, in 20 strains, cyaA was replaced with a peptide transport protein operon. A phylogenetic tree based on partial nucleotide sequences of cyaA revealed that group 2 contains strains of specifically human origin, whereas subgroup 1a contains all but one of the strains from pigs. The human strains showed many PCR-RFLP and sequence variants, confirming the clonal population structure of B. bronchiseptica.

Whole-cell pertussis vaccine potency assays: the Kendrick test and alternative assays

Whole-cell pertussis vaccines are still widely used across the globe and have been shown to produce longer lasting immunity against pertussis infection than acellular pertussis vaccines. Therefore, whole-cell vaccines are likely to continue to be used for the foreseeable future. The intracerebral mouse protection test (Kendrick test) is effective for determining the potency of whole-cell pertussis vaccines and is the only test that has shown a correlation with protection in children. Here we review the Kendrick test in terms of international requirements for vaccine potency and critical technical points to be considered for a successful test including test validity, in-house references and statistical analysis. There are objections to the Kendrick test on animal welfare and technical grounds. Respiratory challenge assays, nitric oxide induction assay and serological assays have been developed and have been proposed as possible methods which might provide alternatives to the Kendrick test. These methods and their limitations are also briefly discussed. Establishment of validated in vitro correlates of protection has yet to be achieved. New technical developments, such as genome sequence and the use of gene microarrays to screen responses triggered by vaccine components may also provide leads to alternative assays to the Kendrick test by identifying biomarkers of protection.

Other Bordetellas, lessons for and from pertussis vaccines

The Bordetella genus comprises nine species of which Bordetella pertussis and B. parapertussis are isolated from humans and are the most studied Bordetella species since they cause whooping cough. They both originate from B. bronchiseptica, which infects several mammals and immune compromised humans, but the intensive use of pertussis vaccines induced changes in B. pertussis and B. parapertussis populations. B. petrii and B. holmesii are other species of unknown reservoir and transmission pattern that have been described in humans. It is still unknown whether these species are pathogens for humans or only opportunistic bacteria but biological diagnosis has confirmed the presence of B. holmesii in human respiratory samples while B. petrii and the four other species have little implications for public health.

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.

Global population structure and evolution of Bordetella pertussis and their relationship with vaccination

Bordetella pertussis causes pertussis, a respiratory disease that is most severe for infants. Vaccination was introduced in the 1950s, and in recent years, a resurgence of disease was observed worldwide, with significant mortality in infants. Possible causes for this include the switch from whole-cell vaccines (WCVs) to less effective acellular vaccines (ACVs), waning immunity, and pathogen adaptation. Pathogen adaptation is suggested by antigenic divergence between vaccine strains and circulating strains and by the emergence of strains with increased pertussis toxin production. We applied comparative genomics to a worldwide collection of 343 B. pertussis strains isolated between 1920 and 2010. The global phylogeny showed two deep branches; the largest of these contained 98% of all strains, and its expansion correlated temporally with the first descriptions of pertussis outbreaks in Europe in the 16th century. We found little evidence of recent geographical clustering of the strains within this lineage, suggesting rapid strain flow between countries. We observed that changes in genes encoding proteins implicated in protective immunity that are included in ACVs occurred after the introduction of WCVs but before the switch to ACVs. Furthermore, our analyses consistently suggested that virulence-associated genes and genes coding for surface-exposed proteins were involved in adaptation. However, many of the putative adaptive loci identified have a physiological role, and further studies of these loci may reveal less obvious ways in which B. pertussis and the host interact. This work provides insight into ways in which pathogens may adapt to vaccination and suggests ways to improve pertussis vaccines. IMPORTANCE Whooping cough is mainly caused by Bordetella pertussis, and current vaccines are targeted against this organism. Recently, there have been increasing outbreaks of whooping cough, even where vaccine coverage is high. Analysis of the genomes of 343 B. pertussis isolates from around the world over the last 100 years suggests that the organism has emerged within the last 500 years, consistent with historical records. We show that global transmission of new strains is very rapid and that the worldwide population of B. pertussis is evolving in response to vaccine introduction, potentially enabling vaccine escape.

Bordetella pertussis pathogenesis: current and future challenges

Pertussis, also known as whooping cough, has recently re-emerged as a major public health threat despite high levels of vaccination against the aetiological agent Bordetella pertussis. In this Review, we describe the pathogenesis of this disease, with a focus on recent mechanistic insights into B. pertussis virulence-factor function. We also discuss the changing epidemiology of pertussis and the challenges facing vaccine development. Despite decades of research, many aspects of B. pertussis physiology and pathogenesis remain poorly understood. We highlight knowledge gaps that must be addressed to develop improved vaccines and therapeutic strategies.

Comparative analyses of a cystic fibrosis isolate of Bordetella bronchiseptica reveal differences in important pathogenic phenotypes

Bordetella bronchiseptica is a Gram-negative bacterium that infects and causes disease in a wide variety of animals. B. bronchiseptica also infects humans, thereby demonstrating zoonotic transmission. An extensive characterization of human B. bronchiseptica isolates is needed to better understand the distinct genetic and phenotypic traits associated with these zoonotic transmission events. Using whole-genome transcriptome and CGH analysis, we report that a B. bronchiseptica cystic fibrosis isolate, T44625, contains a distinct genomic content of virulence-associated genes and differentially expresses these genes compared to the sequenced model laboratory strain RB50, a rabbit isolate. The differential gene expression pattern correlated with unique phenotypes exhibited by T44625, which included lower motility, increased aggregation, hyperbiofilm formation, and an increased in vitro capacity to adhere to respiratory epithelial cells. Using a mouse intranasal infection model, we found that although defective in establishing high bacterial burdens early during the infection process, T44625 persisted efficiently in the mouse nose. By documenting the unique genomic and phenotypic attributes of T44625, this report provides a blueprint for understanding the successful zoonotic potential of B. bronchiseptica and other zoonotic bacteria.

Horizontally acquired divergent O-antigen contributes to escape from cross-immunity in the classical bordetellae

Background

Horizontal gene transfer (HGT) allows for rapid spread of genetic material between species, increasing genetic and phenotypic diversity. Although HGT contributes to adaptation and is widespread in many bacteria, others show little HGT. This study builds on previous work to analyze the evolutionary mechanisms contributing to variation within the locus encoding a prominent antigen of the classical bordetellae.

Results

We observed amongst classical bordetellae discrete regions of the lipopolysaccharide O-antigen locus with higher sequence diversity than the genome average. Regions of this locus had less than 50% sequence similarity, low dN/dS ratios and lower GC content compared to the genome average. Additionally, phylogenetic tree topologies based on genome-wide SNPs were incongruent with those based on genes within these variable regions, suggesting portions of the O-antigen locus may have been horizontally transferred. Furthermore, several predicted recombination breakpoints correspond with the ends of these variable regions. To examine the evolutionary forces that might have selected for this rare example of HGT in bordetellae, we compared in vitro and in vivo phenotypes associated with different O-antigen types. Antibodies against O1- and O2-serotypes were poorly cross-reactive, and did not efficiently kill or mediate clearance of alternative O-type bacteria, while a distinct and poorly immunogenic O-antigen offered no protection against colonization.

Conclusions

This study suggests that O-antigen variation was introduced to the classical bordetellae via HGT through recombination. Additionally, genetic variation may be maintained within the O-antigen locus because it can provide escape from immunity to different O-antigen types, potentially allowing for the circulation of different Bordetella strains within the same host population.

Enhanced de novo assembly of high throughput pyrosequencing data using whole genome mapping

Despite major advances in next-generation sequencing, assembly of sequencing data, especially data from novel microorganisms or re-emerging pathogens, remains constrained by the lack of suitable reference sequences. De novo assembly is the best approach to achieve an accurate finished sequence, but multiple sequencing platforms or paired-end libraries are often required to achieve full genome coverage. In this study, we demonstrated a method to assemble complete bacterial genome sequences by integrating shotgun Roche 454 pyrosequencing with optical whole genome mapping (WGM). The whole genome restriction map (WGRM) was used as the reference to scaffold de novo assembled sequence contigs through a stepwise process. Large de novo contigs were placed in the correct order and orientation through alignment to the WGRM. De novo contigs that were not aligned to WGRM were merged into scaffolds using contig branching structure information. These extended scaffolds were then aligned to the WGRM to identify the overlaps to be eliminated and the gaps and mismatches to be resolved with unused contigs. The process was repeated until a sequence with full coverage and alignment with the whole genome map was achieved. Using this method we were able to achieved 100% WGRM coverage without a paired-end library. We assembled complete sequences for three distinct genetic components of a clinical isolate of Providencia stuartii: a bacterial chromosome, a novel bla NDM-1 plasmid, and a novel bacteriophage, without separately purifying them to homogeneity.

Minor modifications to the phosphate groups and the C3' acyl chain length of lipid A in two Bordetella pertussis strains, BP338 and 18-323, independently affect Toll-like receptor 4 protein activation

Lipopolysaccharides (LPS) of Bordetella pertussis are important modulators of the immune system. Interaction of the lipid A region of LPS with the Toll-like receptor 4 (TLR4) complex causes dimerization of TLR4 and activation of downstream nuclear factor κB (NFκB), which can lead to inflammation. We have previously shown that two strains of B. pertussis, BP338 (a Tohama I-derivative) and 18-323, display two differences in lipid A structure. 1) BP338 can modify the 1- and 4'-phosphates by the addition of glucosamine (GlcN), whereas 18-323 cannot, and 2) the C3' acyl chain in BP338 is 14 carbons long, but only 10 or 12 carbons long in 18-323. In addition, BP338 lipid A can activate TLR4 to a greater extent than 18-323 lipid A. Here we set out to determine the genetic reasons for the differences in these lipid A structures and the contribution of each structural difference to the ability of lipid A to activate TLR4. We show that three genes of the lipid A GlcN modification (Lgm) locus, lgmA, lgmB, and lgmC (previously locus tags BP0399-BP0397), are required for GlcN modification and a single amino acid difference in LpxA is responsible for the difference in C3' acyl chain length. Furthermore, by introducing lipid A-modifying genes into 18-323 to generate isogenic strains with varying penta-acyl lipid A structures, we determined that both modifications increase TLR4 activation, although the GlcN modification plays a dominant role. These results shed light on how TLR4 may interact with penta-acyl lipid A species.

Differential expression of alpha 4 integrins on effector memory T helper cells during Bordetella infections. Delayed responses in Bordetella pertussis

Bordetella pertussis (B. pertussis) is the causative agent of whooping cough, a respiratory disease that is reemerging worldwide. Mechanisms of selective lymphocyte trafficking to the airways are likely to be critical in the immune response to this pathogen. We compared murine infection by B. pertussis, B. parapertussis, and a pertussis toxin-deleted B. pertussis mutant (BpΔPTX) to test the hypothesis that effector memory T-helper cells (emTh) display an altered pattern of trafficking receptor expression in B. pertussis infection due to a defect in imprinting. Increased cell recruitment to the lungs at 5 days post infection (p.i.) with B. parapertussis, and to a lesser extent with BpΔPTX, coincided with an increased frequency of circulating emTh cells expressing the mucosal-associated trafficking receptors α4β7 and α4β1 while a reduced population of these cells was observed in B. pertussis infection. These cells were highly evident in the blood and lungs in B. pertussis infection only at 25 days p.i. when B. parapertussis and BpΔPTX infections were resolved. Although at 5 days p.i., an equally high percentage of lung dendritic cells (DCs) from all infections expressed maturation markers, this expression persisted only in B. pertussis infection at 25 days p.i. Furthermore, at 5 days p.i with B. pertussis, lung DCs migration to draining lymph nodes may be compromised as evidenced by decreased frequency of CCR7(+) DCs, inhibited CCR7-mediated in vitro migration, and fewer DCs in lung draining lymph nodes. Lastly, a reduced frequency of allogeneic CD4(+) cells expressing α4β1 was detected following co-culture with lung DCs from B. pertussis-infected mice, suggesting a defect in DC imprinting in comparison to the other infection groups. The findings in this study suggest that B. pertussis may interfere with imprinting of lung-associated trafficking receptors on T lymphocytes leading to extended survival in the host and a prolonged course of disease.