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

Deciphering Bordetella pertussis epidemiology through culture-independent multiplex amplicon and metagenomic sequencing

Whooping cough (pertussis) has re-emerged despite high vaccine coverage in Australia and many other countries worldwide, partly attributable to genetic adaptation of the causative organism, Bordetella pertussis, to vaccines. Therefore, genomic surveillance has become essential to monitor circulating strains for these genetic changes. However, increasing uptake of PCR for the diagnosis of pertussis has affected the availability of cultured isolates for typing. In this study, we evaluated the use of targeted multiplex PCR (mPCR) amplicon sequencing and shotgun metagenomic sequencing for culture-independent typing of B. pertussis directly from respiratory swabs. We developed a nine-target mPCR amplicon assay that could accurately type major lineages [ptxP3/non-ptxpP3, fim3A/B, fhaB3/non-fhaB3, and epidemic lineages (ELs) 1-5] circulating in Australia. Validation using DNA from isolates and 178 residual specimens collected in 2010-2012 (n = 87) and 2019 (n = 91) showed that mPCR amplicon sequencing was highly sensitive with a limit of detection of 4.6 copies [IS481 cycle threshold (Ct) 27.3]. Shotgun metagenomic sequencing was successful in genotyping B. pertussis in 84% of clinical specimens with PCR Ct < 24 and was concordant with mPCR typing results. The results revealed an expansion of EL4 strains from 2010 to 2012 to 2019 in Australia and identified unrecognized co-circulating cases of Bordetella holmesii. This study provides valuable insight into the circulating lineages in Australia prior to the COVID-19 pandemic during which border closure and other interventions reduced pertussis cases to an all-time low, and paves the way for the genomic surveillance of B. pertussis in the era of culture-independent PCR-based diagnosis.

IMPORTANCE

In this paper, we evaluated the use of targeted multiplex PCR (mPCR) amplicon sequencing and shotgun metagenomic sequencing for culture-independent typing of Bordetella pertussis directly in respiratory swabs. We first developed a novel targeted mPCR amplicon sequencing assay that can type major circulating lineages and validated its accuracy and sensitivity on 178 DNA extracts from clinical swabs. We also demonstrate the feasibility of using deep metagenomic sequencing for determining strain lineage and markers of virulence, vaccine adaptation, macrolide resistance, and co-infections. Our culture-independent typing methods applied to clinical specimens revealed the expansion of a major global epidemic lineage in Australia (termed EL4) just prior to the COVID-19 pandemic. It also detected cases of previously hidden co-infections from another Bordetella species called Bordetella holmesii. These findings offer valuable insight into the circulating pertussis lineages in Australia prior to the COVID-19 pandemic during which border closure and other interventions reduced pertussis cases to an all-time low. It also provides comparative data for future surveillance as pertussis resurgence after the COVID-19 pandemic has been reported this year in Australia and many other countries. Overall, our paper demonstrates the utility, sensitivity, and specificity of mPCR amplicon and metagenomic sequencing-based culture-independent typing of B. pertussis, which not only paves the way for culture-independent genomic surveillance of B. pertussis but also for other pathogens in the era of PCR-based diagnosis.

Pertussis vaccines, epidemiology and evolution

Pertussis, which is caused by Bordetella pertussis, has plagued humans for at least 800 years, is highly infectious and can be fatal in the unvaccinated, especially very young infants. Although the rollout of whole-cell pertussis (wP) vaccines in the 1940s and 1950s was associated with a drastic drop in incidence, concerns regarding the reactogenicity of wP vaccines led to the development of a new generation of safer, acellular (aP) vaccines that have been adopted mainly in high-income countries. Over the past 20 years, some countries that boast high aP coverage have experienced a resurgence in pertussis, which has led to substantial debate over the basic immunology, epidemiology and evolutionary biology of the bacterium. Controversy surrounds the duration of natural immunity and vaccine-derived immunity, the ability of vaccines to prevent transmission and severe disease, and the impact of evolution on evading vaccine immunity. Resolving these issues is made challenging by incomplete detection of pertussis cases, the absence of a serological marker of immunity, modest sequencing of the bacterial genome and heterogeneity in diagnostic methods of surveillance. In this Review, we lay out the complexities of contemporary pertussis and, where possible, propose a parsimonious explanation for apparently incongruous observations.

Insertion sequence elements and unique symmetrical genomic regions mediate chromosomal inversions in Streptococcus pyogenes

Chromosomal inversions are a phenomenon in many bacterial species, often across the axis of replication. Inversions have been shown to alter gene expression, changing persistence of colonisation and infection following environmental stresses. In Streptococcus pyogenes, inversions have been reported. However, frequency and molecular markers of inversions have not been systematically examined. Here, 249 complete S.pyogenes genomes were analysed using a pangenomic core gene synteny framework to identify sequences associated with inversions. 47% of genomes (118/249) contained at least one inversion, from 23 unique inversion locations. Chromosomal locations enabling inversions were usually associated with mobile elements (insertion sequences n = 9 and prophages n = 7). Two insertion sequences, IS1548 and IS1239, accounted for >80% of insertion sequences and were the only insertion sequences associated with inversions. The most observed inversion location (n = 104 genomes, 88% of genomes with an inversion) occurs between two conserved regions encoding rRNAs, tRNAs and sigma factor genes. The regions are symmetrically placed around the origin of replication forming a unique chromosomal structure in S. pyogenes, relative to other streptococci. Cataloging of the chromosomal location and frequency of inversions can direct dissection of phenotypic changes following chromosomal inversions. The framework used here can be transferred to other bacterial species to characterise chromosomal inversions.

Novel 4-chlorophenoxyacetate dioxygenase-mediated phenoxyalkanoic acid herbicides initial catabolism in Cupriavidus sp. DL-D2

Microbial metabolism is an important driving force for the elimination of 4-chlorophenoxyacetic acid residues in the environment. The α-Ketoglutarate-dependent dioxygenase (TfdA) or 2,4-D oxygenase (CadAB) catalyzes the cleavage of the aryl ether bond of 4-chlorophenoxyacetic acid to 4-chlorophenol, which is one of the important pathways for the initial metabolism of 4-chlorophenoxyacetic acid by microorganisms. However, strain Cupriavidus sp. DL-D2 could utilize 4-chlorophenoxyacetic acid but not 4-chlorophenol for growth. This scarcely studied degradation pathway may involve novel enzymes that has not yet been characterized. Here, a gene cluster (designated cpd) responsible for the catabolism of 4-chlorophenoxyacetic acid in strain DL-D2 was cloned and identified, and the dioxygenase CpdA/CpdB responsible for the initial degradation of 4-chlorophenoxyacetic acid was successfully expressed, which could catalyze the conversion of 4-chlorphenoxyacetic acid to 4-chlorocatechol. Then, an aromatic cleavage enzyme CpdC further converts 4-chlorocatechol into 3-chloromuconate. The results of substrate degradation experiments showed that CpdA/CpdB could also degrade 3-chlorophenoxyacetic acid and phenoxyacetic acid, and homologous cpd gene clusters were widely discovered in microbial genomes. Our findings revealed a novel degradation mechanism of 4-chlorophenoxyacetic acid at the molecular level.

Comparison of Bordetella species identification among differing rt-PCR assays in the United States

In the United States, the general laboratory method for diagnosing pertussis, caused by Bordetella pertussis, is real-time PCR (rt-PCR) targeting insertion sequence 481 (IS481). Other Bordetella species (parapertussis, holmesii, and bronchiseptica) can also cause a pertussis-like syndrome, and some commercial laboratory assays include the insertion sequence 1001 (pIS1001) that can detect B. parapertussis/B. bronchiseptica (BppBb). Because IS481 exists in B. pertussis and B. holmesii, current commercial assays cannot differentiate these two species. We used a multiplex rt-PCR assay containing species-specific targets to Bordetella to evaluate clinical specimens detected as B. pertussis/B. holmesii (BpBh) or BppBb by commercial laboratories. A sample of 3,984 clinical specimens positive for IS481 or pIS1001 from two commercial laboratories during 2012-2019 were re-tested at CDC. Agreement of Bordetella species between the CDC and commercial laboratory assays, and the proportion of commercial laboratory specimens that were non-B. pertussis by CDC's assay was assessed. Overall agreement in Bordetella species detection and identification between the CDC and commercial lab assays was 85%. Agreement for identifying B. pertussis was 87% for 3,663 BpBh specimens and 98% for identifying B. parapertussis in 310 BppBb specimens. CDC's assay detected B. holmesii in 55/3,984 (1.4%) specimens. Most discrepant results (410/490, 82%) were BpBh specimens interpreted as indeterminate B. pertussis at CDC. We found a small portion of B. holmesii in a sample of IS481-positive clinical specimens originally identified by commercial laboratory rt-PCR assays, suggesting that commercial PCR assays are a reliable diagnostic tool for correctly identifying Bordetella species in most patients with suspected pertussis.

IMPORTANCE

When testing specimens collected from patients with suspected pertussis, large-scale commercial laboratories in the United States employ an IS481-based assay that cannot differentiate between Bordetella pertussis and Bordetella holmseii. The level of B. holmesii causing pertussis-like illness in the United States is not well-understood given that only B. pertussis is nationally notifiable. After re-testing with a multiplex, species-specific rt-PCR assay, our data show low levels of B. holmesii identified in a sample of IS481-positive clinical specimens originally identified by commercial laboratory rt-PCR assays. These results reinforce the validity of large-scale commercial rt-PCR testing as a reliable diagnostic tool for pertussis in the United States.

Cytochrome oxidase requirements in Bordetella reveal insights into evolution towards life in the mammalian respiratory tract

Little is known about oxygen utilization during infection by bacterial respiratory pathogens. The classical Bordetella species, including B. pertussis, the causal agent of human whooping cough, and B. bronchiseptica, which infects nearly all mammals, are obligate aerobes that use only oxygen as the terminal electron acceptor for electron transport-coupled oxidative phosphorylation. B. bronchiseptica, which occupies many niches, has eight distinct cytochrome oxidase-encoding loci, while B. pertussis, which evolved from a B. bronchiseptica-like ancestor but now survives exclusively in and between human respiratory tracts, has only three functional cytochrome oxidase-encoding loci: cydAB1, ctaCDFGE1, and cyoABCD1. To test the hypothesis that the three cytochrome oxidases encoded within the B. pertussis genome represent the minimum number and class of cytochrome oxidase required for respiratory infection, we compared B. bronchiseptica strains lacking one or more of the eight possible cytochrome oxidases in vitro and in vivo. No individual cytochrome oxidase was required for growth in ambient air, and all three of the cytochrome oxidases conserved in B. pertussis were sufficient for growth in ambient air and low oxygen. Using a high-dose, large-volume persistence model and a low-dose, small-volume establishment of infection model, we found that B. bronchiseptica producing only the three B. pertussis-conserved cytochrome oxidases was indistinguishable from the wild-type strain for infection. We also determined that CyoABCD1 is sufficient to cause the same level of bacterial burden in mice as the wild-type strain and is thus the primary cytochrome oxidase required for murine infection, and that CydAB1 and CtaCDFGE1 fulfill auxiliary roles or are important for aspects of infection we have not assessed, such as transmission. Our results shed light on the environment at the surface of the ciliated epithelium, respiration requirements for bacteria that colonize the respiratory tract, and the evolution of virulence in bacterial pathogens.

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.

Genome-based prediction of cross-protective, HLA-DR-presented epitopes as putative vaccine antigens for multiple Bordetella species

IMPORTANCE

Pertussis, caused by Bordetella pertussis, can cause debilitating respiratory symptoms, so whole-cell pertussis vaccines (wPVs) were introduced in the 1940s. However, reactogenicity of wPV necessitated the development of acellular pertussis vaccines (aPVs) that were introduced in the 1990s. Since then, until the COVID-19 pandemic began, reported pertussis incidence was increasing, suggesting that aPVs do not induce long-lasting immunity and may not effectively prevent transmission. Additionally, aPVs do not provide protection against other Bordetella species that are observed during outbreaks. The significance of this work is in determining potential new vaccine antigens for multiple Bordetella species that are predicted to elicit long-term immune responses. Genome-based approaches have aided the development of novel vaccines; here, these methods identified Bordetella vaccine candidates that may be cross-protective and predicted to induce strong memory responses. These targets can lead to an improved vaccine with a strong safety profile while also strengthening the longevity of the immune response.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Genomic characterization of Bordetella pertussis in South Africa, 2015-2019

Pertussis remains a public health concern in South Africa, with an increase in reported cases and outbreaks in recent years. Whole genome sequencing was performed on 32 Bordetella pertussis isolates sourced from three different surveillance programmes in South Africa between 2015 and 2019. Genome sequences were characterized using multilocus sequence typing, vaccine antigen genes (ptxP, ptxA, ptxB, prn and fimH) and overall genome structure. All isolates were sequence type 2 and harboured the pertussis toxin promoter allele ptxP3. The dominant genotype was ptxP3-ptxA1-ptxB2-prn2-fimH2 (31/32, 96.9 %), with no pertactin-deficient or other mutations in vaccine antigen genes identified. Amongst 21 isolates yielding closed genome assemblies, eight distinct genome structures were detected, with 61.9 % (13/21) of the isolates exhibiting three predominant structures. Increases in case numbers are probably not due to evolutionary changes in the genome but possibly due to other factors such as the cyclical nature of B. pertussis disease, waning immunity due to the use of acellular vaccines and/or population immunity gaps.

Bordetella spp. block eosinophil recruitment to suppress the generation of early mucosal protection

Bordetella spp. are respiratory pathogens equipped with immune evasion mechanisms. We previously characterized a Bordetella bronchiseptica mutant (RB50ΔbtrS) that fails to suppress host responses, leading to rapid clearance and long-lasting immunity against reinfection. This work revealed eosinophils as an exclusive requirement for RB50ΔbtrS clearance. We also show that RB50ΔbtrS promotes eosinophil-mediated B/T cell recruitment and inducible bronchus-associated lymphoid tissue (iBALT) formation, with eosinophils being present throughout iBALT for Th17 and immunoglobulin A (IgA) responses. Finally, we provide evidence that XCL1 is critical for iBALT formation but not maintenance, proposing a novel role for eosinophils as facilitators of adaptive immunity against B. bronchiseptica. RB50ΔbtrS being incapable of suppressing eosinophil effector functions illuminates active, bacterial targeting of eosinophils to achieve successful persistence and reinfection. Overall, our discoveries contribute to understanding cellular mechanisms for use in future vaccines and therapies against Bordetella spp. and extension to other mucosal pathogens.

Pasteurian Contributions to the Study of Bordetella pertussis Toxins

As a tribute to Louis Pasteur on the occasion of the 200th anniversary of his birth, this article summarizes the main contributions of scientists from Pasteur Institutes to the current knowledge of toxins produced by Bordetella pertussis. The article therefore focuses on publications authored by researchers from Pasteur Institutes and is not intended as a systematic review of B. pertussis toxins. Besides identifying B. pertussis as the causative agent of whooping cough, Pasteurians have made several major contributions with respect to the structure–function relationship of the Bordetella lipo-oligosaccharide, adenylyl cyclase toxin and pertussis toxin. In addition to contributing to the understanding of these toxins’ mechanisms at the molecular and cellular levels and their role in pathogenesis, scientists at Pasteur Institutes have also exploited potential applications of the gathered knowledge of these toxins. These applications range from the development of novel tools to study protein–protein interactions over the design of novel antigen delivery tools, such as prophylactic or therapeutic vaccine candidates against cancer and viral infection, to the development of a live attenuated nasal pertussis vaccine. This scientific journey from basic science to applications in the field of human health matches perfectly with the overall scientific objectives outlined by Louis Pasteur himself.

Impact of Salmonella genome rearrangement on gene expression

In addition to nucleotide variation, many bacteria also undergo changes at a much larger scale via rearrangement of their genome structure (GS) around long repeat sequences. These rearrangements result in genome fragments shifting position and/or orientation in the genome without necessarily affecting the underlying nucleotide sequence. To date, scalable techniques have not been applied to GS identification, so it remains unclear how extensive this variation is and the extent of its impact upon gene expression. However, the emergence of multiplexed, long-read sequencing overcomes the scale problem, as reads of several thousand bases are routinely produced that can span long repeat sequences to identify the flanking chromosomal DNA, allowing GS identification. Genome rearrangements were generated in Salmonella enterica serovar Typhi through long-term culture at ambient temperature. Colonies with rearrangements were identified via long-range PCR and subjected to long-read nanopore sequencing to confirm genome variation. Four rearrangements were investigated for differential gene expression using transcriptomics. All isolates with changes in genome arrangement relative to the parent strain were accompanied by changes in gene expression. Rearrangements with similar fragment movements demonstrated similar changes in gene expression. The most extreme rearrangement caused a large imbalance between the origin and terminus of replication and was associated with differential gene expression as a factor of distance moved toward or away from the origin of replication. Genome structure variation may provide a mechanism through which bacteria can quickly adapt to new environments and warrants routine assessment alongside traditional nucleotide-level measures of variation.

Replication-associated inversions are the dominant form of bacterial chromosome structural variation

The structural arrangements of bacterial chromosomes vary widely between closely related species and can result in significant phenotypic outcomes. The appearance of large-scale chromosomal inversions that are symmetric relative to markers for the origin of replication (OriC) has been previously observed; however, the overall prevalence of replication-associated structural rearrangements (RASRs) in bacteria and their causal mechanisms are currently unknown. Here, we systematically identify the locations of RASRs in species with multiple complete-sequenced genomes and investigate potential mediating biological mechanisms. We found that 247 of 313 species contained sequences with at least one large (>50 Kb) inversion in their sequence comparisons, and the aggregated inversion distances away from symmetry were normally distributed with a mean of zero. Many inversions that were offset from dnaA were found to be centered on a different marker for the OriC Instances of flanking repeats provide evidence that breaks formed during the replication process could be repaired to opposing positions. We also found a strong relationship between the later stages of replication and the range in distance variation from symmetry.

Genome Characteristic of Bordetella parapertussis Isolated from Iran

Pertussis also known as whooping cough is a respiratory infection in humans particularly with severe symptoms in infants and usually caused by Bordetella pertussis. However, Bordetella parapertussis can also cause a similar clinical syndrome. During 2012 to 2015, from nasal swabs sent from different provinces to the pertussis reference laboratory of Pasture Institute of Iran for pertussis confirmation, seven B. parapertussis isolates were identified by bacterial culture, biochemical tests, and the presence of IS1001 insertion in the genome. The expression of pertactin (Prn) as one the major virulence factor for bacterial adhesion was investigated using western blot. Moreover, the genomic characteristic of one recently collected isolate, IRBP134, from a seven-month infant was investigated using Illumina NextSeq sequencing protocol. The results revealed the genome with G+C content 65% and genome size 4.7 Mbp. A total of 81 single nucleotide polymorphisms and 13 short insertions and deletions were found in the genome compared to the B. parapertussis 12822 as a reference genome showing ongoing evolutionary changes. A phylogeny relationship of IRBP134 was also investigated using global B. parapertussis available genomes.

Characterization of Bordetella pertussis Strains Isolated from India

Despite high level vaccination and the availability of two different types of vaccines, whole cell (wP) and acellular vaccines (aP), the resurgence of pertussis has been reported in many countries. Antigenic variation within circulating and vaccine strains is the most documented reason reported for the resurgence of pertussis. Research on genetic divergence among circulating and vaccine strains has largely been reported in countries using aP vaccines. There are inadequate data available for antigenic variation in B. pertussis from wP-using countries. India has used wP for more than 40 years in their primary immunization program. The present study reports five clinical isolates of B. pertussis from samples of pediatric patients with pertussis symptoms observed in India. Genotypic and phenotypic characterization of clinical isolates were performed by serotyping, genotyping, whole genome analyses and comparative genomics. All clinical isolates showed serotype 1, 2 and 3 based on the presence of fimbriae 2 and 3. Genotyping showed genetic similarities in allele types for five aP genes within vaccine strains and clinical isolates reported from India. The presence of the ptxP3 genotype was observed in two out of five clinical isolates. Whole-genome sequencing was performed for clinical isolates using the hybrid strategy of combining Illumina (short reads) and oxford nanopore (long reads) sequencing strategies. Clinical isolates (n = 5) and vaccine strains (n = 7) genomes of B. pertussis from India were compared with 744 B. pertussis closed genomes available in the public databases. The phylogenomic comparison of B. pertussis genomes reported from India will be advantageous in better understanding pertussis resurgence reported globally with respect to pathogen adaptation.

Analysis of the Taxonomy, Synteny, and Virulence Factors for Soft Rot Pathogen Pectobacterium aroidearum in Amorphophallus konjac Using Comparative Genomics

Bacterial soft rot is a devastating disease for a wide range of crops, vegetables, and ornamental plants including konjac (Amorphophallus konjac). However, the pangenome and genomic plasticity of the konjac soft rot pathogens is little explored. In this study, we reported the complete genome sequences of 11 bacterial isolates that can cause typical soft rot symptoms in konjac by in vitro and in vivo pathogenicity tests. Based on in silico DNA-DNA hybridization, average nucleotide identity and phylogenomic analysis, all 11 isolates were determined to be Pectobacterium aroidearum. In addition, synteny analysis of these genomes revealed considerable chromosomal inversions, one of which is triggered by homologous recombination of ribose operon. Pangenome analysis and COG enrichment analysis showed that the pangenome of P. aroidearum is open and that accessory genes are enriched in replication, recombination, and repair. Variations in type IV secretion system and type VI secretion system were found, while plant cell wall degrading enzymes were conserved. Furthermore, sequence analyses also provided evidence for the presence of a type V secretion system in Pectobacterium. These findings advance our understanding of the pathogenicity determinants, genomic plasticity, and evolution of P. aroidearum.

A comprehensive resource for Bordetella genomic epidemiology and biodiversity studies

The genus Bordetella includes bacteria that are found in the environment and/or associated with humans and other animals. A few closely related species, including Bordetella pertussis, are human pathogens that cause diseases such as whooping cough. Here, we present a large database of Bordetella isolates and genomes and develop genotyping systems for the genus and for the B. pertussis clade. To generate the database, we merge previously existing databases from Oxford University and Institut Pasteur, import genomes from public repositories, and add 83 newly sequenced B. bronchiseptica genomes. The public database currently includes 2582 Bordetella isolates and their provenance data, and 2085 genomes ( https://bigsdb.pasteur.fr/bordetella/ ). We use core-genome multilocus sequence typing (cgMLST) to develop genotyping systems for the whole genus and for B. pertussis, as well as specific schemes to define antigenic, virulence and macrolide resistance profiles. Phylogenetic analyses allow us to redefine evolutionary relationships among known Bordetella species, and to propose potential new species. Our database provides an expandable resource for genotyping of environmental and clinical Bordetella isolates, thus facilitating evolutionary and epidemiological research on whooping cough and other Bordetella infections.

Tatajuba: exploring the distribution of homopolymer tracts

Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialized bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and help predict their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source licence GNU GPLv3.

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.

Comparative genomics of Bordetella pertussis isolates from New Zealand, a country with an uncommonly high incidence of whooping cough

Whooping cough, the respiratory disease caused by Bordetella pertussis, has undergone a wide-spread resurgence over the last several decades. Previously, we developed a pipeline to assemble the repetitive B. pertussis genome into closed sequences using hybrid nanopore and Illumina sequencing. Here, this sequencing pipeline was used to conduct a more high-throughput, longitudinal screen of 66 strains isolated between 1982 and 2018 in New Zealand. New Zealand has a higher incidence of whooping cough than many other countries; usually at least twice as many cases per 100000 people as the USA and UK and often even higher, despite similar rates of vaccine uptake. To the best of our knowledge, these strains are the first New Zealand B. pertussis isolates to be sequenced. The analyses here show that, on the whole, genomic trends in New Zealand B. pertussis isolates, such as changing allelic profile in vaccine-related genes and increasing pertactin deficiency, have paralleled those seen elsewhere in the world. At the same time, phylogenetic comparisons of the New Zealand isolates with global isolates suggest that a number of strains are circulating in New Zealand, which cluster separately from other global strains, but which are closely related to each other. The results of this study add to a growing body of knowledge regarding recent changes to the B. pertussis genome, and are the first genetic investigation into B. pertussis isolates from New Zealand.

Bilateral symmetry of linear streptomycete chromosomes

Here, we characterize an uncommon set of telomeres from Streptomyces rimosus ATCC 10970, the parental strain of a lineage of one of the earliest-discovered antibiotic producers. Following the closure of its genome sequence, we compared unusual telomeres from this organism with the other five classes of replicon ends found amongst streptomycetes. Closed replicons of streptomycete chromosomes were organized with respect to their phylogeny and physical orientation, which demonstrated that different telomeres were not associated with particular clades and are likely shared amongst different strains by plasmid-driven horizontal gene transfer. Furthermore, we identified a ~50 kb origin island with conserved synteny that is located at the core of all streptomycete chromosomes and forms an axis around which symmetrical chromosome inversions can take place. Despite this chromosomal bilateral symmetry, a bias in parS sites to the right of oriC is maintained across the family Streptomycetaceae and suggests that the formation of ParB/parS nucleoprotein complexes on the right replichore is a conserved feature in streptomycetes. Consequently, our studies reveal novel features of linear bacterial replicons that, through their manipulation, may lead to improvements in growth and productivity of this important industrial group of bacteria.

Pathogenicity and virulence of Bordetella pertussis and its adaptation to its strictly human host

The highly contagious whooping cough agent Bordetella pertussis has evolved as a human-restricted pathogen from a progenitor which also gave rise to Bordetella parapertussis and Bordetella bronchiseptica. While the latter colonizes a broad range of mammals and is able to survive in the environment, B. pertussis has lost its ability to survive outside its host through massive genome decay. Instead, it has become a highly successful human pathogen by the acquisition of tightly regulated virulence factors and evolutionary adaptation of its metabolism to its particular niche. By the deployment of an arsenal of highly sophisticated virulence factors it overcomes many of the innate immune defenses. It also interferes with vaccine-induced adaptive immunity by various mechanisms. Here, we review data from invitro, human and animal models to illustrate the mechanisms of adaptation to the human respiratory tract and provide evidence of ongoing evolutionary adaptation as a highly successful human pathogen.

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.

Pertussis among patients with clinically compatible illness in the Amhara Regional State, Ethiopia

BACKGROUND

Pertussis is an acute respiratory tract disease caused by Bordetella pertussis. In 2014, 24.1 million pertussis cases, resulting in 160,700 deaths, were estimated to have occurred worldwide. This study aimed to determine the epidemiology of pertussis among patients with clinically compatible illness who visited selected hospitals in the Amhara Regional State of Ethiopia.

METHODS

A cross-sectional study design was used to review pertussis patients with clinically compatible illness. Nasopharyngeal swabs were collected from 515 patients from July 2018 through February 2019. DNA was extracted from all nasopharyngeal swabs and samples were analyzed using real-time (RT-) PCR. Crude and adjusted odds ratios with corresponding 95% confidence intervals were estimated using bivariable and multivariable logistic regression analysis, respectively.

RESULTS

The overall prevalence of Bordetella species among the study participants was 156 of 515 (30.3%) [95% CI = 26.4-34.6] as determined by Bordetella RT-PCR, including: 65 (41.7%) B. pertussis, 89 (57.1%) indeterminate B. pertussis, one (0.6%) Bordetella holmesii and one (0.6%) Bordetella parapertussis.

CONCLUSIONS

This study found that pertussis is potentially endemic and a common health problem among patients visiting health institutions in the Amhara Regional State of Ethiopia. More data regarding pertussis in Ethiopia could inform development of effective prevention strategies.

A Mutation Upstream of the rplN-rpsD Ribosomal Operon Downregulates Bordetella pertussis Virulence Factor Production without Compromising Bacterial Survival within Human Macrophages

The BvgS/BvgA two-component system controls expression of ∼550 genes of Bordetella pertussis, of which, ∼245 virulence-related genes are positively regulated by the BvgS-phosphorylated transcriptional regulator protein BvgA (BvgA∼P). We found that a single G-to-T nucleotide transversion in the 5'-untranslated region (5'-UTR) of the rplN gene enhanced transcription of the ribosomal protein operon and of the rpoA gene and provoked global dysregulation of B. pertussis genome expression. This comprised overproduction of the alpha subunit (RpoA) of the DNA-dependent RNA polymerase, downregulated BvgA and BvgS protein production, and impaired production and secretion of virulence factors by the mutant. Nonetheless, the mutant survived like the parental bacteria for >2 weeks inside infected primary human macrophages and persisted within infected mouse lungs for a longer period than wild-type B. pertussis These observations suggest that downregulation of virulence factor production by bacteria internalized into host cells may enable persistence of the whooping cough agent in the airways.IMPORTANCE We show that a spontaneous mutation that upregulates transcription of an operon encoding ribosomal proteins and causes overproduction of the downstream-encoded α subunit (RpoA) of RNA polymerase causes global effects on gene expression levels and proteome composition of Bordetella pertussis Nevertheless, the resulting important downregulation of the BvgAS-controlled expression of virulence factors of the whooping cough agent did not compromise its capacity to persist for prolonged periods inside primary human macrophage cells, and it even enhanced its capacity to persist in infected mouse lungs. These observations suggest that the modulation of BvgAS-controlled expression of virulence factors may occur also during natural infections of human airways by Bordetella pertussis and may possibly account for long-term persistence of the pathogen within infected cells of the airways.

Bordetella Type III Secretion Injectosome and Effector Proteins

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