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

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

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

Recent increase in pertussis incidence in South Korea: An age-period-cohort analysis

OBJECTIVES

Pertussis or whooping cough—one of the most contagious diseases—is caused by the Gram-negative bacterium Bordetella pertussis. Despite a high vaccination rate, Korea recently experienced a resurgence of pertussis. This study explores patterns and possible explanations for this resurgence through an age-period-cohort analysis.

METHODS

Using secondary data from the infectious disease portal of the Korea Disease Control and Prevention Agency and the Korea Statistical Information Service of Statistics Korea, this study analyzed the incidence of pertussis in Korea to determine which factors contributed to the recent increase using an age-period-cohort model.

RESULTS

Analysis of the age effect indicated that the age group most vulnerable to pertussis was 0-year to 2-year-olds. Analysis of the period effect showed a sharp increase in the incidence rate after 2016. Analysis of the cohort effect showed a significant decrease in incidence beginning with the 1955 birth cohort, with the risk increasing again with the 2000s birth cohort.

CONCLUSIONS

Previous studies have suggested 3 main possible explanations for our results. First, the increased incidence rate can be attributed to contact rates. Second, the rate of immunity through natural exposure has decreased due to the low number of circulating pathogens, in turn affecting the trend of infection. Lastly, variations in pathogens may have also contributed to the increase in incidence. Given that the most significant increase in incidence was observed among infants younger than 1 year old, sufficient maternal immunity must be prioritized to provide passive immunity to newborns via the placenta.

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

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

The epidemic of erythromycin-resistant Bordetella pertussis with limited genome variation associated with pertussis resurgence in China

BACKGROUND

The resurgence of Bordetella pertussis infections leading to whooping cough is a concern in many parts of the world. The number of pertussis cases in China has increased significantly since 2013.

RESEARCH DESIGN AND METHODS

In this study, whole-genome sequencing analysis was performed for 388 B. pertussis strains isolated in China from the 1970s to 2018, combining 594 published strains from around the world.

RESULTS

This study revealed that lineage V diverged about 50 years ago in China, while lineage IV is dominant in the other countries. It also revealed that the erythromycin-resistant sub-lineages Va, Vb, and Vc with limited genomic variation emerged 11 ~ 12 years ago. These three sub-lineages were identified after the co-purified acellular vaccines (cp-ACVs) completely replaced the previous whole cell vaccines (WCVs) after the national immunization program of 2012. It suggests that the cp-ACVs cannot induce immunity that is potent enough to restrict the spread of the lineage V, antibiotic abuse further favors the spread of this lineage in China.

CONCLUSIONS

These findings demand a reassessment of the immunization strategy and development of new vaccines in China to stop the resurgence and drug resistance of B. pertussis.

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.

Engineering of an Iranian Bordetella pertussis strain producing inactive pertussis toxin

Introduction. Differences between the genomic and virulence profile of Bordetella pertussis circulating strains and vaccine strains are considered as one of the important reasons for the resurgence of whooping cough (pertussis) in the world. Genetically inactivated B. pertussis is one of the new strategies to generate live-attenuated vaccines against whooping cough.Aim. The aim of this study was to construct a B. pertussis strain based on a predominant profile of circulating Iranian isolates that produces inactivated pertussis toxin (PTX).Methodology. The B. pertussis strain BPIP91 with predominant genomic and virulence pattern was selected from the biobank of the Pasteur Institute of Iran. A BPIP91 derivative with R9K and E129G alterations in the S1 subunit of PTX (S1mBPIP91) was constructed by the site-directed mutagenesis and homologous recombination. Genetic stability and antigen expression of S1mBPIP91 were tested by serially in vitro passages and immunoblot analyses, respectively. The reduction in toxicity of S1mBPIP91 was determined by Chinese hamster ovary (CHO) cell clustering.Results. All constructs and S1mBPIP91 were confirmed via restriction enzyme analysis and DNA sequencing. The engineered mutations in S1mBPIP91 were stable after 20 serial in vitro passages. The production of virulence factors was also confirmed in S1mBPIP91. The CHO cell-clustering test demonstrated the reduction in PTX toxicity in S1mBPIP91.Conclusion. A B. pertussis of the predominant genomic and virulence lineage in Iran was successfully engineered to produce inactive PTX. This attenuated strain will be useful to further studies to develop both whole cell and acellular pertussis vaccines.

How Genomics Is Changing What We Know About the Evolution and Genome of Bordetella pertussis

The evolution of Bordetella pertussis from a common ancestor similar to Bordetella bronchiseptica has occurred through large-scale gene loss, inactivation and rearrangements, largely driven by the spread of insertion sequence element repeats throughout the genome. B. pertussis is widely considered to be monomorphic, and recent evolution of the B. pertussis genome appears to, at least in part, be driven by vaccine-based selection. Given the recent global resurgence of whooping cough despite the wide-spread use of vaccination, a more thorough understanding of B. pertussis genomics could be highly informative. In this chapter we discuss the evolution of B. pertussis, including how vaccination is changing the circulating B. pertussis population at the gene-level, and how new sequencing technologies are revealing previously unknown levels of inter- and intra-strain variation at the genome-level.

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.

Epidemiology of transmissible diseases: Array hybridization and next generation sequencing as universal nucleic acid-mediated typing tools

The magnitude of interest in the epidemiology of transmissible human diseases is reflected in the vast number of tools and methods developed recently with the expressed purpose to characterize and track evolutionary changes that occur in agents of these diseases over time. Within the past decade a new suite of such tools has become available with the emergence of the so-called "omics" technologies. Among these, two are exponents of the ongoing genomic revolution. Firstly, high-density nucleic acid probe arrays have been proposed and developed using various chemical and physical approaches. Via hybridization-mediated detection of entire genes or genetic polymorphisms in such genes and intergenic regions these so called "DNA chips" have been successfully applied for distinguishing very closely related microbial species and strains. Second and even more phenomenal, next generation sequencing (NGS) has facilitated the assessment of the complete nucleotide sequence of entire microbial genomes. This technology currently provides the most detailed level of bacterial genotyping and hence allows for the resolution of microbial spread and short-term evolution in minute detail. We will here review the very recent history of these two technologies, sketch their usefulness in the elucidation of the spread and epidemiology of mostly hospital-acquired infections and discuss future developments.

The History of Bordetella pertussis Genome Evolution Includes Structural Rearrangement

Despite high pertussis vaccine coverage, reported cases of whooping cough (pertussis) have increased over the last decade in the United States and other developed countries. Although Bordetella pertussis is well known for its limited gene sequence variation, recent advances in long-read sequencing technology have begun to reveal genomic structural heterogeneity among otherwise indistinguishable isolates, even within geographically or temporally defined epidemics. We have compared rearrangements among complete genome assemblies from 257 B. pertussis isolates to examine the potential evolution of the chromosomal structure in a pathogen with minimal gene nucleotide sequence diversity. Discrete changes in gene order were identified that differentiated genomes from vaccine reference strains and clinical isolates of various genotypes, frequently along phylogenetic boundaries defined by single nucleotide polymorphisms. The observed rearrangements were primarily large inversions centered on the replication origin or terminus and flanked by IS481, a mobile genetic element with >240 copies per genome and previously suspected to mediate rearrangements and deletions by homologous recombination. These data illustrate that structural genome evolution in B. pertussis is not limited to reduction but also includes rearrangement. Therefore, although genomes of clinical isolates are structurally diverse, specific changes in gene order are conserved, perhaps due to positive selection, providing novel information for investigating disease resurgence and molecular epidemiology.IMPORTANCE Whooping cough, primarily caused by Bordetella pertussis, has resurged in the United States even though the coverage with pertussis-containing vaccines remains high. The rise in reported cases has included increased disease rates among all vaccinated age groups, provoking questions about the pathogen's evolution. The chromosome of B. pertussis includes a large number of repetitive mobile genetic elements that obstruct genome analysis. However, these mobile elements facilitate large rearrangements that alter the order and orientation of essential protein-encoding genes, which otherwise exhibit little nucleotide sequence diversity. By comparing the complete genome assemblies from 257 isolates, we show that specific rearrangements have been conserved throughout recent evolutionary history, perhaps by eliciting changes in gene expression, which may also provide useful information for molecular epidemiology.

Emerging Bordetella pertussis Strains Induce Enhanced Signaling of Human Pattern Recognition Receptors TLR2, NOD2 and Secretion of IL-10 by Dendritic Cells

Vaccines against pertussis have been available for more than 60 years. Nonetheless, this highly contagious disease is reemerging even in countries with high vaccination coverage. Genetic changes of Bordetella pertussis over time have been suggested to contribute to the resurgence of pertussis, as these changes may favor escape from vaccine-induced immunity. Nonetheless, studies on the effects of these bacterial changes on the immune response are limited. Here, we characterize innate immune recognition and activation by a collection of genetically diverse B. pertussis strains isolated from Dutch pertussis patients before and after the introduction of the pertussis vaccines. For this purpose, we used HEK-Blue cells transfected with human pattern recognition receptors TLR2, TLR4, NOD2 and NOD1 as a high throughput system for screening innate immune recognition of more than 90 bacterial strains. Physiologically relevant human monocyte derived dendritic cells (moDC), purified from peripheral blood of healthy donors were also used. Findings indicate that, in addition to inducing TLR2 and TLR4 signaling, all B. pertussis strains activate the NOD-like receptor NOD2 but not NOD1. Furthermore, we observed a significant increase in TLR2 and NOD2, but not TLR4, activation by strains circulating after the introduction of pertussis vaccines. When using moDC, we observed that the recently circulating strains induced increased activation of these cells with a dominant IL-10 production. In addition, we observed an increased expression of surface markers including the regulatory molecule PD-L1. Expression of PD-L1 was decreased upon blocking TLR2. These in vitro findings suggest that emerging B. pertussis strains have evolved to dampen the vaccine-induced inflammatory response, which would benefit survival and transmission of this pathogen. Understanding how this disease has resurged in a highly vaccinated population is crucial for the design of improved vaccines against pertussis.

The Pertussis resurgence: putting together the pieces of the puzzle

Pertussis incidence is rising in almost every country where acellular pertussis (aP) vaccines have been introduced, and is occurring across all age groups from infancy to adulthood. The key question is why? While several known factors such as waning of immunity, detection bias due to more sensitive tests and higher awareness of the disease among practitioners, and evolutionary shifts among B. pertussis all likely contribute, collectively, these do not adequately explain the existing epidemiologic data, suggesting that additional factors also contribute. Key amongst these is recent data indicating that the immune responses induced by aP vaccines differ fundamentally from those induced by the whole cell pertussis (wP) vaccines, and do not lead to mucosal immunity. If so, it appears likely that differences in how the two categories of vaccines work, may be pivotal to our overall understanding of the pertussis resurgence.

Better colonisation of newly emerged Bordetella pertussis in the co-infection mouse model study

Molecular epidemiological data indicates that the resurgence of pertussis (whooping cough) in populations with high vaccine coverage is associated with genomic adaptation of Bordetella pertussis, the causative agent of the disease, to vaccine selection pressure. We have previously shown that in the period after the introduction of acellular pertussis vaccine (ACV), the majority of circulating strains in Australia switched to single nucleotide polymorphism (SNP) cluster I (carrying ptxP3/prn2), replacing SNP cluster II (carrying ptxP1/prn3). In this study, we carried out an in vivo competition assay using a mouse model infected with SNP cluster I and II B. pertussis strains from Australia. We found that the SNP cluster I strain colonised better than the SNP cluster II strain, in both naïve and immunised mice, suggesting that SNP cluster I strains had better fitness regardless of immunisation status of the host, consistent with SNP cluster I strains replacing SNP cluster II. Nevertheless, we found that ACV enhanced clearance of both SNP cluster I and II strains from the mouse respiratory tract.

Bordetella pertussis epidemiology and evolution in the light of pertussis resurgence

Whooping cough, or pertussis, is resurgent in many countries world-wide. This is linked to switching from the use of whole cell vaccines to acellular vaccines in developed countries. Current evidence suggests that this has resulted in the earlier waning of vaccine-induced immunity, an increase in asymptomatic infection with concomitant increases in transmission and increased selection pressure for Bordetellapertussis variants that are better able to evade vaccine-mediated immunity than older isolates. This review discusses recent findings in B. pertussis epidemiology and evolution in the light of pertussis resurgence, and highlights the important role for genomics-based studies in monitoring B. pertussis adaptation.

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.

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.

A Functional Tricarboxylic Acid Cycle Operates during Growth of Bordetella pertussis on Amino Acid Mixtures as Sole Carbon Substrates

It has been claimed that citrate synthase, aconitase and isocitrate dehydrogenase activities are non-functional in Bordetella pertussis and that this might explain why this bacterium’s growth is sometimes associated with accumulation of polyhydroxybutyrate (PHB) and/or free fatty acids. However, the sequenced genome includes the entire citric acid pathway genes. Furthermore, these genes were expressed and the corresponding enzyme activities detected at high levels for the pathway when grown on a defined medium imitating the amino acid content of complex media often used for growth of this pathogenic microorganism. In addition, no significant PHB or fatty acids could be detected. Analysis of the carbon balance and stoichiometric flux analysis based on specific rates of amino acid consumption, and estimated biomass requirements coherent with the observed growth rate, clearly indicate that a fully functional tricarboxylic acid cycle operates in contrast to previous reports.

Protecting Newborns Against Pertussis: Treatment and Prevention Strategies

Pertussis is a potentially severe respiratory disease, which affects all age groups from young infants to older adults and is responsible for an estimated 195,000 deaths occurred globally in 2008. Active research is ongoing to better understand the pathogenesis, immunology, and diagnosis of pertussis. For diagnosis, molecular assays (e.g., polymerase chain reaction) for detection of Bordetella pertussis have become more widely available and support improved outbreak detection. In children, pertussis vaccines have been incorporated into routine immunization schedules and deployed for pertussis outbreak control. Lower levels of vaccine coverage are now being observed in communities where vaccine hesitancy is rising. Additionally, recognition that newborn babies are at risk of pertussis in the USA and UK has led to recommendations to immunize pregnant women. Among adolescents and older adults in the USA, Tetanus Toxoid, Reduced Diphtheria Toxoid and Acellular pertussis (Tdap) Vaccines are recommended, but substantial individual- and system-level barriers exist that will make achieving national Healthy People 2020 targets for immunization challenging. Current antimicrobial regimens for pertussis are focused on reducing the severity of disease, reducing rates of sequelae, and minimizing transmission of infection to susceptible individuals. Continued surveillance for pertussis will be important to identify opportunities for reducing young infants' exposure and reducing the impact of outbreaks among school-aged children. Laboratory-based surveillance for newly emerging strains of B. pertussis will be important to identify strains that may evade protection elicited by currently available vaccines. Efforts to develop new-generation pertussis vaccines should be considered now in anticipation of vaccine development programs, which may require ten or more years to deliver a licensed vaccine.

Bordetella pertussis evolution in the (functional) genomics era

The incidence of whooping cough caused by Bordetella pertussis in many developed countries has risen dramatically in recent years. This has been linked to the use of an acellular pertussis vaccine. In addition, it is thought that B. pertussis is adapting under acellular vaccine mediated immune selection pressure, towards vaccine escape. Genomics-based approaches have revolutionized the ability to resolve the fine structure of the global B. pertussis population and its evolution during the era of vaccination. Here, we discuss the current picture of B. pertussis evolution and diversity in the light of the current resurgence, highlight import questions raised by recent studies in this area and discuss the role that functional genomics can play in addressing current knowledge gaps.

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.

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.

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.

Recent trends of antigenic variation in Bordetella pertussis isolates in Korea

Pertussis is a representative vaccine-preventable disease. However, there have been recent outbreaks in countries where even higher vaccination against the disease. One reason is the emergence of antigenic variants, which are different to vaccine type. In Korea, reported cases have rapidly increased since 2009. Therefore, we analyzed genotype of strains isolated in 2011-2012 by multilocus sequence typing method. As expected, the genotype profiles of tested genes dramatically changed. The major sequence type changed from ST1 to ST2, and new sequence type (ST8) appeared. In the minimum spanning tree, recent isolates belonging to the ACC-I-ST3 subgroup were detected that were composed of ST2, ST3, and ST6. In particular, the ST2 frequency increased to 81%. The novel ST8 was linked to the increased frequency of ST2. In addition, toxic strains carrying the ptxP3 promoter type were confirmed. This ptxP3 type emerged from 2009 and its frequency had increased to 100% in 2012. Based on these results, it can be inferred that the genotypic changes in the currently circulating strains are strongly associated with the recent increasing of pertussis in Korea. Therefore, the surveillance system should be strengthened, and genetic characterization of the isolates should be expanded to the whole genome sequence level.

Age related differences in dynamics of specific memory B cell populations after clinical pertussis infection

For a better understanding of the maintenance of immune mechanisms to Bordetella pertussis (Bp) in relation to age, we investigated the dynamic range of specific B cell responses in various age-groups at different time points after a laboratory confirmed pertussis infection. Blood samples were obtained in a Dutch cross sectional observational study from symptomatic pertussis cases. Lymphocyte subpopulations were phenotyped by flowcytometry before and after culture. Memory B (B_mem) cells were differentiated into IgG antibody secreting cells (ASC) by polyclonal stimulation and detected by an ELISPOT assay specific for pertussis antigens pertussis toxin (Ptx), filamentous haemagglutinin (FHA) and pertactin (Prn). Bp antigen specific IgG concentrations in plasma were determined using multiplex technology. The majority of subjects having experienced a clinical pertussis episode demonstrated high levels of both Bp specific IgG and B_mem cell levels within the first 6 weeks after diagnosis. Significantly lower levels were observed thereafter. Waning of cellular and humoral immunity to maintenance levels occurred within 9 months after antigen encounter. Age was found to determine the maximum but not base-line frequencies of B_mem cell populations; higher levels of B_mem cells specific for Ptx and FHA were reached in adults and (pre-) elderly compared to under-fours and schoolchildren in the first 6 weeks after Bp exposure, whereas not in later phases. This age effect was less obvious for specific IgG levels. Nonetheless, subjects' levels of specific B_mem cells and specific IgG were weakly correlated. This is the first study to show that both age and closeness to last Bp encounter impacts the size of Bp specific B_mem cell and plasma IgG levels.

Genome-wide gene expression analysis of Bordetella pertussis isolates associated with a resurgence in pertussis: elucidation of factors involved in the increased fitness of epidemic strains

Bordetella pertussis (B. pertussis) is the causative agent of whooping cough, which is a highly contagious disease in the human respiratory tract. Despite vaccination since the 1950s, pertussis remains the most prevalent vaccine-preventable disease in developed countries. A recent resurgence pertussis is associated with the expansion of B. pertussis strains with a novel allele for the pertussis toxin (ptx) promoter ptxP3 in place of resident ptxP1 strains. The recent expansion of ptxP3 strains suggests that these strains carry mutations that have increased their fitness. Compared to the ptxP1 strains, ptxP3 strains produce more Ptx, which results in increased virulence and immune suppression. In this study, we investigated the contribution of gene expression changes of various genes on the increased fitness of the ptxP3 strains. Using genome-wide gene expression profiling, we show that several virulence genes had higher expression levels in the ptxP3 strains compared to the ptxP1 strains. We provide the first evidence that wildtype ptxP3 strains are better colonizers in an intranasal mouse infection model. This study shows that the ptxP3 mutation and the genetic background of ptxP3 strains affect fitness by contributing to the ability to colonize in a mouse infection model. These results show that the genetic background of ptxP3 strains with a higher expression of virulence genes contribute to increased fitness.

Re-emergence of pertussis: what are the solutions?

Whooping cough, due to Bordetella pertussis and Bordetella parapertussis, is an important cause of childhood morbidity and mortality. Despite widespread pertussis immunization in childhood, there are an estimated 50 million cases and 300,000 deaths due to pertussis globally each year. Infants who are too young to be vaccinated, children who are partially vaccinated and fully-vaccinated persons with waning immunity are especially vulnerable to disease. Since pertussis is one of the vaccine-preventable diseases on the rise, additional vaccine approaches are needed. These approaches include vaccination of newborns, additional booster doses for older adolescents and adults, and immunization of pregnant women with existing vaccines. Innovative new vaccines are also being studied. Each of these options will be discussed and their potential impact on pertussis control assessed.

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.

The evolutionary consequences of blood-stage vaccination on the rodent malaria Plasmodium chabaudi

A candidate malaria vaccine promoted the evolution of more virulent malaria parasites in mice.

Malaria vaccine developers are concerned that antigenic escape will erode vaccine efficacy. Evolutionary theorists have raised the possibility that some types of vaccine could also create conditions favoring the evolution of more virulent pathogens. Such evolution would put unvaccinated people at greater risk of severe disease. Here we test the impact of vaccination with a single highly purified antigen on the malaria parasite Plasmodium chabaudi evolving in laboratory mice. The antigen we used, AMA-1, is a component of several candidate malaria vaccines currently in various stages of trials in humans. We first found that a more virulent clone was less readily controlled by AMA-1-induced immunity than its less virulent progenitor. Replicated parasites were then serially passaged through control or AMA-1 vaccinated mice and evaluated after 10 and 21 rounds of selection. We found no evidence of evolution at the ama-1 locus. Instead, virulence evolved; AMA-1-selected parasites induced greater anemia in naïve mice than both control and ancestral parasites. Our data suggest that recombinant blood stage malaria vaccines can drive the evolution of more virulent malaria parasites.

Vaccination can drive the evolution of pathogens. Most obviously, molecules targeted by vaccine-induced immunity can change. Such evolution makes vaccines less effective. A different possibility is that more virulent pathogens are favored in vaccinated hosts. In that case, vaccination would create pathogens that cause more harm to unvaccinated individuals. To test this idea, we studied a rodent malaria parasite in laboratory mice immunized with a component of malaria vaccines currently in human trials. We found that a more virulent parasite clone was less well controlled by vaccine-induced immunity than was its less virulent ancestor. We then passaged parasites through sham- or vaccinated mice to study how the parasites might evolve after multiple rounds of infection of mouse hosts. The parasite molecule targeted by the vaccine did not change during this process. Instead, the parasites became more virulent if they evolved in vaccinated hosts. Our data suggest that some vaccines can drive the evolution of more virulent parasites.

High confidence prediction of essential genes in Burkholderia cenocepacia

Background

Essential genes are absolutely required for the survival of an organism. The identification of essential genes, besides being one of the most fundamental questions in biology, is also of interest for the emerging science of synthetic biology and for the development of novel antimicrobials. New antimicrobial therapies are desperately needed to treat multidrug-resistant pathogens, such as members of the Burkholderia cepacia complex.

Methodology/Principal Findings

We hypothesize that essential genes may be highly conserved within a group of evolutionary closely related organisms. Using a bioinformatics approach we determined that the core genome of the order Burkholderiales consists of 649 genes. All but two of these identified genes were located on chromosome 1 of Burkholderia cenocepacia. Although many of the 649 core genes of Burkholderiales have been shown to be essential in other bacteria, we were also able to identify a number of novel essential genes present mainly, or exclusively, within this order. The essentiality of some of the core genes, including the known essential genes infB, gyrB, ubiB, and valS, as well as the so far uncharacterized genes BCAL1882, BCAL2769, BCAL3142 and BCAL3369 has been confirmed experimentally in B. cenocepacia.

Conclusions/Significance

We report on the identification of essential genes using a novel bioinformatics strategy and provide bioinformatics and experimental evidence that the large majority of the identified genes are indeed essential. The essential genes identified here may represent valuable targets for the development of novel antimicrobials and their detailed study may shed new light on the functions required to support life.

BpsR modulates Bordetella biofilm formation by negatively regulating the expression of the Bps polysaccharide

Bordetella bacteria are Gram-negative respiratory pathogens of animals, birds, and humans. A hallmark feature of some Bordetella species is their ability to efficiently survive in the respiratory tract even after vaccination. Bordetella bronchiseptica and Bordetella pertussis form biofilms on abiotic surfaces and in the mouse respiratory tract. The Bps exopolysaccharide is one of the critical determinants for biofilm formation and the survival of Bordetella in the murine respiratory tract. In order to gain a better understanding of regulation of biofilm formation, we sought to study the mechanism by which Bps expression is controlled in Bordetella. Expression of bpsABCD (bpsA-D) is elevated in biofilms compared with levels in planktonically grown cells. We found that bpsA-D is expressed independently of BvgAS. Subsequently, we identified an open reading frame (ORF), BB1771 (designated here bpsR), that is located upstream of and in the opposite orientation to the bpsA-D locus. BpsR is homologous to the MarR family of transcriptional regulators. Measurement of bpsA and bpsD transcripts and the Bps polysaccharide levels from the wild-type and the ΔbpsR strains suggested that BpsR functions as a repressor. Consistent with enhanced production of Bps, the bpsR mutant displayed considerably more structured biofilms. We mapped the bpsA-D promoter region and showed that purified BpsR protein specifically bound to the bpsA-D promoter. Our results provide mechanistic insights into the regulatory strategy employed by Bordetella for control of the production of the Bps polysaccharide and biofilm formation.

Differences in the genomic content of Bordetella pertussis isolates before and after introduction of pertussis vaccines in four European countries

Resurgence of pertussis has been observed in many countries with high vaccination coverage and clonal expansion of certain Bordetella pertussis strains has been associated with recent epidemics in Europe. It is known that vaccinations have selected strains which are different from those used for vaccine production. However, little is known about the differences in genomic content of strains circulating before the vaccination was introduced. In this study, we compared the genomes of 39 vaccine strains and old clinical isolates (isolated 1941-1984) collected from Finland (n = 5), Poland (n = 14), Serbia (n = 10) and the UK (n = 10). The analysis included genotyping, pulsed-field gel electrophoresis (PFGE) and comparative genomic hybridisation (CGH). Compared to the strain Tohama I, the European isolates analyzed have lost three major regions of difference (RD3, 5 and 29). However, difference in frequency of the absent RDs 3 (BP0910A-BP0934), 5 (BP1135-BP1141) or 29 (BP1225) was observed among isolates from the four countries. Of the isolates with absent RD5, half had also a duplicated region in the genome. All four RDs (RD22 (BB0535-BB0541), 23 (BB0916-BB0921), 24 (BB1140-BB1158) and 26 (BB4880-BB4888)) absent in Tohama I were present in majority of the tested isolates. Results obtained from PFGE analysis correlated well with those of CGH. Recently a novel pertussis toxin promoter allele (ptxP3) was described. Isolates with ptxP3 have replaced resident ptxP1 isolates in the countries where this was investigated. When the recent isolates, collected in 2000-2004, selected from the four countries were examined, the ptxP3 allele was found in all countries except Poland. Our result indicates that at least three clusters of B. pertussis circulated in Europe in pre- and early vaccine era and their genomes were distinct from that of the reference strain Tohama I. Although progressive gene loss occurs in B. pertussis population with time, difference in frequency of the lost genes were observed among isolates from the four countries. The observed differences in genomic content might be vaccine-driven.

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.