Machine learning to attribute the source of Campylobacter infections in the United States: A retrospective analysis of national surveillance data.

OBJECTIVES: Integrating pathogen genomic surveillance with bioinformatics can enhance public health responses by identifying risk and guiding interventions. This study focusses on the two predominant Campylobacter species, which are commonly found in the gut of birds and mammals and often infect humans via contaminated food. Rising incidence and antimicrobial resistance (AMR) are a global concern, and there is an urgent need to quantify the main routes to human infection. METHODS: During routine US national surveillance (2009-2019), 8856 Campylobacter genomes from human infections and 16,703 from possible sources were sequenced. Using machine learning and probabilistic models, we target genetic variation associated with host adaptation to attribute the source of human infections and estimate the importance of different disease reservoirs. RESULTS: Poultry was identified as the primary source of human infections, responsible for an estimated 68% of cases, followed by cattle (28%), and only a small contribution from wild birds (3%) and pork sources (1%). There was also evidence of an increase in multidrug resistance, particularly among isolates attributed to chickens. CONCLUSIONS: National surveillance and source attribution can guide policy, and our study suggests that interventions targeting poultry will yield the greatest reductions in campylobacteriosis and spread of AMR in the US. DATA AVAILABILITY: All sequence reads were uploaded and shared on NCBI's Sequence Read Archive (SRA) associated with BioProjects; PRJNA239251 (CDC / PulseNet surveillance), PRJNA287430 (FSIS surveillance), PRJNA292668 & PRJNA292664 (NARMS) and PRJNA258022 (FDA surveillance). Publicly available genomes, including reference genomes and isolates sampled worldwide from wild birds are associated with BioProject accessions: PRJNA176480, PRJNA177352, PRJNA342755, PRJNA345429, PRJNA312235, PRJNA415188, PRJNA524300, PRJNA528879, PRJNA529798, PRJNA575343, PRJNA524315 and PRJNA689604. Contiguous assemblies of all genome sequences compared are available at Mendeley data (assembled C. coli genomes doi: 10.17632/gxswjvxyh3.1; assembled C. jejuni genomes doi: 10.17632/6ngsz3dtbd.1) and individual project and accession numbers can be found in Supplementary tables S1 and S2, which also includes pubMLST identifiers for assembled genomes. Figshare (10.6084/m9.figshare.20279928). Interactive phylogenies are hosted on microreact separately for C. jejuni (https://microreact.org/project/pascoe-us-cjejuni) and C. coli (https://microreact.org/project/pascoe-us-ccoli).

Development of the Pneumococcal Genome Library, a core genome multilocus sequence typing scheme, and a taxonomic life identification number barcoding system to investigate and define pneumococcal population structure.

Investigating the genomic epidemiology of major bacterial pathogens is integral to understanding transmission, evolution, colonization, disease, antimicrobial resistance and vaccine impact. Furthermore, the recent accumulation of large numbers of whole genome sequences for many bacterial species enhances the development of robust genome-wide typing schemes to define the overall bacterial population structure and lineages within it. Using the previously published data, we developed the Pneumococcal Genome Library (PGL), a curated dataset of 30 976 genomes and contextual data for carriage and disease pneumococci recovered between 1916 and 2018 in 82 countries. We leveraged the size and diversity of the PGL to develop a core genome multilocus sequence typing (cgMLST) scheme comprised of 1222 loci. Finally, using multilevel single-linkage clustering, we stratified pneumococci into hierarchical clusters based on allelic similarity thresholds and defined these with a taxonomic life identification number (LIN) barcoding system. The PGL, cgMLST scheme and LIN barcodes represent a high-quality genomic resource and fine-scale clustering approaches for the analysis of pneumococcal populations, which support the genomic epidemiology and surveillance of this leading global pathogen.

Profiling IgG and IgA antibody responses during vaccination and infection in a high-risk gonorrhoea population.

Development of a vaccine against gonorrhoea is a global priority, driven by the rise in antibiotic resistance. Although Neisseria gonorrhoeae (Ng) infection does not induce substantial protective immunity, highly exposed individuals may develop immunity against re-infection with the same strain. Retrospective epidemiological studies have shown that vaccines containing Neisseria meningitidis (Nm) outer membrane vesicles (OMVs) provide a degree of cross-protection against Ng infection. We conducted a clinical trial (NCT04297436) of 4CMenB (Bexsero, GSK), a licensed Nm vaccine containing OMVs and recombinant antigens, comprising a single arm, open label study of two doses with 50 adults in coastal Kenya who have high exposure to Ng. Data from a Ng antigen microarray established that serum IgG and IgA reactivities against the gonococcal homologs of the recombinant antigens in the vaccine peaked at 10 but had declined by 24 weeks. For most reactive OMV-derived antigens, the reverse was the case. A cohort of similar individuals with laboratory-confirmed gonococcal infection were compared before, during, and after infection: their reactivities were weaker and differed from the vaccinated cohort. We conclude that the cross-protection of the 4CMenB vaccine against gonorrhoea could be explained by cross-reaction against a diverse selection of antigens derived from the OMV component.

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae.

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria.

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

The genomic epidemiology of Neisseria meningitidis carriage from a randomised controlled trial of 4CMenB vaccination in an asymptomatic adolescent population.

Background: Oropharyngeal carriage of Neisseria meningitidis is frequent during adolescence, representing a major source of invasive meningococcal disease. This study examined the impact of a serogroup B vaccination (Bexsero, GSK 4CMenB) programme on adolescent N. meningitidis carriage using genomic data. Methods: A total 34,489 oropharyngeal samples were collected as part of a state-wide cluster randomised-controlled trial in South Australia during 2017 and 2018 (NCT03089086). Samples were screened for the presence of N. meningitidis DNA by porA PCR prior to culture. Whole genome sequencing was performed on all 1772 N. meningitidis culture isolates and their genomes were analysed. Findings: Unencapsulated meningococci were predominant at baseline (36.3% of isolates), followed by MenB (31.0%), and MenY (20.5%). Most MenB were ST-6058 from hyperinvasive cc41/44, or ST-32 and ST-2870 from cc32. For MenY, ST-23 and ST-1655 from cc23 were prevalent. Meningococcal carriage was mostly unchanged due to the vaccination programme; however, a significant reduction in ST-53 capsule-null meningococci prevalence was observed in 2018 compared to 2017 (OR = 0.52; 95% CI: 0.30-0.87, p = 0.0106). This effect was larger in the vaccinated compared to the control group (OR = 0.37; 95% CI: 0.12-0.98, p = 0.0368). Interpretation: While deployment of the 4CMenB vaccination did not alter the carriage of hyperinvasive MenB in the vaccinated population, it altered the carriage of other N. meningitidis sequence types following the vaccination program. Our findings suggest 4CMenB vaccination is unlikely to reduce transmission of hyperinvasive N. meningitidis strains and therefore ongoing targeted vaccination is likely a more effective public health intervention. Funding: This work was funded by GlaxoSmithKline Biologicals SA.

Longitudinal study of meningococcal carriage in adolescents and young adults in South Australia 2017-2020.

BACKGROUND: This analysis investigated longitudinal changes in meningococcal carriage in adolescents in South Australia over 4 years. METHODS: Data from the "B Part of It" study, which included a state-wide cluster randomized controlled trial in secondary-school students (n = 34,489 in 2017 and 2018) and serial cross-sectional studies in school leavers aged 17-25 years (n = 4028 in 2019-2020). Individuals had oropharyngeal swabs collected annually. This study included two unique cohorts: (1) individuals enrolled in 2019, with three consecutive annual swabs taken in 2017, 2018 and 2019; and (2) individuals enrolled in 2020, with swabs taken in 2017, 2018, and 2020. Disease-associated N. meningitidis genogroups were identified using PCR and whole genome sequencing. Univariate analysis identified risk factors for recurrent carriage (≥2). RESULTS: Among school leavers, 50 (1.7%, total n = 2980) had carriage detected at successive visits. In participants with meningococcal carriage at successive visits, 38/50 (76.0%) had the same genogroup detected by porA PCR. Of those, 19 had the same MLST type and demonstrated minimal variation, indicating they most likely had sustained carriage of the same isolate (range 226 to 490 days, mean duration 352 [SD 51] days). In the 2019 school leaver cohort, 6.7% acquired carriage in their first year out of school compared to 3.3% in their final school year. Compared to single carriage detection, recurrent carriage was potentially more likely in older adolescents (16 compared to ≤15 years; OR = 1.97 (95%CI 1.0, 3.86); p = 0.048). CONCLUSION: Whilst carriage is typically transient, some adolescents/young adults may have persistent carriage and are likely to be an important group in the transmission of meningococci.

The metabolic, virulence and antimicrobial resistance profiles of colonising Streptococcus pneumoniae shift after PCV13 introduction in urban Malawi.

Streptococcus pneumoniae causes substantial mortality among children under 5-years-old worldwide. Polysaccharide conjugate vaccines (PCVs) are highly effective at reducing vaccine serotype disease, but emergence of non-vaccine serotypes and persistent nasopharyngeal carriage threaten this success. We investigated the hypothesis that following vaccine, adapted pneumococcal genotypes emerge with the potential for vaccine escape. We genome sequenced 2804 penumococcal isolates, collected 4-8 years after introduction of PCV13 in Blantyre, Malawi. We developed a pipeline to cluster the pneumococcal population based on metabolic core genes into "Metabolic genotypes" (MTs). We show that S. pneumoniae population genetics are characterised by emergence of MTs with distinct virulence and antimicrobial resistance (AMR) profiles. Preliminary in vitro and murine experiments revealed that representative isolates from emerging MTs differed in growth, haemolytic, epithelial infection, and murine colonisation characteristics. Our results suggest that in the context of PCV13 introduction, pneumococcal population dynamics had shifted, a phenomenon that could further undermine vaccine control and promote spread of AMR.

Evolution and exchange of plasmids in pathogenic Neisseria.

IMPORTANCE: Horizontal gene transfer (HGT) is a major influence in driving the spread of antimicrobial resistance (AMR) in many bacteria. A conjugative plasmid which is widespread in Neisseria gonorrhoeae, pConj, prevented the use of tetracycline/doxycycline for treating gonococcal infection. Here, we show that pConj evolved in the related pathogen, Neisseria meningitidis, and has been acquired by the gonococcus from the meningococcus on multiple occasions. Following its initial acquisition, pConj spread to different gonococcal lineages; changes in the plasmid's conjugation machinery associated with another transfer event limit spread in the gonococcal populations. Our findings have important implications for the use of doxycycline to prevent bacterial sexually transmitted disease which is likely to exacerbate the spread of AMR through HGT in pathogenic bacteria.

Gut microbiota of the critically endangered Saiga antelope across two wild populations in a year without mass mortality.

The Saiga are migratory antelopes inhabiting the grasslands of Eurasia. Over the last century, Saiga have been pushed to the brink of extinction by mass mortality events and intense poaching. Yet, despite the high profile of the Saiga as an animal of conservation concern, little is known of its biology. In particular, the gut microbiota of Saiga has not been studied, despite its potential importance in health. Here, we characterise the gut microbiota of Saiga from two geographically distinct populations in Kazakhstan and compare it with that of other antelope species. We identified a consistent gut microbial diversity and composition among individuals and across two Saiga populations during a year without die-offs, with over 85% of bacterial genera being common to both populations despite vast geographic separation. We further show that the Saiga gut microbiota resembled that of five other antelopes. The putative causative agent of Saiga mass die-offs, Pasteurella multocida, was not detected in the Saiga microbiota. Our findings provide the first description of the Saiga gut microbiota, generating a baseline for future work investigating the microbiota's role in health and mass die-offs, and supporting the conservation of this critically endangered species.

Acinetobacter baumannii IC2 and IC5 Isolates with Co-Existing blaOXA-143-like and blaOXA-72 and Exhibiting Strong Biofilm Formation in a Mexican Hospital.

Acinetobacter baumannii is an opportunistic pathogen responsible for healthcare-associated infections (HAIs) and outbreaks. Antimicrobial resistance mechanisms and virulence factors allow it to survive and spread in the hospital environment. However, the molecular mechanisms of these traits and their association with international clones are frequently unknown in low- and middle-income countries. Here, we analyze the phenotype and genotype of seventy-six HAIs and outbreak-causing A. baumannii isolates from a Mexican hospital over ten years, with special attention to the carbapenem resistome and biofilm formation. The isolates belonged to the global international clone (IC) 2 and the Latin America endemic IC5 and were predominantly extensively drug-resistant (XDR). Oxacillinases were identified as a common source of carbapenem resistance. We noted the presence of the blaOXA-143-like family (not previously described in Mexico), the blaOXA-72 and the blaOXA-398 found in both ICs. A low prevalence of efflux pump overexpression activity associated with carbapenem resistance was observed. Finally, strong biofilm formation was found, and significant biofilm-related genes were identified, including bfmRS, csuA/BABCDE, pgaABCD and ompA. This study provides a comprehensive profile of the carbapenem resistome of A. baumannii isolates belonging to the same pulse type, along with their significant biofilm formation capacity. Furthermore, it contributes to a better understanding of their role in the recurrence of infection and the endemicity of these isolates in a Mexican hospital.

There are three major Neisseria gonorrhoeae ß-lactamase plasmid variants which are associated with specific lineages and carry distinct TEM alleles.

Neisseria gonorrhoeae is a significant threat to global health with an estimated incidence of over 80 million cases each year and high levels of antimicrobial resistance. The gonococcal ß-lactamase plasmid, pbla, carries the TEM ß-lactamase, which requires only one or two amino acid changes to become an extended-spectrum ß-lactamase (ESBL); this would render last resort treatments for gonorrhoea ineffective. Although pbla is not mobile, it can be transferred by the conjugative plasmid, pConj, found in N. gonorrhoeae. Seven variants of pbla have been described previously, but little is known about their frequency or distribution in the gonococcal population. We characterised sequences of pbla variants and devised a typing scheme, Ng_pblaST that allows their identification from whole genome short-read sequences. We implemented Ng_pblaST to assess the distribution of pbla variants in 15 532 gonococcal isolates. This demonstrated that only three pbla variants commonly circulate in gonococci, which together account for >99 % of sequences. The pbla variants carry different TEM alleles and are prevalent in distinct gonococcal lineages. Analysis of 2758 pbla-containing isolates revealed the co-occurrence of pbla with certain pConj types, indicating co-operativity between pbla and pConj variants in the spread of plasmid-mediated AMR in N. gonorrhoeae. Understanding the variation and distribution of pbla is essential for monitoring and predicting the spread of plasmid-mediated ß-lactam resistance in N. gonorrhoeae.

Genome-Wide Association Studies Identify an Association of Transferrin Binding Protein B Variation and Invasive Serogroup Y Meningococcal Disease in Older Adults.

BACKGROUND: Neisseria meningitidis serogroup Y, especially ST-23 clonal complex (Y:cc23), represents a larger proportion of invasive meningococcal disease (IMD) in older adults compared to younger individuals. This study explored the meningococcal genetic variation underlying this association. METHODS: Maximum-likelihood phylogenies and the pangenome were analyzed using whole-genome sequence (WGS) data from 200 Y:cc23 isolates in the Neisseria PubMLST database. Genome-wide association studies (GWAS) were performed on WGS data from 250 Y:cc23 isolates from individuals with IMD aged ≥65 years versus < 65 years. RESULTS: Y:cc23 meningococcal variants did not cluster by age group or disease phenotype in phylogenetic analyses. Pangenome comparisons found no differences in presence or absence of genes in IMD isolates from the different age groups. GWAS identified differences in nucleotide polymorphisms within the transferrin-binding protein B (tbpB) gene in isolates from individuals ≥65 years of age. TbpB structure modelling suggests these may impact binding of human transferrin. CONCLUSIONS: These data suggest differential iron scavenging capacity amongst Y:cc23 meningococci isolated from older compared to younger patients. Iron acquisition is essential for many bacterial pathogens including the meningococcus. These polymorphisms may facilitate colonization, thereby increasing the risk of disease in vulnerable older people with altered nasopharyngeal microbiomes and nutritional status.

Ribosomal MLST nucleotide identity (rMLST-NI), a rapid bacterial species identification method: application to Klebsiella and Raoultella genomic species validation.

Bacterial genomics is making an increasing contribution to the fields of medicine and public health microbiology. Consequently, accurate species identification of bacterial genomes is an important task, particularly as the number of genomes stored in online databases increases rapidly and new species are frequently discovered. Existing database entries require regular re-evaluation to ensure that species annotations are consistent with the latest species definitions. We have developed an automated method for bacterial species identification that is an extension of ribosomal multilocus sequence typing (rMLST). The method calculates an 'rMLST nucleotide identity' (rMLST-NI) based on the nucleotides present in the protein-encoding ribosomal genes derived from bacterial genomes. rMLST-NI was used to validate the species annotations of 11839 publicly available Klebsiella and Raoultella genomes based on a comparison with a library of type strain genomes. rMLST-NI was compared with two whole-genome average nucleotide identity methods (OrthoANIu and FastANI) and the k-mer based Kleborate software. The results of the four methods agreed across a dataset of 11839 bacterial genomes and identified a small number of entries (n=89) with species annotations that required updating. The rMLST-NI method was 3.5 times faster than Kleborate, 4.5 times faster than FastANI and 1600 times faster than OrthoANIu. rMLST-NI represents a fast and generic method for species identification using type strains as a reference.

Emerging patterns of fluoroquinolone resistance in Campylobacter jejuni in the UK [1998-2018].

Campylobacter jejuni (C.jejuni) is the most common causative agent of bacterial food poisoning worldwide and is known to be genetically highly diverse. C. jejuni is increasingly resistant to fluoroquinolone antibiotics, but very few studies have investigated variant-specific patterns of resistance across time. Here we use statistical modelling and clustering techniques to investigate patterns of fluoroquinolone resistance amongst 10,359 UK isolates from human disease sampled over 20 years. We observed six distinct patterns of fluoroquinolone sensitivity/resistance in C. jejuni across time, grouping by clonal complex (CC). Some CCs were fully resistant, some shifted from susceptible to resistant following a sigmoidal shape, and some remained susceptible over time. Our findings indicate that the fluoroquinolone resistance patterns of C. jejuni are complicated and cannot be analysed as a single species but divided into variant dynamics so that the factors driving resistance can be thoroughly investigated.

Impact of meningococcal ACWY conjugate vaccines on pharyngeal carriage in adolescents: evidence for herd protection from the UK MenACWY programme.

OBJECTIVE: Serogroup W and Y invasive meningococcal disease increased globally from 2000 onwards. Responding to a rapid increase in serogroup W clonal complex 11 (W:cc11) invasive meningococcal disease, the UK replaced an adolescent booster dose of meningococcal C conjugate vaccine with quadrivalent MenACWY conjugate vaccine in 2015. By 2018, the vaccine coverage in the eligible school cohorts aged 14 to 19 years was 84%. We assessed the impact of the MenACWY vaccination programme on meningococcal carriage. METHODS: An observational study of culture-defined oropharyngeal meningococcal carriage prevalence before and after the start of the MenACWY vaccination programme in UK school students, aged 15 to 19 years, using two cross-sectional studies: 2014 to 2015 "UKMenCar4" and 2018 "Be on the TEAM" (ISRCTN75858406). RESULTS: A total of 10 625 participants preimplementation and 13 438 postimplementation were included. Carriage of genogroups C, W, and Y (combined) decreased from 2.03 to 0.71% (OR 0.34 [95% CI 0.27-0.44], p < 0.001). Carriage of genogroup B meningococci did not change (1.26% vs 1.23% [95% CI 0.77-1.22], p = 0.80) and genogroup C remained rare (n = 7/10 625 vs 17/13 438, p = 0.135). The proportion of serogroup positive isolates (i.e. those expressing capsule) decreased for genogroup W by 53.8% (95% CI -5.0 - 79.8, p = 0.016) and for genogroup Y by 30.1% (95% CI 8.946·3, p = 0.0025). DISCUSSION: The UK MenACWY vaccination programme reduced carriage acquisition of genogroup and serogroup Y and W meningococci and sustained low levels of genogroup C carriage. These data support the use of quadrivalent MenACWY conjugate vaccine for indirect (herd) protection.

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.

A Dual Barcoding Approach to Bacterial Strain Nomenclature: Genomic Taxonomy of Klebsiella pneumoniae Strains.

Sublineages (SLs) within microbial species can differ widely in their ecology and pathogenicity, and their precise definition is important in basic research and for industrial or public health applications. Widely accepted strategies to define SLs are currently missing, which confuses communication in population biology and epidemiological surveillance. Here, we propose a broadly applicable genomic classification and nomenclature approach for bacterial strains, using the prominent public health threat Klebsiella pneumoniae as a model. Based on a 629-gene core genome multilocus sequence typing (cgMLST) scheme, we devised a dual barcoding system that combines multilevel single linkage (MLSL) clustering and life identification numbers (LINs). Phylogenetic and clustering analyses of >7,000 genome sequences captured population structure discontinuities, which were used to guide the definition of 10 infraspecific genetic dissimilarity thresholds. The widely used 7-gene multilocus sequence typing (MLST) nomenclature was mapped onto MLSL SLs (threshold: 190 allelic mismatches) and clonal group (threshold: 43) identifiers for backwards nomenclature compatibility. The taxonomy is publicly accessible through a community-curated platform (https://bigsdb.pasteur.fr/klebsiella), which also enables external users' genomic sequences identification. The proposed strain taxonomy combines two phylogenetically informative barcode systems that provide full stability (LIN codes) and nomenclatural continuity with previous nomenclature (MLSL). This species-specific dual barcoding strategy for the genomic taxonomy of microbial strains is broadly applicable and should contribute to unify global and cross-sector collaborative knowledge on the emergence and microevolution of bacterial pathogens.

A restatement of the natural science evidence base regarding the source, spread and control of Campylobacter species causing human disease.

Food poisoning caused by Campylobacter (campylobacteriosis) is the most prevalent bacterial disease associated with the consumption of poultry, beef, lamb and pork meat and unpasteurized dairy products. A variety of livestock industry, food chain and public health interventions have been implemented or proposed to reduce disease prevalence, some of which entail costs for producers and retailers. This paper describes a project that set out to summarize the natural science evidence base relevant to campylobacteriosis control in as policy-neutral terms as possible. A series of evidence statements are listed and categorized according to the nature of the underlying information. The evidence summary forms the appendix to this paper and an annotated bibliography is provided in the electronic supplementary material.

Application of a Neisseria meningitidis antigen microarray to identify candidate vaccine proteins from a human Phase I clinical trial.

Meningococcal meningitis is a rare but serious condition affecting mainly children and young adults. Outer membrane vesicles (OMV) from Neisseria meningitidis have been used successfully as vaccines against the disease, although they only provide protection against a limited number of the many existing variants. There have been many attempts to identify suitable protein antigens for use in defined vaccines that provide broad protection against the disease, such as that leading to the development of the four component 4CMenB vaccine. We previously reported the use of a protein antigen microarray to screen for IgG antibodies in sera derived from human recipients of an OMV-based vaccine, as part of a Phase I clinical trial. Here, we show that computational methods can be used to cluster antigens that elicit similar responses in the same individuals. Fitting of IgG antibody binding data to 4,005 linear regressions identified pairs of antigens that exhibited significant correlations. Some were from the same antigens in different quaternary states, whilst others might be correlated for functional or immunological reasons. We also conducted statistical analyses to examine correlations between individual serum bactericidal antibody (SBA) titres and IgG reactivity against specific antigens. Both Kendall's tau and Spearman's rank correlation coefficient statistics identified specific antigens that correlated with log(SBA) titre in five different isolates. The principal antigens identified were PorA and PorB, RmpM, OpcA, and the type IV pilus assembly secretin, PilQ. Other minor antigens identified included a lipoprotein, two proteins from the BAM complex and the efflux channel MtrE. Our results suggest that consideration of the entire antigen composition, and allowance for potential interaction between antigens, could be valuable in designing future meningococcal vaccines. Such an approach has the advantages that it uses data derived from human, rather than animal, immunization and that it avoids the need to screen individual antigens.