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

Genome-Based Analysis of Aspergillus niger Aggregate Species from China and Their Potential for Fumonisin B2 and Ochratoxin A Production

Based on entire genome sequencing, this study focused on the classification of Aspergillus niger aggregation species and investigated their potential for fumonisin B₂ (FB₂) and ochratoxin A (OTA) production. In the current study, 22 strains were used, namely 17 A. niger strains, four A. welwitschiae strains, and one A. lacticoffeatus (a synonym of A. niger) strain. Traditional multigene phylogenetic analysis, average nucleotide identity analysis (ANI), and the whole-genome single-nucleotide polymorphism (SNP) analyses were used to reconfirm the taxonomic status of A. niger, A. welwitschiae, and A. lacticoffeatus. The ability of A. niger to produce FB₂ and OTA on five culture substrates was determined, and the association between FB₂ and OTA gene clusters and toxin-producing abilities was explored. The results revealed that the ANI method could distinguish A. niger from A. welwitschiae, with an ANI value of < 98%. The SNP-based phylogenetic analysis suggested that A. niger and A. welwitschiae were two independent phylogenetic species. The ANI, SNP, and multigene phylogenetic analysis supported previous findings that A. lacticoffeatus was a synonymous species of A. niger. Aspergillus niger strains exhibited the varied potential of producing FB₂ and OTA on different culture media. The A. niger genome sequence analysis revealed no significant difference in fumonisin gene clusters between FB₂-nonproducing isolates and FB₂-producing isolates, and the integrity of the ochratoxin biosynthesis genes cluster was clearly associated with OTA production. In conclusion, gene sequencing can be useful in assessing A. niger's ability to produce OTA, but it cannot reliably predict its ability to produce FB₂.

Rapid and accurate SNP genotyping of clonal bacterial pathogens with BioHansel

Hierarchical genotyping approaches can provide insights into the source, geography and temporal distribution of bacterial pathogens. Multiple hierarchical SNP genotyping schemes have previously been developed so that new isolates can rapidly be placed within pre-computed population structures, without the need to rebuild phylogenetic trees for the entire dataset. This classification approach has, however, seen limited uptake in routine public health settings due to analytical complexity and the lack of standardized tools that provide clear and easy ways to interpret results. The BioHansel tool was developed to provide an organism-agnostic tool for hierarchical SNP-based genotyping. The tool identifies split k-mers that distinguish predefined lineages in whole genome sequencing (WGS) data using SNP-based genotyping schemes. BioHansel uses the Aho-Corasick algorithm to type isolates from assembled genomes or raw read sequence data in a matter of seconds, with limited computational resources. This makes BioHansel ideal for use by public health agencies that rely on WGS methods for surveillance of bacterial pathogens. Genotyping results are evaluated using a quality assurance module which identifies problematic samples, such as low-quality or contaminated datasets. Using existing hierarchical SNP schemes for Mycobacterium tuberculosis and Salmonella Typhi, we compare the genotyping results obtained with the k-mer-based tools BioHansel and SKA, with those of the organism-specific tools TBProfiler and genotyphi, which use gold-standard reference-mapping approaches. We show that the genotyping results are fully concordant across these different methods, and that the k-mer-based tools are significantly faster. We also test the ability of the BioHansel quality assurance module to detect intra-lineage contamination and demonstrate that it is effective, even in populations with low genetic diversity. We demonstrate the scalability of the tool using a dataset of ~8100 S. Typhi public genomes and provide the aggregated results of geographical distributions as part of the tool’s output. BioHansel is an open source Python 3 application available on PyPI and Conda repositories and as a Galaxy tool from the public Galaxy Toolshed. In a public health context, BioHansel enables rapid and high-resolution classification of bacterial pathogens with low genetic diversity.

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

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

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

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

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

Whole-genome analysis of probiotic product isolates reveals the presence of genes related to antimicrobial resistance, virulence factors, and toxic metabolites, posing potential health risks

Background

Safety issues of probiotic products have been reported frequently in recent years. Ten bacterial strains isolated from seven commercial probiotic products on market were evaluated for their safety, by whole-genome analysis.

Results

We found that the bacterial species of three probiotic products were incorrectly labeled. Furthermore, six probiotic product isolates (PPS) contained genes for the production of toxic metabolites, while another three strains contained virulence genes, which might pose a potential health risk. In addition, three of them have drug-resistance genes, among which two strains potentially displayed multidrug resistance. One isolate has in silico predicted transferable genes responsible for toxic metabolite production, and they could potentially transfer to human gut microflora or environmental bacteria. Isolates of Lactobacillus rhamnosus and Bifidobacterium animalis subsp. lactis are associated with low risk for human consumption. Based on a comparative genome analysis, we found that the isolated Enterococcus faecium TK-P5D clustered with a well-defined probiotic strain, while E. faecalis TK-P4B clustered with a pathogenic strain.

Conclusions

Our work clearly illustrates that whole-genome analysis is a useful method to evaluate the quality and safety of probiotic products. Regulatory quality control and stringent regulations on probiotic products are needed to ensure safe consumption and protect human health.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07539-9.

Rapid and simple SNP genotyping for Bordetella pertussis epidemic strain MT27 based on a multiplexed single-base extension assay

Multilocus variable-number tandem repeat analysis (MLVA) is widely used for genotyping of Bordetella pertussis, the causative bacteria for pertussis. However, MLVA genotyping is losing its discriminate power because prevalence of the epidemic MT27 strain (MLVA-27) is increasing worldwide. To address this, we developed a single nucleotide polymorphism (SNP) genotyping method for MT27 based on multiplexed single-base extension (SBE) assay. A total of 237 MT27 isolates collected in Japan during 1999–2018 were genotyped and classified into ten SNP genotypes (SG1 to SG10) with a Simpson’s diversity index (DI) of 0.79 (95% CI 0.76–0.82). Temporal trends showed a marked increase in the genotypic diversity in the 2010s: Simpson’s DI was zero in 1999–2004, 0.16 in 2005–2009, 0.83 in 2010–2014, and 0.76 in 2015–2018. This indicates that the SNP genotyping is applicable to the recently circulating MT27 strain. Additionally, almost all outbreak-associated MT27 isolates were classified into the same SNP genotypes for each outbreak. Multiplexed SBE assay allows for rapid and simple genotyping, indicating that the SNP genotyping can potentially be a useful tool for subtyping the B. pertussis MT27 strain in routine surveillance and outbreak investigations.

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.

New Data on Vaccine Antigen Deficient Bordetella pertussis Isolates

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

Genetic Variation of Bordetella pertussis in Austria

In Austria, vaccination coverage against Bordetella pertussis infections during infancy is estimated at around 90%. Within the last years, however, the number of pertussis cases has increased steadily, not only in children but also in adolescents and adults, indicating both insufficient herd immunity and vaccine coverage. Waning immunity in the host and/or adaptation of the bacterium to the immunised hosts could contribute to the observed re-emergence of pertussis. In this study we therefore addressed the genetic variability in B. pertussis strains from several Austrian cities. Between the years 2002 and 2008, 110 samples were collected from Vienna (n = 32), Linz (n = 63) and Graz (n = 15) by nasopharyngeal swabs. DNA was extracted from the swabs, and bacterial sequence polymorphisms were examined by MLVA (multiple-locus variable number of tandem repeat analysis) (n = 77), by PCR amplification and conventional Sanger sequencing of the polymorphic regions of the prn (pertactin) gene (n = 110), and by amplification refractory mutation system quantitative PCR (ARMS-qPCR) (n = 110) to directly address polymorphisms in the genes encoding two pertussis toxin subunits (ptxA and ptxB), a fimbrial adhesin (fimD), tracheal colonisation factor (tcfA), and the virulence sensor protein (bvgS). Finally, the ptxP promoter region was screened by ARMS-qPCR for the presence of the ptxP3 allele, which has been associated with elevated production of pertussis toxin. The MLVA analysis revealed the highest level of polymorphisms with an absence of MLVA Type 29, which is found outside Austria. Only Prn subtypes Prn1/7, Prn2 and Prn3 were found with a predominance of the non-vaccine type Prn2. The analysis of the ptxA, ptxB, fimD, tcfA and bvgS polymorphisms showed a genotype mixed between the vaccine strain Tohama I and a clinical isolate from 2006 (L517). The major part of the samples (93%) displayed the ptxP3 allele. The consequences for the vaccination strategy are discussed.

Studying Bordetella pertussis populations by use of SNPeX, a simple high-throughput single nucleotide polymorphism typing method

Large outbreaks of pertussis occur despite vaccination. A first step in the analyses of outbreaks is strain typing. However, the typing of Bordetella pertussis, the causative agent of pertussis, is problematic because the available assays are insufficiently discriminatory, not unequivocal, time-consuming, and/or costly. Here, we describe a single nucleotide primer extension assay for the study of B. pertussis populations, SNPeX (single nucleotide primer extension), which addresses these problems. The assay is based on the incorporation of fluorescently labeled dideoxynucleotides (ddNTPs) at the 3' end of allele-specific poly(A)-tailed primers and subsequent analysis with a capillary DNA analyzer. Each single nucleotide polymorphism (SNP) primer has a specific length, and as a result, up to 20 SNPs can be determined in one SNPeX reaction. Importantly, PCR amplification of target DNA is not required. We selected 38 SNPeX targets from the whole-genome sequencing data of 74 B. pertussis strains collected from across the world. The SNPeX-based phylogenetic trees preserved the general tree topology of B. pertussis populations based on whole-genome sequencing, with a minor loss of details. We envisage a strategy whereby SNP types (SnpTs) are quickly identified with the SNPeX assay during an outbreak, followed by whole-genome sequencing (WGS) of a limited number of isolates representing predominant SnpTs and the incorporation of novel SNPs in the SNPeX assay. The flexibility of the SNPeX assay allows the method to evolve along with the pathogen, making it a promising method for studying outbreaks of B. pertussis and other pathogens.

Analysis of Bordetella pertussis clinical isolates circulating in European countries during the period 1998-2012

Despite more than 50 years of vaccination, pertussis is still an endemic disease, with regular epidemic outbreaks. With the exception of Poland, European countries have replaced whole-cell vaccines (WCVs) by acellular vaccines (ACVs) in the 1990s. Worldwide, antigenic divergence in vaccine antigens has been found between vaccine strains and circulating strains. In this work, 466 Bordetella pertussis isolates collected in the period 1998–2012 from 13 European countries were characterised by multi-locus antigen sequence typing (MAST) of the pertussis toxin promoter (ptxP) and of the genes coding for proteins used in the ACVs: pertussis toxin (Ptx), pertactin (Prn), type 2 fimbriae (Fim2) and type 3 fimbriae (Fim3). Isolates were further characterised by fimbrial serotyping, multi-locus variable-number tandem repeat analysis (MLVA) and pulsed-field gel electrophoresis (PFGE). The results showed a very similar B. pertussis population for 12 countries using ACVs, while Poland, which uses a WCV, was quite distinct, suggesting that ACVs and WCVs select for different B. pertussis populations. This study forms a baseline for future studies on the effect of vaccination programmes on B. pertussis populations.

Bordetella pertussis naturally occurring isolates with altered lipooligosaccharide structure fail to fully mature human dendritic cells

Bordetella pertussis is a Gram-negative bacterium and the causative agent of whooping cough. Despite high vaccination coverage, outbreaks are being increasingly reported worldwide. Possible explanations include adaptation of this pathogen, which may interfere with recognition by the innate immune system. Here, we describe innate immune recognition and responses to different B. pertussis clinical isolates. By using HEK-Blue cells transfected with different pattern recognition receptors, we found that 3 out of 19 clinical isolates failed to activate Toll-like receptor 4 (TLR4). These findings were confirmed by using the monocytic MM6 cell line. Although incubation with high concentrations of these 3 strains resulted in significant activation of the MM6 cells, it was found to occur mainly through interaction with TLR2 and not through TLR4. When using live bacteria, these 3 strains also failed to activate TLR4 on HEK-Blue cells, and activation of MM6 cells or human monocyte-derived dendritic cells was significantly lower than activation induced by the other 16 strains. Mass spectrum analysis of the lipid A moieties from these 3 strains indicated an altered structure of this molecule. Gene sequence analysis revealed mutations in genes involved in lipid A synthesis. Findings from this study indicate that B. pertussis isolates that do not activate TLR4 occur naturally and that this phenotype may give this bacterium an advantage in tempering the innate immune response and establishing infection. Knowledge on the strategies used by this pathogen in evading the host immune response is essential for the improvement of current vaccines or for the development of new ones.

Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing

We explored genetic variation by sequencing a selection of 84 tomato accessions and related wild species representative of the Lycopersicon, Arcanum, Eriopersicon and Neolycopersicon groups, which has yielded a huge amount of precious data on sequence diversity in the tomato clade. Three new reference genomes were reconstructed to support our comparative genome analyses. Comparative sequence alignment revealed group-, species- and accession-specific polymorphisms, explaining characteristic fruit traits and growth habits in the various cultivars. Using gene models from the annotated Heinz 1706 reference genome, we observed differences in the ratio between non-synonymous and synonymous SNPs (dN/dS) in fruit diversification and plant growth genes compared to a random set of genes, indicating positive selection and differences in selection pressure between crop accessions and wild species. In wild species, the number of single-nucleotide polymorphisms (SNPs) exceeds 10 million, i.e. 20-fold higher than found in most of the crop accessions, indicating dramatic genetic erosion of crop and heirloom tomatoes. In addition, the highest levels of heterozygosity were found for allogamous self-incompatible wild species, while facultative and autogamous self-compatible species display a lower heterozygosity level. Using whole-genome SNP information for maximum-likelihood analysis, we achieved complete tree resolution, whereas maximum-likelihood trees based on SNPs from ten fruit and growth genes show incomplete resolution for the crop accessions, partly due to the effect of heterozygous SNPs. Finally, results suggest that phylogenetic relationships are correlated with habitat, indicating the occurrence of geographical races within these groups, which is of practical importance for Solanum genome evolution studies.

Investigations into the emergence of pertactin-deficient Bordetella pertussis isolates in six European countries, 1996 to 2012

Pathogen adaptation has been proposed to contribute to the resurgence of pertussis. A striking recent example is the emergence of isolates deficient in the vaccine component pertactin (Prn). This study explores the emergence of such Prn-deficient isolates in six European countries. During 2007 to 2009, 0/83 isolates from the Netherlands, 0/18 from the United Kingdom, 0/17 Finland, 0/23 Denmark, 4/99 Sweden and 5/20 from Norway of the isolates collected were Prn-deficient. In the Netherlands and Sweden, respectively 4/146 and 1/8 were observed in a later period (2010–12). The Prn-deficient isolates were genetically diverse and different mutations were found to inactivate the prn gene. These are indications that Prn-deficiency is subject to positive selective pressure. We hypothesise that the switch from whole cell to acellular pertussis vaccines has affected the balance between ‘costs and benefits’ of Prn production by Bordetella pertussis to the extent that isolates that do not produce Prn are able to expand. The absence of Prn-deficient isolates in some countries may point to ways to prevent or delay the spread of Prn-deficient strains. In order to substantiate this hypothesis, trends in the European B. pertussis population should be monitored continuously.

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.

Bacterial strain typing

Over the course of the past several decades, rapid advancements in molecular technologies have revolutionized the practice of public health microbiology, and have fundamentally changed the nature, accuracy, and timeliness of laboratory data for outbreak investigation and response. Whole-genome sequencing, in particular, is becoming an increasingly feasible and cost-effective approach for near real-time high-resolution strain typing, genomic characterization, and comparative analyses. This review discusses the current state of the art in bacterial strain typing for outbreak investigation and infectious disease surveillance, and the impact of emerging genomic technologies on the field of public health microbiology.

Differentially expressed genes in Bordetella pertussis strains belonging to a lineage which recently spread globally

Pertussis is a highly contagious, acute respiratory disease in humans caused by the Gram-negative pathogen Bordetella pertussis. Pertussis has resurged in the face of intensive vaccination and this has coincided with the emergence of strains carrying a particular allele for the pertussis toxin promoter, ptxP3, which is associated with higher levels of pertussis toxin (Ptx) production. Within 10 to 20 years, ptxP3 strains have nearly completely replaced the previously dominant ptxP1 strains resulting in a worldwide selective sweep. In order to identify B. pertussis genes associated with the selective sweep, we compared the expression of genes in ptxP1 and ptxP3 strains that are under control of the Bordetella master virulence regulatory locus (bvgASR). The BvgAS proteins comprise a two component sensory transduction system which is regulated by temperature, nicotinic acid and sulfate. By increasing the sulfate concentration, it is possible to change the phase of B. pertussis from virulent to avirulent. Until recently, the only distinctive phenotype of ptxP3 strains was a higher Ptx production. Here we identify additional phenotypic differences between ptxP1 and ptxP3 strains which may have contributed to its global spread by comparing global transcriptional responses under sulfate-modulating conditions. We show that ptxP3 strains are less sensitive to sulfate-mediated gene suppression, resulting in an increased production of the vaccine antigens pertactin (Prn) and Ptx and a number of other virulence genes, including a type III secretion toxin, Vag8, a protein involved in complement resistance, and lpxE involved in lipid A modification. Furthermore, enhanced expression of the vaccine antigens Ptx and Prn by ptxP3 strains was confirmed at the protein level. Identification of genes differentially expressed between ptxP1 and ptxP3 strains may elucidate how B. pertussis has adapted to vaccination and allow the improvement of pertussis vaccines by identifying novel vaccine candidates.

Pertussis in the Netherlands, is the current vaccination strategy sufficient to reduce disease burden in young infants?

BACKGROUND

Pertussis has resurged in the Netherlands since 1996. Several measures, i.e. acceleration of the schedule, introduction of a preschool acellular pertussis booster and change from an infant whole cell to an acellular pertussis combination vaccine were implemented in the National Immunisation Programme to decrease disease burden, in particular among very young infants who have the highest morbidity and mortality of pertussis. Nevertheless, a large outbreak occurred in 2011-2012.

METHODS

1996-2010 was divided in 3-year-periods to assess the impact of the measures taken, using notifications and hospitalisations. These results were compared with 2011-2012. Mean Incidence rates (IRs) per 100,000 were calculated.

RESULTS

Although the measures taken resulted in decreased IRs among the targeted age groups after implementation, overall mean IRs of notifications increased from 32 (1996-2004) to 37 (2005-2010) and 63 (2011-2012). Young infants, not yet vaccinated, did not benefit; during the 2011-2012 outbreak, IR in 0-2-month-olds amounted to 259.6. IR among persons over 9 years of age increased from 6.8 (1996-1999) to 59.1 (2011-2012) For hospitalisations overall mean IRs decreased from 1.95 per 100,000 (1997-2004) to 0.88 (2005-2010) and 0.76 (2011).

CONCLUSION

The measures taken reduced IRs of notifications and hospitalisations among groups eligible for vaccination, but had no effect on the increasing IRs in adolescents and adults. This trend is also observed in other countries. The high IRs in 2012 in adolescents and adults probably resulted in increased transmission to infants, who are at risk for contracting severe pertussis. Therefore, additional measures to protect this group should be considered.

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.

Pertussis resurgence: waning immunity and pathogen adaptation - two sides of the same coin

Pertussis or whooping cough has persisted and resurged in the face of vaccination and has become one of the most prevalent vaccine-preventable diseases in Western countries. The high circulation rate of Bordetella pertussis poses a threat to infants that have not been (completely) vaccinated and for whom pertussis is a severe, life-threatening, disease. The increase in pertussis is mainly found in age groups in which immunity has waned and this has resulted in the perception that waning immunity is the main or exclusive cause for the resurgence of pertussis. However, significant changes in B. pertussis populations have been observed after the introduction of vaccinations, suggesting a role for pathogen adaptation in the persistence and resurgence of pertussis. These changes include antigenic divergence with vaccine strains and increased production of pertussis toxin. Antigenic divergence will affect both memory recall and the efficacy of antibodies, while higher levels of pertussis toxin may increase suppression of the innate and acquired immune system. We propose these adaptations of B. pertussis have decreased the period in which pertussis vaccines are effective and thus enhanced the waning of immunity. We plead for a more integrated approach to the pertussis problem which includes the characteristics of the vaccines, the B. pertussis populations and the interaction between the two.

Development of a temperature-switch PCR-based SNP typing method for Mycobacterium ulcerans

Mycobacterium ulcerans (M. ulcerans), the causative agent of the devastating skin disease Buruli ulcer (BU), is characterized by an extremely low level of genetic diversity. Recently, we have reported the first discrimination of closely related M. ulcerans variants in the BU endemic Densu River Valley of Ghana. In the study real-time PCR-based single nucleotide polymorphism (SNP) typing at 89 predefined loci revealed the presence of ten M. ulcerans haplotypes circulating in the BU endemic region. Here we describe the development of temperature-switch PCR (TSP) assays that allow distinguishing these haplotypes by conventional agarose gel-based analysis of the PCR products. After validation of the accuracy of typing results, the TSP assays were successfully established in a reference laboratory in Ghana. Development of the cost-effective and rapid TSP-based genetic fingerprinting method will thus allow investigating the spread of M. ulcerans clones by regular genetic monitoring in BU endemic countries.

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.

Direct molecular typing of Bordetella pertussis from clinical specimens submitted for diagnostic quantitative (real-time) PCR

Molecular typing of Bordetella pertussis is routinely performed on bacterial isolates, but not on DNA extracted from nasopharyngeal aspirates or pernasal swabs submitted for diagnostic real-time PCR (qPCR). We investigated whether these DNA extracts were suitable for multilocus variable-number tandem repeat analysis (MLVA) and DNA sequence-based typing. We analysed all the available qPCR-positive samples received by our laboratory from patients <1 year of age between January 2008 and August 2010. Eighty-one per cent (106/131) of these generated a complete MLVA profile. This rose to 92 % (105/114) if only samples positive for both of the two targets used for the B. pertussis PCR (insertion element IS481 and pertussis toxin promoter ptxP) were analysed. Sequence-based typing of the pertactin, pertussis toxin S1 subunit and pertussis promoter regions (prn, ptxA and ptxP) was attempted on 89 of the DNA extracts that had generated a full MLVA profile. Eighty-three (93 %) of these produced complete sequences for all three targets. Comparison of molecular typing data from the 89 extracts with those from 111 contemporary bacterial isolates showed that the two sources yielded the same picture of the B. pertussis population [dominated by the MLVA-27 prn(2) ptxA(1) ptxP(3) clonal type]. There was no significant difference in MLVA type distribution or diversity between the two sample sets. This suggests that clinical extracts can be used in place of, or to complement, bacterial cultures for typing purposes (at least, in this age group). With small modifications to methodology, generating MLVA and sequence-based typing data from qPCR-positive clinical DNA extracts is likely to generate a complete dataset in the majority of samples from the <1 year age group. Its success with samples from older subjects remains to be seen. However, our data suggest that it is suitable for inclusion in molecular epidemiological studies of the B. pertussis population or as a tool in outbreak investigations.

Molecular methods for pathogen and microbial community detection and characterization: current and potential application in diagnostic microbiology

Clinical microbiology laboratories worldwide have historically relied on phenotypic methods (i.e., culture and biochemical tests) for detection, identification and characterization of virulence traits (e.g., antibiotic resistance genes, toxins) of human pathogens. However, limitations to implementation of molecular methods for human infectious diseases testing are being rapidly overcome allowing for the clinical evaluation and implementation of diverse technologies with expanding diagnostic capabilities. The advantages and limitation of molecular techniques including real-time polymerase chain reaction, partial or whole genome sequencing, molecular typing, microarrays, broad-range PCR and multiplexing will be discussed. Finally, terminal restriction fragment length polymorphism (T-RFLP) and deep sequencing are introduced as technologies at the clinical interface with the potential to dramatically enhance our ability to diagnose infectious diseases and better define the epidemiology and microbial ecology of a wide range of complex infections.

Selection and emergence of pertussis toxin promoter ptxP3 allele in the evolution of Bordetella pertussis

Evolutionary studies using single nucleotide polymorphisms (SNPs) have separated Bordetella pertussis isolates into six major clusters, with recent isolates forming cluster I. The expansion of cluster I isolates was characterised by changes in genes encoding antigenic components in acellular vaccines, including pertactin (Prn). Here, we determined the initial emergence of the pertussis toxin promoter allele, ptxP3, from an evolutionary perspective. This allele was previously shown in a study from the Netherlands to be associated with increased pertussis toxin production as a result of a single base mutation in the ptxP. The ptxP region of 313 worldwide isolates was sequenced, including 208 isolates from Australia collected over a 40 year period. Eight alleles were identified, of which only two predominated: ptxP1 and ptxP3. One novel allele was also found. ptxP3 was only found in SNP cluster I of B. pertussis and its emergence is concurrent with the change to the non-vaccine prn2 allele. Our results suggest that the globally distributed cluster I of B. pertussis has the ability to evade vaccine induced selection pressure.

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.