Genomic plasticity in natural populations of Bordetella pertussis

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

Genomic Surveillance and Improved Molecular Typing of Bordetella pertussis Using wgMLST

Multilocus sequence typing (MLST) provides allele-based characterization of bacterial pathogens in a standardized framework. However, classical MLST schemes for Bordetella pertussis, the causative agent of whooping cough, seldom reveal diversity among the small number of gene targets and thereby fail to delineate population structure. To improve the discriminatory power of allele-based molecular typing of B. pertussis, we have developed a whole-genome MLST (wgMLST) scheme from 225 reference-quality genome assemblies. Iterative refinement and allele curation resulted in a scheme of 3,506 coding sequences and covering 81.4% of the B. pertussis genome. This wgMLST scheme was further evaluated with data from a convenience sample of 2,389 B. pertussis isolates sequenced on Illumina instruments, including isolates from known outbreaks and epidemics previously characterized by existing molecular assays, as well as replicates collected from individual patients. wgMLST demonstrated concordance with whole-genome single nucleotide polymorphism (SNP) profiles, accurately resolved outbreak and sporadic cases in a retrospective comparison, and clustered replicate isolates collected from individual patients during diagnostic confirmation. Additionally, a reanalysis of isolates from two statewide epidemics using wgMLST reconstructed the population structures of circulating strains with increased resolution, revealing new clusters of related cases. Comparison with an existing core genome (cgMLST) scheme highlights the stable gene content of this bacterium and forms the initial foundation for necessary standardization. These results demonstrate the utility of wgMLST for improving B. pertussis characterization and genomic surveillance during the current pertussis disease resurgence.

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.

Changes in Predominance of Pulsed-Field Gel Electrophoresis Profiles of Bordetella pertussis Isolates, United States, 2000-2012

These changes are concurrrent with other recent molecular changes and may be contributing to US pertussis reemergence.

To clarify the characteristics of circulating Bordetella pertussis isolates, we used pulsed-field gel electrophoresis (PFGE) to analyze 5,262 isolates collected in the United States during 2000–2012. We found 199 PFGE profiles; 5 profiles accounted for 72% of isolates. The most common profile, CDC013, accounted for 35%–46% of isolates tested from 2000–2009; however, the proportion of isolates of this profile rapidly decreased in 2010. Profile CDC237, first seen in 2009, increased rapidly and accounted for 29% of 2012 isolates. No location bias was observed among profiles during 2000–2010, but differences were observed among isolates from different states during 2012. Predominant profiles match those observed in recent European PFGE studies. PFGE profile changes are concurrent with other recent molecular changes in B. pertussis and may be contributing to the reemergence of pertussis in the United States. Continued PFGE monitoring is critical for understanding the changing epidemiology of pertussis.

Genome Structural Diversity among 31 Bordetella pertussis Isolates from Two Recent U.S. Whooping Cough Statewide Epidemics

During 2010 and 2012, California and Vermont, respectively, experienced statewide epidemics of pertussis with differences seen in the demographic affected, case clinical presentation, and molecular epidemiology of the circulating strains. To overcome limitations of the current molecular typing methods for pertussis, we utilized whole-genome sequencing to gain a broader understanding of how current circulating strains are causing large epidemics. Through the use of combined next-generation sequencing technologies, this study compared de novo, single-contig genome assemblies from 31 out of 33 Bordetella pertussis isolates collected during two separate pertussis statewide epidemics and 2 resequenced vaccine strains. Final genome architecture assemblies were verified with whole-genome optical mapping. Sixteen distinct genome rearrangement profiles were observed in epidemic isolate genomes, all of which were distinct from the genome structures of the two resequenced vaccine strains. These rearrangements appear to be mediated by repetitive sequence elements, such as high-copy-number mobile genetic elements and rRNA operons. Additionally, novel and previously identified single nucleotide polymorphisms were detected in 10 virulence-related genes in the epidemic isolates. Whole-genome variation analysis identified state-specific variants, and coding regions bearing nonsynonymous mutations were classified into functional annotated orthologous groups. Comprehensive studies on whole genomes are needed to understand the resurgence of pertussis and develop novel tools to better characterize the molecular epidemiology of evolving B. pertussis populations. IMPORTANCE Pertussis, or whooping cough, is the most poorly controlled vaccine-preventable bacterial disease in the United States, which has experienced a resurgence for more than a decade. Once viewed as a monomorphic pathogen, B. pertussis strains circulating during epidemics exhibit diversity visible on a genome structural level, previously undetectable by traditional sequence analysis using short-read technologies. For the first time, we combine short- and long-read sequencing platforms with restriction optical mapping for single-contig, de novo assembly of 31 isolates to investigate two geographically and temporally independent U.S. pertussis epidemics. These complete genomes reshape our understanding of B. pertussis evolution and strengthen molecular epidemiology toward one day understanding the resurgence of pertussis.

Laboratory adaptation of Bordetella pertussis is associated with the loss of type three secretion system functionality

Although Bordetella pertussis contains and transcribes loci encoding type III secretion system (TTSS) homologues, expression of TTSS-associated proteins has been reported only for non-laboratory-adapted Irish clinical isolates. Here we confirm such a result for clinical isolates obtained from patients treated in Argentinean hospitals. Moreover, we demonstrate that the expression of TTSS-associated proteins is independent both of the year in which the isolate was obtained and of the types of polymorphic alleles for other virulence factors but is dependent on environmental growth conditions. Interestingly, we observed that TTSS-associated protein expression is lost after successive in vitro passages but becomes operative again when bacteria come into contact with the host. This in vivo activation of TTSS expression was observed not only for clinical isolates previously adapted to the laboratory after successive in vitro passages but also for vaccine strains that did not express the system in vitro. The reversibility of TTSS expression, demonstrated by its switching off-on when the bacterium comes into contact with the host, appears to be an adaptive response of this pathogen.

Complete genome sequence of Bordetella pertussis CS, a Chinese pertussis vaccine strain

Bordetella pertussis is the causative agent of pertussis. Here, we report the genome sequence of Bordetella pertussis strain CS, isolated from an infant patient in Beijing and widely used as a vaccine strain for production of an acellular pertussis vaccine in China.

Architecture of Burkholderia cepacia complex sigma70 gene family: evidence of alternative primary and clade-specific factors, and genomic instability

BACKGROUND

The Burkholderia cepacia complex (Bcc) groups bacterial species with beneficial properties that can improve crop yields or remediate polluted sites but can also lead to dramatic human clinical outcomes among cystic fibrosis (CF) or immuno-compromised individuals. Genome-wide regulatory processes of gene expression could explain parts of this bacterial duality. Transcriptional sigma70 factors are components of these processes. They allow the reversible binding of the DNA-dependent RNA polymerase to form the holoenzyme that will lead to mRNA synthesis from a DNA promoter region. Bcc genome-wide analyses were performed to investigate the major evolutionary trends taking place in the sigma70 family of these bacteria.

RESULTS

Twenty sigma70 paralogous genes were detected in the Burkholderia cenocepacia strain J2315 (Bcen-J2315) genome, of which 14 were of the ECF (extracytoplasmic function) group. Non-ECF paralogs were related to primary (rpoD), alternative primary, stationary phase (rpoS), flagellin biosynthesis (fliA), and heat shock (rpoH) factors. The number of sigma70 genetic determinants among this genome was of 2,86 per Mb. This number is lower than the one of Pseudomonas aeruginosa, a species found in similar habitats including CF lungs. These two bacterial groups showed strikingly different sigma70 family architectures, with only three ECF paralogs in common (fecI-like, pvdS and algU). Bcen-J2315 sigma70 paralogs showed clade-specific distributions. Some paralogs appeared limited to the ET12 epidemic clone (ecfA2), particular Bcc species (sigI), the Burkholderia genus (ecfJ, ecfF, and sigJ), certain proteobacterial groups (ecfA1, ecfC, ecfD, ecfE, ecfG, ecfL, ecfM and rpoS), or were broadly distributed in the eubacteria (ecfI, ecfK, ecfH, ecfB, and rpoD-, rpoH-, fliA-like genes). Genomic instability of this gene family was driven by chromosomal inversion (ecfA2), recent duplication events (ecfA and RpoD), localized (ecfG) and large scale deletions (sigI, sigJ, ecfC, ecfH, and ecfK), and a phage integration event (ecfE).

CONCLUSION

The Bcc sigma70 gene family was found to be under strong selective pressures that could lead to acquisition/deletion, and duplication events modifying its architecture. Comparative analysis of Bcc and Pseudomonas aeruginosa sigma70 gene families revealed distinct evolutionary strategies, with the Bcc having selected several alternative primary factors, something not recorded among P. aeruginosa and only previously reported to occur among the actinobacteria.

Significant gene order and expression differences in Bordetella pertussis despite limited gene content variation

Bordetella pertussis, an obligate human pathogen and the agent of whooping cough, is a clonal species, despite the dynamic selection pressures imposed by host immunity and vaccine usage. Because the generation of variation is critical for species evolution, we employed a variety of approaches to examine features of B. pertussis genetic variation. We found a high level of conservation of gene content among 137 B. pertussis strains with different geographical, temporal, and epidemiological associations, using comparative genomic hybridization. The limited number of regions of difference were frequently located adjacent to copies of the insertion element IS481, which is present in high numbers in the B. pertussis chromosome. This repeated sequence appears to provide targets for homologous recombination, resulting in deletion of intervening sequences. Using subtractive hybridization, we searched for previously undetected genes in diverse clinical isolates but did not detect any new genes, indicating that gene acquisition is rare in B. pertussis. In contrast, we found evidence of altered gene order in the several strains that were examined and again found an association of IS481 with sites of rearrangement. Finally, we compared whole-genome expression profiles of different strains and found significant changes in transcript abundance, even in the same strain after as few as 12 laboratory passages. This combination of approaches provides a detailed picture of a pathogenic species with little gene loss or gain but with the capacity to generate variation by rearranging its chromosome and altering gene expression. These findings have broad implications for host adaptation by microbial pathogens.

Characterization of adenylate cyclase-hemolysin gene duplication in a Bordetella pertussis isolate

We describe a clinical isolate of Bordetella pertussis, the agent responsible for whooping cough, composed of at least two clones harboring one or two copies of the cya locus encoding one of the major toxins, adenylate cyclase-hemolysin. No difference was observed between the two clones in murine and cellular models, probably due to the high instability of the cya locus duplication.

Reference system for characterization of Bordetella pertussis pulsed-field gel electrophoresis profiles

Pulsed-field gel electrophoresis (PFGE) has been used as an epidemiological tool for surveillance studies of Bordetella pertussis since the early 1990s. To date there is no standardized procedure for comparison of results, and therefore it has been difficult to directly compare PFGE results between laboratories. We propose a profile-based reference system for PFGE characterization of B. pertussis strain variation and to establish traceability of B. pertussis PFGE results. We initially suggest 35 Swedish reference strains as reference material for PFGE traceability. This reference material is deposited at the Culture Collection of the University of Gothenburg, Gothenburg, Sweden. Altogether, 1,810 Swedish clinical isolates from between 1970 and 2003 were studied, together with the Swedish Pw vaccine strain, six reference strains, and two U.S. isolates. Our system provides evidence that profiles obtained by using only one enzyme, i.e., XbaI, give enough data to analyze the epidemiological relationship between them. Characterization with one enzyme is far less labor intensive, yielding results in half the time than when a two-enzyme procedure is used. Also, we can see that there is a correlation between PFGE profile and pertactin type. One common PFGE profile, BpSR11 (n = 455), showed 100% prn2 and 100% Fim3 when analyzed for pertactin type and serotype. On the other hand, strains with the same profile may express various serotypes when isolated over longer periods of time. Subculturing of the same isolate eight times or lyophilization caused no change in PFGE profile.

Chromosomal constraints in Gram-positive bacteria revealed by artificial inversions

We used artificial chromosome inversions to investigate the chromosomal constraints that preserve genome organization in the Gram-positive bacterium Lactococcus lactis. Large inversions, 80-1260 kb in length, disturbing the symmetry of the origin and terminus of the replication axis to various extents, were constructed using the site-specific Cre-loxP recombination system. These inversions were all mechanistically feasible and fell into various classes according to stability and effect on cell fitness. The L. lactis chromosome supports only to some extent unbalance in length of its replication arms. The location of detrimental inversions allowed identification of two constrained chromosomal regions: a large domain covering one fifth of the genome that encompasses the origin of replication (Ori domain), and a smaller domain located at the opposite of the chromosome (Ter domain).

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

Pathogens of the bacterial genus Bordetella cause respiratory disease in humans and animals. Although virulence and host specificity vary across the genus, the genetic determinants of this diversity remain unidentified. To identify genes that may underlie key phenotypic differences between these species and clarify their evolutionary relationships, we performed a comparative analysis of genome content in 42 Bordetella strains by hybridization of genomic DNA to a microarray representing the genomes of three Bordetella species and by subtractive hybridization. Here we show that B. pertussis and B. parapertussis are predominantly differentiated from B. bronchiseptica by large, species-specific regions of difference, many of which encode or direct synthesis of surface structures, including lipopolysaccharide O antigen, which may be important determinants of host specificity. The species also exhibit sequence diversity at a number of surface protein-encoding loci, including the fimbrial major subunit gene, fim2. Gene loss, rather than gene acquisition, accompanied by the proliferation of transposons, has played a fundamental role in the evolution of the pathogenic bordetellae and may represent a conserved evolutionary mechanism among other groups of microbial pathogens.

Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica

Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica are closely related Gram-negative beta-proteobacteria that colonize the respiratory tracts of mammals. B. pertussis is a strict human pathogen of recent evolutionary origin and is the primary etiologic agent of whooping cough. B. parapertussis can also cause whooping cough, and B. bronchiseptica causes chronic respiratory infections in a wide range of animals. We sequenced the genomes of B. bronchiseptica RB50 (5,338,400 bp; 5,007 predicted genes), B. parapertussis 12822 (4,773,551 bp; 4,404 genes) and B. pertussis Tohama I (4,086,186 bp; 3,816 genes). Our analysis indicates that B. parapertussis and B. pertussis are independent derivatives of B. bronchiseptica-like ancestors. During the evolution of these two host-restricted species there was large-scale gene loss and inactivation; host adaptation seems to be a consequence of loss, not gain, of function, and differences in virulence may be related to loss of regulatory or control functions.

Comparison of bipA alleles within and across Bordetella species

The Bordetella BvgAS signal transduction system controls the expression of at least three phenotypic phases, the Bvg(+) or virulent phase, the Bvg(-) or avirulent phase, and the Bvg(i) or Bvg intermediate phase, which has been hypothesized to be important for transmission. bipA, the first identified Bvg(i)-phase gene, encodes a protein with similarity to the well-characterized bacterial adhesins intimin and invasin. Proteins encoded by the bipA genes present in Bordetella pertussis Tohama I and Bordetella bronchiseptica RB50 differ in the number of 90-amino-acid repeats which they possess and in the sequence of the C-terminal domain. To investigate the possibility that bipA alleles segregate according to host specificity and to gain insight into the role of BipA and the Bvg(i) phase in the Bordetella infectious cycle, we compared bipA alleles across members of the B. bronchiseptica cluster, which includes both human-infective (B. pertussis and B. parapertussis(hu)) and non-human-infective (B. bronchiseptica and B. parapertussis(ov)) strains. bipA genes were present in most, but not all, strains. All bipA genes present in B. bronchiseptica strains were identical to bipA of RB50 (at least with regard to the DNA sequence of the 3' C-terminal-domain-encoding region, the number of 90-amino-acid repeats encoded, and expression patterns). Although all bipA genes present in the other Bordetella strains were identical in the 3' C-terminal-domain-encoding region to bipA of B. pertussis Tohama I, they varied in the number of 90-amino-acid repeats that they encoded and in expression level. Notably, the genes present in B. parapertussis(hu) strains were pseudogenes, and the genes present in B. parapertussis(ov) strains were expressed at significantly reduced levels compared with the levels in B. pertussis and B. bronchiseptica strains. Our results indicate that there is a correlation between specific bipA alleles and specific hosts. They also support the hypothesis that both horizontal gene transfer and fine-tuning of gene expression patterns contribute to the evolution of host adaptation in lineages of the B. bronchiseptica cluster.

Tandem tetramer-based microsatellite fingerprinting for typing of Proteus mirabilis strains

Two microsatellite tandem repeated tetramers, (GACA)(4) and (CAAT)(4), were used for Proteus mirabilis strain differentiation. The microsatellite-based PCR tests were applied for the examination of interstrain diversity for 87 P. mirabilis strains. Forty-six of the investigated strains were clinical isolates (5 were hospital isolates and 39 were outpatient clinic isolates); 42 strains were derived from the Kauffmann-Perch collection of laboratory strains. Fingerprinting done with the tetramers had a high discrimination ability [0.992 and 0.940 for (GACA)(4) and (CAAT)(4), respectively]. The distributions of clinical isolates among well-defined laboratory strains, determined by numerical analysis (unweighted pair-group method with arithmetic averages; Dice similarity coefficient), proved their genetic similarity to reference strains in the Kauffmann-Perch collection. This analysis also indicated that it is possible to estimate some phenotypic properties of P. mirabilis clinical isolates solely on the basis of microsatellite fingerprinting.

Two-dimensional display and whole genome comparison of bacterial pathogen genomes of high G+C DNA content

High-resolution comparison of bacterial genomes facilitates the identification of the genetic changes responsible for clinically relevant phenotypes. For this purpose we have established a method for the display and comparison of high G+C bacterial genomes in two dimensions. Here we describe the application of two-dimensional bacterial genomic display to resolve the genomes of Bordetella pertussis, Mycobacterium avium and Mycobacterium tuberculosis, and its utility in strain comparison and detection of insertion and substitution mutations.

Use of temporal temperature gradient gel electrophoresis to identify flaA and fim3 sequence types in Bordetella bronchiseptica

AIMS

The aim of this study was to develop an approach to detect variation in the flaA and fim3 genes amongst animal isolates of Bordetella bronchiseptica using temporal temperature gradient gel electrophoresis (TTGE).

METHODS AND RESULTS

Amplicons representing three flaA and two fim3 sequence types were subjected to TTGE analysis. It was possible to distinguish clearly between each of the sequence types using TTGE.

CONCLUSION

This suggests that TTGE could be a useful tool for studying the epidemiology of B. bronchiseptica.

SIGNIFICANCE AND IMPACT OF THE STUDY

PCR amplification coupled to TTGE offers a general method for the rapid screening of large numbers of microbial strains for variations in gene sequences.

Evolutionary trends in the genus Bordetella

The genus Bordetella comprises seven species with pathogenic potential for different host organisms. This article attempts to review our current knowledge about the systematics and evolution of this important group of pathogens, their relationship to environmental microorganisms and about molecular mechanisms of host adaptation.

Evaluating genome dynamics: the constraints on rearrangements within bacterial genomes

Inversions and translocations distinguish the genomes of closely related bacterial species, but most of these rearrangements preserve the relationship between the rearranged fragments and the axis of chromosome replication. Within species, such rearrangements are found less frequently, except in the case of clinical isolates of human pathogens, where rearrangements are very frequent.