Search

The Bengal Bay clone (ST772) is a community-associated and multidrug-resistant Staphylococcus aureus lineage first isolated from Bangladesh and India in 2004. In this study, we showed that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally, resulting in small-scale community and health care outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug resistance of health care-associated S. aureus lineages. This study demonstrates the importance of whole-genome sequencing for the surveillance of highly antibiotic-resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world., The evolution and global transmission of antimicrobial resistance have been well documented for Gram-negative bacteria and health care-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. In this study, we traced the recent origins and global spread of a multidrug-resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole-genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data show that the clone emerged on the Indian subcontinent in the early 1960s and disseminated rapidly in the 1990s. Short-term outbreaks in community and health care settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the emergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth, and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional health care-associated clones with the epidemiological transmission of community-associated methicillin-resistant S. aureus (MRSA). Our study demonstrates the importance of whole-genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.

Restriction-Modification systems (RMS) are one of the main mechanisms of defence against foreign DNA invasion and can have an important role in the regulation of gene expression. The obligate human pathogen Neisseria gonorrhoeae carries one of the highest loads of RMS in its genome; between 13 to 15 of the three main types. Previous work has described their organization in the reference genome FA1090 and has experimentally inferred the associated methylated motifs. Here, we studied the structure of RMS and target methylated motifs in 25 gonococcal strains sequenced with Single Molecule Real-Time (SMRT) technology, which provides data on DNA modification. The results showed a variable picture of active RMS in different strains, with phase variation switching the activity of Type III RMS, and both the activity and specificity of a Type I RMS. Interestingly, the Dam methylase was found in place of the NgoAXI endonuclease in two of the strains, despite being previously thought to be absent in the gonococcus. We also identified the real methylation target of NgoAX as 5’-GCAGA-3’, different from that previously described. Results from this work give further insights into the diversity and dynamics of RMS and methylation patterns in N. gonorrhoeae.

Recent advances in the scale and diversity of population genomic datasets for bacteria now provide the potential for genome-wide patterns of co-evolution to be studied at the resolution of individual bases. Here we describe a new statistical method, genomeDCA, which uses recent advances in computational structural biology to identify the polymorphic loci under the strongest co-evolutionary pressures. We apply genomeDCA to two large population data sets representing the major human pathogens Streptococcus pneumoniae (pneumococcus) and Streptococcus pyogenes (group A Streptococcus). For pneumococcus we identified 5,199 putative epistatic interactions between 1,936 sites. Over three-quarters of the links were between sites within the pbp2x, pbp1a and pbp2b genes, the sequences of which are critical in determining non-susceptibility to beta-lactam antibiotics. A network-based analysis found these genes were also coupled to that encoding dihydrofolate reductase, changes to which underlie trimethoprim resistance. Distinct from these antibiotic resistance genes, a large network component of 384 protein coding sequences encompassed many genes critical in basic cellular functions, while another distinct component included genes associated with virulence. The group A Streptococcus (GAS) data set population represents a clonal population with relatively little genetic variation and a high level of linkage disequilibrium across the genome. Despite this, we were able to pinpoint two RNA pseudouridine synthases, which were each strongly linked to a separate set of loci across the chromosome, representing biologically plausible targets of co-selection. The population genomic analysis method applied here identifies statistically significantly co-evolving locus pairs, potentially arising from fitness selection interdependence reflecting underlying protein-protein interactions, or genes whose product activities contribute to the same phenotype. This discovery approach greatly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work., Author summary Epistatic interactions between polymorphisms in DNA are recognized as important drivers of evolution in numerous organisms. Study of epistasis in bacteria has been hampered by the lack of densely sampled population genomic data, suitable statistical models and inference algorithms sufficiently powered for extremely high-dimensional parameter spaces. We introduce the first model-based method for genome-wide epistasis analysis and use two of the largest available bacterial population genome data sets on Streptococcus pneumoniae (the pneumococcus) and Streptococcus pyogenes (group A Streptococcus) to demonstrate its potential for biological discovery. Our approach reveals interacting networks of resistance, virulence and core machinery genes in the pneumococcus, which highlights putative candidates for novel drug targets. We also discover a number of plausible targets of co-selection in S. pyogenes linked to RNA pseudouridine synthases. Our method significantly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work.

Neisseria meningitidis is the leading cause of epidemic meningitis in the “meningitis belt” of Africa, where clonal waves of colonization and disease are observed. Point mutations and horizontal gene exchange lead to constant diversification of meningococcal populations during clonal spread. Maintaining a high genomic diversity may be an evolutionary strategy of meningococci that increases chances of fixing occasionally new highly successful “fit genotypes”. We have performed a longitudinal study of meningococcal carriage and disease in northern Ghana by analyzing cerebrospinal fluid samples from all suspected meningitis cases and monitoring carriage of meningococci by twice yearly colonization surveys. In the framework of this study, we observed complete replacement of an A: sequence types (ST)-2859 clone by a W: ST-2881 clone. However, after a gap of 1 year, A: ST-2859 meningococci re-emerged both as colonizer and meningitis causing agent. Our whole genome sequencing analyses compared the A population isolated prior to the W colonization and disease wave with the re-emerging A meningococci. This analysis revealed expansion of one clone differing in only one nonsynonymous SNP from several isolates already present in the original A: ST-2859 population. The colonization bottleneck caused by the competing W meningococci thus resulted in a profound reduction in genomic diversity of the A meningococcal population.

The sexually transmitted pathogen Neisseria gonorrhoeae is regarded as being on the way to becoming an untreatable superbug. Despite its clinical importance, little is known about its emergence and evolution, and how this corresponds with the introduction of antimicrobials. We present a genome-based phylogeographical analysis of 419 gonococcal isolates from across the globe. Results indicate that modern gonococci originated in Europe or Africa, possibly as late as the sixteenth century and subsequently disseminated globally. We provide evidence that the modern gonococcal population has been shaped by antimicrobial treatment of sexually transmitted infections as well as other infections, leading to the emergence of two major lineages with different evolutionary strategies. The well-described multidrug-resistant lineage is associated with high rates of homologous recombination and infection in high-risk sexual networks. A second, multisusceptible lineage is more associated with heterosexual networks, with potential implications for infection control.

14 páginas, 6 figuras, Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants., This work was funded by projects of the European Research Council (638553-TB-ACCELERATE) and Ministerio de Economía y Competitividad (Spanish government) research grant SAF2016-77346-R (to I.C.); BFU2014-58656-R and BFU2017-89594-R from Ministerio de Economía y Competitividad (Spanish government) and PROMETEO/2016/122 from Generalitat Valenciana (to F.G.-C); and the Swiss National Science Foundation (grants 310030_166687, IZRJZ3_164171, IZLSZ3_170834, and CRSII5_177163), the European Research Council (309540-EVODRTB), and SystemsX.ch (to S.G.). Á.C.-O. is the recipient of an FPU fellowship from Ministerio de Ciencia, Innovación y Universidades FPU13/00913 (Spanish government). L.S.-B. was funded by Wellcome grant number 098051.

Burkholderia is a multi-talented genus of Gram-negative bacteria, which in recent years has become increasingly recognised as a promising source of bioactive natural products. Metabolite profiling of Burkholderia gladioli BCC0238 showed that it produces the asymmetric lipopeptidiolide antibiotic icosalide A1, originally isolated from a fungus. Comparative bioinformatics analysis of several genome-sequenced B. gladioli isolates identified a gene encoding a nonribosomal peptide synthase (NRPS) with an unusual architecture that was predicted to be responsible for icosalide biosynthesis. Inactivation of this gene in B. gladioli BCC0238 abolished icosalide production. PCR analysis and sequencing of total DNA from the original fungal icosalide A1 producer revealed it has a B. gladioli strain associated with it that harbours an NRPS with an identical architecture to that responsible for icosalide A1 assembly in B. gladioli BCC0238. Sequence analysis of the icosalide NRPS indicated that it contains two chain-initiating condensation (CI) domains. One of these is appended to the N-terminus of module 1 - a common architecture for NRPSs involved in lipopeptide assembly. The other is embedded in module 3, immediately downstream of a putative chain-elongating condensation domain. Analysis of the reactions catalysed by a tridomain construct from module 3 of the NRPS using intact protein mass spectrometry showed that the embedded CI domain initiates assembly of a second lipopeptide chain, providing key insights into the mechanism for asymmetric diolide assembly. Open Access article ispartof: Chemical Science vol:10 issue:21 ispartof: location:England status: published

Group A Streptococcus (GAS; Streptococcus pyogenes) is a bacterial pathogen for which a commercial vaccine for humans is not available. Employing the advantages of high-throughput DNA sequencing technology to vaccine design, we have analyzed 2,083 globally sampled GAS genomes. The global GAS population structure reveals extensive genomic heterogeneity driven by homologous recombination and overlaid with high levels of accessory gene plasticity. We identified the existence of more than 290 clinically associated genomic phylogroups across 22 countries, highlighting challenges in designing vaccines of global utility. To determine vaccine candidate coverage, we investigated all of the previously described GAS candidate antigens for gene carriage and gene sequence heterogeneity. Only 15 of 28 vaccine antigen candidates were found to have both low naturally occurring sequence variation and high (>99%) coverage across this diverse GAS population. This technological platform for vaccine coverage determination is equally applicable to prospective GAS vaccine antigens identified in future studies., info:eu-repo/semantics/published

The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease., Claire Chewapreecha et al. combine 753 newly sequenced Thai Burkholderia pseudomallei isolates with 258 Australian isolates to identify genes associated with either clinical or environmental strains. They find 47 genes, 12 of which replicate in both locations, that may provide clues to the strategy used by this microbe to adapt to survive in wide range of ecological niches, including human hosts.

Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans. The disease is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Currently, the major hurdles facing disease control are incomplete understandings of both the mode of transmission and environmental reservoirs of M. ulcerans. As decades of spasmodic environmental sampling surveys have not brought us much closer to overcoming these hurdles, the Buruli ulcer research community has recently switched to using comparative genomics. The significance of our research is in how we used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our approach shows that it might be possible to use bacterial population genomics to assess the impact of health interventions, providing valuable feedback for managers of disease control programs in areas where health surveillance infrastructure is poor., Buruli ulcer is a neglected tropical disease of skin and subcutaneous tissue caused by infection with the pathogen Mycobacterium ulcerans. Many critical issues for disease control, such as understanding the mode of transmission and identifying source reservoirs of M. ulcerans, are still largely unknown. Here, we used genomics to reconstruct in detail the evolutionary trajectory and dynamics of M. ulcerans populations at a central African scale and at smaller geographical village scales. Whole-genome sequencing (WGS) data were analyzed from 179 M. ulcerans strains isolated from all Buruli ulcer foci in the Democratic Republic of the Congo, The Republic of Congo, and Angola that have ever yielded positive M. ulcerans cultures. We used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our phylogeographic analysis revealed one almost exclusively predominant sublineage of M. ulcerans that arose in Central Africa and proliferated in its different regions of endemicity during the Age of Discovery. We observed how the best sampled endemic hot spot, the Songololo territory, became an area of endemicity while the region was being colonized by Belgium (1880s). We furthermore identified temporal parallels between the observed past population fluxes of M. ulcerans from the Songololo territory and the timing of health policy changes toward control of the Buruli ulcer epidemic in that region. These findings suggest that an intervention based on detecting and treating human cases in an area of endemicity might be sufficient to break disease transmission chains, irrespective of other reservoirs of the bacterium. IMPORTANCE Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans. The disease is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Currently, the major hurdles facing disease control are incomplete understandings of both the mode of transmission and environmental reservoirs of M. ulcerans. As decades of spasmodic environmental sampling surveys have not brought us much closer to overcoming these hurdles, the Buruli ulcer research community has recently switched to using comparative genomics. The significance of our research is in how we used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our approach shows that it might be possible to use bacterial population genomics to assess the impact of health interventions, providing valuable feedback for managers of disease control programs in areas where health surveillance infrastructure is poor.

Gonorrhoea infections are on the increase and strains that are resistant to all antimicrobials used to treat the disease have been found worldwide. These observations encouraged the World Health Organization to include Neisseria gonorrhoeae on their list of high-priority organisms in need of new treatments. Fortunately, concurrent resistance to both antimicrobials used in dual therapy is still rare. The fight against antimicrobial resistance (AMR) must begin from an understanding of how it evolves and spreads in sexual networks. Genome-based analyses have allowed the study of the gonococcal population dynamics and transmission, giving a novel perspective on AMR gonorrhoea. Here, we will review past, present and future treatment options for gonorrhoea and explain how genomics is helping to increase our understanding of the changing AMR and transmission landscape. This article contains data hosted by Microreact.

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology.However, current whole-genome clustering and multi-locus typing approaches do not fully exploit core and accessory genomic variation, and cannot both automatically identify, and subsequently expand, clusters of significantly-similar isolates in large datasets and across species.Here we describe PopPUNK (Population Partitioning Using Nucleotide K-mers; https://poppunk.readthedocs.io/en/latest/). software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering.Variable-length k-mer comparisons are used to distinguish isolates’ divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using both simulated data and real datasets from ten taxonomically-widespread species. Connections between closely-related isolates of the same strain are robustly identified, despite variation in the discontinuous pairwise distance distributions that reflects species’ diverse evolutionary patterns. PopPUNK can process 103-104 genomes as single batch, with minimal memory use and runtimes up to 200-fold faster than existing methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to re-analyse all genomes de novo.This facilitates real-time surveillance with stable cluster naming and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualisation and online publication is streamlined through automatic output of results to multiple platforms.PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

[Background]: Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen., [Methods]: Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes., [Results]: Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p, [Conclusion]: Our data indicate potential zoonotic transmission of bTB in Ghana and hence calls for intensified public education on bTB, especially among risk groups., This work was supported by the Wellcome Trust Intermediate Fellowship awarded to DYM (Grant Number 097134/Z/11/Z).

The sequencing and comparative analysis of a collection of bacterial genomes from a single species or lineage of interest can lead to key insights into its evolution, ecology or epidemiology. The tool of choice for such a study is often to build a phylogenetic tree, and more specifically when possible a dated phylogeny, in which the dates of all common ancestors are estimated. Here we propose a new Bayesian methodology to construct dated phylogenies which is specifically designed for bacterial genomics. Unlike previous Bayesian methods aimed at building dated phylogenies, we consider that the phylogenetic relationships between the genomes have been previously evaluated using a standard phylogenetic method, which makes our methodology much faster and scalable. This two-steps approach also allows us to directly exploit existing phylogenetic methods that detect bacterial recombination, and therefore to account for the effect of recombination in the construction of a dated phylogeny. We analysed many simulated datasets in order to benchmark the performance of our approach in a wide range of situations. Furthermore, we present applications to three different real datasets from recent bacterial genomic studies. Our methodology is implemented in a R package called BactDating which is freely available for download athttps://github.com/xavierdidelot/BactDating.

e21546 Background: Four independent groups have demonstrated that the pre-treatment gut microbiome of cancer patients impacts the subsequent response to Immune Checkpoint Inhibitor (ICIs) therapy [1-4]. However, the patient’s outcome was linked to different bacteria in each study, which has limited the development of drug response biomarkers and clinic-first design of novel microbiome-based therapeutics. Methods: The Cambridge (UK) MELRESIST study includes a cohort of advanced melanoma patients receiving approved ICIs. Pretreatment stool samples from MELRESIST were analysed by Microbiotica using shotgun metagenomic sequencing. Microbiotica’s platform comprises the world’s leading Reference Genome Database to give the most comprehensive and precise mapping of the gut microbiome. Results: MELRESIST samples showed an overall difference in the microbiome composition between advanced melanoma patients and healthy donors, but not between patients who did or did not respond to ICIs. However, we did identify a discrete microbiome signature that differentiated responders and non-responders with an accuracy of 93%. We extended this signature by reanalysing three published melanoma cohorts [1-3] using the Microbiotica platform, and a propriety bioinformatic model. The resultant bacterial signature was very accurate at predicting response in all 4 published studies combined (91%), and each cohort individually (82-100%). We validated the model using independent validation cohorts and the signature using lung and renal cancer studies [4]. At the core of our microbiome signature was 9 bacteria most significantly associated with ICI efficacy. All 9 were overrepresented in patients who responded to immunotherapy suggesting high abundance of these bacteria is a central driver of ICI response. A consortium comprised of all 9 strains had very potent anti-tumor efficacy in a cancer syngeneic mouse model. The bacteria also demonstrate multiple interactions with primary human immune cells in vitro leading to dendritic cells activation, Cytotoxic T lymphocyte activation and tumor cell killing. These validate the potential of this consortium as a novel therapy for use in combination with ICIs. Conclusions: We have identified a unique microbiome signature predictive of ICI response in 4 independent melanoma cancer cohorts. This removes a major challenge to the field, and could represent a new highly accurate biomarker with clinical application. Nine core bacteria appear to be driving response, and demonstrate anti-tumor activity in vivo and in vitro. This consortium holds great potential as a co-therapy with ICIs. References:1 Matson V et al, Science (2018) 359:104; 2 Gopalakrishnan V et al, Science (2018) 359:97; 3 Frankel AE et al, Neoplasia (2017) 19:848; 4 Routy B et al, Science (2018) 359:91.

This paper contributes theoretical understanding to the learning processes adopted by firms for successful internationalization. Drawing from the internationalization process and organizational learning research, our longitudinal case study examines a firm’s learning processes as it grows internationally. In different phases and in different areas the firm learned at different paces, gradual and steep, and by responding retrospectively to past failings, and prospectively to anticipated challenges. From case data we identify and develop four constructs of internationalization learning: Gradual Retrospective, Steep Retrospective, Gradual Prospective, and Steep Prospective learning. The appropriateness of each depends on the urgency towards addressing challenges, and the firm’s learning capabilities. Firms can develop and steepen retrospective learning capability, enabling rapid identification and response to internationalization problems. Firms can also employ prospective learning to plan the development of knowledge capabilities ahead of future internationalization activity, to accelerate their learning, and increase the speed and effectiveness of internationalization.

Bacterial genomes vary extensively in terms of both gene content and gene sequence. This plasticity hampers the use of traditional SNP-based methods for identifying all genetic associations with phenotypic variation. Here we introduce a computationally scalable and widely applicable statistical method (SEER) for the identification of sequence elements that are significantly enriched in a phenotype of interest. SEER is applicable to tens of thousands of genomes by counting variable-length k-mers using a distributed string-mining algorithm. Robust options are provided for association analysis that also correct for the clonal population structure of bacteria. Using large collections of genomes of the major human pathogens Streptococcus pneumoniae and Streptococcus pyogenes, SEER identifies relevant previously characterized resistance determinants for several antibiotics and discovers potential novel factors related to the invasiveness of S. pyogenes. We thus demonstrate that our method can answer important biologically and medically relevant questions., Plasticity and clonal population structure in bacterial genomes can hinder traditional SNP-based genetic association studies. Here, Corander and colleagues present a method to identify variable-length sequence elements enriched in a phenotype of interest, and demonstrate its use in human pathogens.

Sequence-based typing (SBT), analogous to multilocus sequence typing (MLST), is the current “gold standard” typing method for investigation of legionellosis outbreaks caused by Legionella pneumophila . However, as common sequence types (STs) cause many infections, some investigations remain unresolved. In this study, various whole-genome sequencing (WGS)-based methods were evaluated according to published guidelines, including (i) a single nucleotide polymorphism (SNP)-based method, (ii) extended MLST using different numbers of genes, (iii) determination of gene presence or absence, and (iv) a kmer-based method. L. pneumophila serogroup 1 isolates ( n = 106) from the standard “typing panel,” previously used by the European Society for Clinical Microbiology Study Group on Legionella Infections (ESGLI), were tested together with another 229 isolates. Over 98% of isolates were considered typeable using the SNP- and kmer-based methods. Percentages of isolates with complete extended MLST profiles ranged from 99.1% (50 genes) to 86.8% (1,455 genes), while only 41.5% produced a full profile with the gene presence/absence scheme. Replicates demonstrated that all methods offer 100% reproducibility. Indices of discrimination range from 0.972 (ribosomal MLST) to 0.999 (SNP based), and all values were higher than that achieved with SBT (0.940). Epidemiological concordance is generally inversely related to discriminatory power. We propose that an extended MLST scheme with ∼50 genes provides optimal epidemiological concordance while substantially improving the discrimination offered by SBT and can be used as part of a hierarchical typing scheme that should maintain backwards compatibility and increase discrimination where necessary. This analysis will be useful for the ESGLI to design a scheme that has the potential to become the new gold standard typing method for L. pneumophila .

Objectives Gonorrhoea and MDR Neisseria gonorrhoeae remain public health concerns globally. Enhanced, quality-assured, gonococcal antimicrobial resistance (AMR) surveillance is essential worldwide. The WHO global Gonococcal Antimicrobial Surveillance Programme (GASP) was relaunched in 2009. We describe the phenotypic, genetic and reference genome characteristics of the 2016 WHO gonococcal reference strains intended for quality assurance in the WHO global GASP, other GASPs, diagnostics and research worldwide. Methods The 2016 WHO reference strains (n = 14) constitute the eight 2008 WHO reference strains and six novel strains. The novel strains represent low-level to high-level cephalosporin resistance, high-level azithromycin resistance and a porA mutant. All strains were comprehensively characterized for antibiogram (n = 23), serovar, prolyliminopeptidase, plasmid types, molecular AMR determinants, N. gonorrhoeae multiantigen sequence typing STs and MLST STs. Complete reference genomes were produced using single-molecule PacBio sequencing. Results The reference strains represented all available phenotypes, susceptible and resistant, to antimicrobials previously and currently used or considered for future use in gonorrhoea treatment. All corresponding resistance genotypes and molecular epidemiological types were described. Fully characterized, annotated and finished references genomes (n = 14) were presented. Conclusions The 2016 WHO gonococcal reference strains are intended for internal and external quality assurance and quality control in laboratory investigations, particularly in the WHO global GASP and other GASPs, but also in phenotypic (e.g. culture, species determination) and molecular diagnostics, molecular AMR detection, molecular epidemiology and as fully characterized, annotated and finished reference genomes in WGS analysis, transcriptomics, proteomics and other molecular technologies and data analysis.

OBJECTIVES: To elucidate the genome-based epidemiology and phylogenomics of azithromycin-resistant (MIC >2 mg/L) Neisseria gonorrhoeae strains collected in 2009-14 in Europe and clarify the azithromycin resistance mechanisms. METHODS: Seventy-five azithromycin-resistant (MIC 4 to >256 mg/L) N. gonorrhoeae isolates collected in 17 European countries during 2009-14 were examined using antimicrobial susceptibility testing and WGS. RESULTS: Thirty-six N. gonorrhoeae multi-antigen sequence typing STs and five phylogenomic clades, including 4-22 isolates from several countries per clade, were identified. The azithromycin target mutation A2059G (Escherichia coli numbering) was found in all four alleles of the 23S rRNA gene in all isolates with high-level azithromycin resistance (n = 4; MIC ≥256 mg/L). The C2611T mutation was identified in two to four alleles of the 23S rRNA gene in the remaining 71 isolates. Mutations in mtrR and its promoter were identified in 43 isolates, comprising isolates within the whole azithromycin MIC range. No mutations associated with azithromycin resistance were found in the rplD gene or the rplV gene and none of the macrolide resistance-associated genes [mef(A/E), ere(A), ere(B), erm(A), erm(B), erm(C) and erm(F)] were identified in any isolate. CONCLUSIONS: Clonal spread of relatively few N. gonorrhoeae strains accounts for the majority of the azithromycin resistance (MIC >2 mg/L) in Europe. The four isolates with high-level resistance to azithromycin (MIC ≥256 mg/L) were widely separated in the phylogenomic tree and did not belong to any of the main clades. The main azithromycin resistance mechanisms were the A2059G mutation (high-level resistance) and the C2611T mutation (low- and moderate-level resistance) in the 23S rRNA gene. The present work was supported by grants from theO¨ rebro County Council Research Committee and the Foundation for Medical Research at O¨ rebro University Hospital, Sweden info:eu-repo/semantics/publishedVersion

Microevolution resulted in considerable genotypic variation., Microevolution associated with emergence and expansion of new epidemic clones of bacterial pathogens holds the key to epidemiologic success. To determine microevolution associated with monophasic Salmonella Typhimurium during an epidemic, we performed comparative whole-genome sequencing and phylogenomic analysis of isolates from the United Kingdom and Italy during 2005–2012. These isolates formed a single clade distinct from recent monophasic epidemic clones previously described from North America and Spain. The UK monophasic epidemic clones showed a novel genomic island encoding resistance to heavy metals and a composite transposon encoding antimicrobial drug resistance genes not present in other Salmonella Typhimurium isolates, which may have contributed to epidemiologic success. A remarkable amount of genotypic variation accumulated during clonal expansion that occurred during the epidemic, including multiple independent acquisitions of a novel prophage carrying the sopE gene and multiple deletion events affecting the phase II flagellin locus. This high level of microevolution may affect antigenicity, pathogenicity, and transmission.

Molecularly imprinted polymers (MIPs) templated with either the phytoalexin, (E)-resveratrol, or its structural analog, 3,5-dihydroxy-N-(4-hydroxyphenyl)benzamide, have been used in tandem for the sequential extraction of (E)-resveratrol from aqueous peanut meal extracts in high purity and in near quantitative yields. Re-processing of the (E)-resveratrol-depleted peanut meal extract with the 3,5-dihydroxy-N-(4-hydroxyphenyl)benzamide imprinted MIP yielded additional polyphenolic components, identified as A-type procyanidins. Tandem liquid chromatography-electrospray ionization mass spectrometry confirmed the identity and purity of the isolated products. This study documents the advantages of tandem approaches with MIPs for the solid phase extraction and analysis of multiple bioactive compounds present in complex biomass waste streams.

Chlamydia trachomatis causes sexually transmitted infections and the blinding disease trachoma. Current data on C. trachomatis phylogeny show that there is only a single trachoma-causing clade, which is distinct from the lineages causing urogenital tract (UGT) and lymphogranuloma venerum diseases. Here we report the whole-genome sequences of ocular C. trachomatis isolates obtained from young children with clinical signs of trachoma in a trachoma endemic region of northern Australia. The isolates form two lineages that fall outside the classical trachoma lineage, instead being placed within UGT clades of the C. trachomatis phylogenetic tree. The Australian trachoma isolates appear to be recombinants with UGT C. trachomatis genome backbones, in which loci that encode immunodominant surface proteins (ompA and pmpEFGH) have been replaced by those characteristic of classical ocular isolates. This suggests that ocular tropism and association with trachoma are functionally associated with some sequence variants of ompA and pmpEFGH., Chlamydia trachomatis isolates causing a blinding disease (trachoma) form a single lineage that is different from the lineages causing urogenital infections. Here, Andersson et al. show however that trachoma isolates from Australia are more closely related to urogenital strains than to other trachoma isolates.

Global spread of aggressive mycobacteria Many mycobacteria, in addition to those causing leprosy and tuberculosis, are capable of infecting humans. Some can be particularly dangerous in patients suffering from immunosuppression or chronic disease, such as cystic fibrosis (CF). Bryant et al. observed clusters of near-identical isolates of drug-resistant Mycobacterium abscessus in patients reporting to CF clinics. The similarity of the isolates suggests transmission between patients, rather than environmental acquisition. Although this bacterium is renowned for its environmental resilience, the mechanism for its long-distance transmission among the global CF patient community remains a puzzle. Science , this issue p. 751

The ability to DNA fingerprint Mycobacterium bovis (M. bovis) isolates helped to define the role of wildlife in the persistnace of bovine tuberculosis in New Zealand. DNA fingerprinting results currently help to guide wildlife control measures and also aid in tracing the source of infections that result from movement of livestock. During the last 5 years we have developed the ability to distinguish New Zealand (NZ) M. bovis isolates by comparing the sequences of whole genome sequenced (WGS) M. bovis samples. WGS provides much higher resolution than our other established typing methods and greatly improved the definition of the regional localization of NZ M. bovis types. Three outbreak investigations are described and results demonstrate how WGS analysis has led to the confirmation of epidemiological sourcing of infection, to better definition of the source of new infections by ruling out other possible sources, and has revealed probable wildlife infection in an area considered to be free of infected wildlife. The routine use of WGS analyses for sourcing new M. bovis infections will be an important component of the strategy employed to eradicate bovine TB from NZ livestock and wildlife.

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K-mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates’ divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species’ diverse evolutionary patterns. PopPUNK can process 10(3)–10(4) genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

C. neoformansvar.grubii(C. neoformans) is an environmentally acquired pathogen causing 181 000 HIV-associated deaths each year. We used whole genome sequencing (WGS) to characterise 699 isolates, primarilyC. neoformansfrom HIV-infected patients, from 5 countries in Asia and Africa. We found that 91% of our clinical isolates belonged to one of three highly clonal sub-clades of VNIa, which we have termed VNIa-4, VNIa-5 and VNIa-93. Parsimony analysis revealed frequent, long distance transmissions ofC. neoformans; international transmissions took place on 13% of VNIa-4 branches, and intercontinental transmissions on 7% of VNIa-93 branches. The median length of within sub-clade internal branches was 3-6 SNPs, while terminal branches were 44.5-77.5 SNPs. The short median internal branches were partly driven by the large number (12-15% of internal branches) of polytomies in the within-sub-clade trees. To simultaneously explain our observation of no apparent molecular clock, short internal branches and frequent polytomies we hypothesise thatC. neoformansVNIa spends much of its time in the environment in a quiescent state, while, when it is sampled, it has almost always undergone an extended period of growth. Infections with VNIa-93 were associated with a significantly reduced risk of death by 10 weeks compared with infections with VNIa-4 (Hazard Ratio = 0.45, p = 0.003). We detected a recombination in the mitochondrial sequence of VNIa-5, suggesting that mitochondria could be involved in the propensity of this sub-clade to infect HIV-uninfected patients. These data highlight the insight into the biology and epidemiology of pathogenic fungi which can be gained from WGS data.

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is common, compromising gonorrhoea treatment internationally. Rapid characterisation of AMR strains could ensure appropriate and personalised treatment, and support identification and investigation of gonorrhoea outbreaks in nearly real-time. Whole-genome sequencing is ideal for investigation of emergence and dissemination of AMR determinants, predicting AMR, in the gonococcal population and spread of AMR strains in the human population. The novel, rapid and revolutionary long-read sequencer MinION is a small hand-held device that generates bacterial genomes within one day. However, accuracy of MinION reads has been suboptimal for many objectives and the MinION has not been evaluated for gonococci. In this first MinION study for gonococci, we show that MinION-derived sequences analysed with existing open-access, web-based sequence analysis tools are not sufficiently accurate to identify key gonococcal AMR determinants. Nevertheless, using an in house-developed CLC Genomics Workbench including de novo assembly and optimised BLAST algorithms, we show that 2D ONT-derived sequences can be used for accurate prediction of decreased susceptibility or resistance to recommended antimicrobials in gonococcal isolates. We also show that the 2D ONT-derived sequences are useful for rapid phylogenomic-based molecular epidemiological investigations, and, in hybrid assemblies with Illumina sequences, for producing contiguous assemblies and finished reference genomes.

The sexually transmitted pathogen Neisseria gonorrhoeae is regarded as being on the way to becoming an untreatable superbug. Despite its clinical importance, little is known about its emergence and evolution, and how this corresponds with the introduction of antimicrobials. We present a genome-based phylogeographic analysis of 419 gonococcal isolates from across the globe. Results indicate that modern gonococci originated in Europe or Africa as late as the 16thcentury and subsequently disseminated globally. We provide evidence that the modern gonococcal population has been shaped by antimicrobial treatment of sexually transmitted and other infections, leading to the emergence of two major lineages with different evolutionary strategies. The well-described multi-resistant lineage is associated with high rates of homologous recombination and infection in high-risk sexual networks where antimicrobial treatment is frequent. A second, multi-susceptible lineage associated with heterosexual networks, where asymptomatic infection is more common, was also identified, with potential implications for infection control.

Background The live, temperature-attenuated vaccine strain 1B of Chlamydia abortus, the aetiological agent of ovine enzootic abortion (OEA), has been implicated in cases of vaccine breakdown. The aim of this study was to understand the nature of this attenuation through sequencing of the vaccine parent strain (AB7) and the derived mutant strains 1B and 1H, as well as to clarify the role of the vaccine strain in causing disease through comparative whole genome analysis. Methods Whole genome sequencing was performed on: vaccine parent strain AB7; N-methyl-N′-nitro-N-nitrosoguanidine (NTG)-induced temperature attenuated mutant strain 1B grown from the commercial live vaccines Cevac Chlamydia and Enzovax; strain 1H a reverted NTG mutant; and 5 strains isolated from cases of OEA originating from animals from the original vaccine safety trial (2 strains) or from vaccinated ewes or ewes exposed to vaccinated animals (3 strains). Results We confirmed that AB7 is in a different lineage from the reference strain S26/3. The genome of vaccine strain 1B contains ten single nucleotide polymorphisms (SNPs) created by the NTG treatment, which are identical to those found in strain 1H. The strains from OEA cases also cluster phylogenetically very tightly with these vaccine strains. Conclusions The results show that C. abortus vaccine strain 1B has an identical genome sequence to the non-attenuated “reverted mutant” strain 1H. Thus, the protection of the 1B vaccine is unlikely to be due to the NTG induced SNPs and is more likely caused by the administration of high doses of C. abortus elementary bodies that stimulate protective immunity. Vaccine-identical strains were also isolated from cases of disease, as well as strains which had acquired 1–3 SNPs, including an animal that had not been vaccinated with either of the commercial live OEA vaccines, indicating that the 1B vaccine strain may be circulating and causing disease.

The emergence of Clostridium difficile as a significant human diarrheal pathogen is associated with the production of highly transmissible spores and the acquisition of antimicrobial resistance genes (ARGs) and virulence factors. Unlike the hospital-associated C. difficile RT027 lineage, the community-associated C. difficile RT078 lineage is isolated from both humans and farm animals; however, the geographical population structure and transmission networks remain unknown. Here, we applied whole-genome phylogenetic analysis of 248 C. difficile RT078 strains from 22 countries. Our results demonstrate limited geographical clustering for C. difficile RT078 and extensive coclustering of human and animal strains, thereby revealing a highly linked intercontinental transmission network between humans and animals. Comparative whole-genome analysis reveals indistinguishable accessory genomes between human and animal strains and a variety of antimicrobial resistance genes in the pangenome of C. difficile RT078. Thus, bidirectional spread of C. difficile RT078 between farm animals and humans may represent an unappreciated route disseminating antimicrobial resistance genes between humans and animals. These results highlight the importance of the “One Health” concept to monitor infectious disease emergence and the dissemination of antimicrobial resistance genes.

The Neutron Transport Equation (NTE) describes the flux of neutrons through inhomogeneous fissile medium. Whilst well treated in the nuclear physics literature (cf. [9, 27]), the NTE has had a somewhat scattered treatment in mathematical literature with a variety of different approaches (cf. [8, 25]). Within a probabilistic framework it has somewhat undeservingly received little attention in recent years; nonetheless, probabilistic treatments can be found see for example [19, 26, 24, 29, 4, 3]. In this article our aim is threefold. First we want to introduce a slightly more general setting for the NTE, which gives a more complete picture of the different species of particle and radioactive fluxes that are involved in fission. Second we consolidate the classical c0-semigroup approach to solving the NTE with the method of stochastic representation which involves expectation semigroups. Third we provide the leading asymptotic of our multi-species NTE, which will turn out to be crucial for further stochastic analysis of the NTE in forthcoming work [6, 5]. The methodology used in this paper harmonises the culture of expectation semigroup analysis from the theory of stochastic processes against c0-semigroup theory from functional analysis. In this respect, our presentation is thus part review of existing theory and part presentation of new research results based on generalisation of existing results., Comment: 2 figures

The evolution and global transmission of antimicrobial resistance has been well documented in Gram-negative bacteria and healthcare-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. Here, we trace the recent origins and global spread of a multidrug resistant, community-associatedStaphylococcus aureuslineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data shows that the clone emerged on the Indian subcontinent in the early 1970s and disseminated rapidly in the 1990s. Short-term outbreaks in community and healthcare settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the divergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional healthcare-associated clones with the epidemiological transmission of community-associated MRSA. Our study demonstrates the importance of whole genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.ImportanceThe Bengal Bay clone (ST772) is a community-acquired and multidrug-resistantStaphylococcus aureuslineage first isolated from Bangladesh and India in 2004. In this study, we show that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally resulting in small-scale community and healthcare outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug-resistance of healthcare-associatedS. aureuslineages. This study demonstrates the importance of whole genome sequencing for the surveillance of highly antibiotic resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world.

Background Although Mycoplasma genitalium is a common sexually transmitted pathogen causing clinically distinct diseases both in male and females, few genomes have been sequenced up to now, due mainly to its fastidious nature and slow growth. Hence, we lack a robust phylogenetic framework to provide insights into the population structure of the species. Currently our understanding of the nature and diversity of M. genitalium relies on molecular tests targeting specific genes or regions of the genome and knowledge is limited by a general under-testing internationally. This is set against a background of drug resistance whereby M. genitalium has developed resistance to mainly all therapeutic antimicrobials. Results We sequenced 28 genomes of Mycoplasma genitalium from temporally (1980–2010) and geographically (Europe, Japan, Australia) diverse sources. All the strain showed essentially the same genomic content without any accessory regions found. However, we identified extensive recombination across their genomes with a total of 25 regions showing heightened levels of SNP density. These regions include the MgPar loci, associated with host interactions, as well as other genes that could also be involved in this role. Using these data, we generated a robust phylogeny which shows that there are two main clades with differing degrees of genomic variability. SNPs found in region V of 23S rRNA and parC were consistent with azithromycin/erythromycin and fluoroquinolone resistances, respectively, and with their phenotypic MIC data. Conclusions The sequence data here generated is essential for designing rational approaches to type and track Mycoplasma genitalium as antibiotic resistance increases. It represents a first approach to its population genetics to better appreciate the role of this organism as a sexually transmitted pathogen. Electronic supplementary material The online version of this article (10.1186/s12864-017-4399-6) contains supplementary material, which is available to authorized users.

A major outbreak of whooping cough, or pertussis, occurred in 2012 in the United Kingdom (UK), with nearly 10 000 laboratory-confirmed cases and 14 infant deaths attributed to pertussis. A worldwide resurgence of pertussis has been linked to switch to the use of acellular pertussis vaccines and the evolution of Bordetella pertussis away from vaccine-mediated immunity. We have conducted genomic analyses of multiple strains from the UK outbreak. We show that the UK outbreak was polyclonal in nature, caused by multiple distinct but closely related strains. Importantly, we demonstrate that acellular vaccine antigen-encoding genes are evolving at higher rates than other surface protein-encoding genes. This was true even prior to the introduction of pertussis vaccines but has become more pronounced since the introduction of the current acellular vaccines. The fast evolution of vaccine antigen-encoding genes has serious consequences for the ability of current vaccines to continue to control pertussis.

Mycobacterium africanum(Maf) causes up to half of human tuberculosis in West Africa, but little is known on this pathogen. We compared the genomes of 253Mafclinical isolates from Ghana, including both L5 and L6. We found that the genomic diversity of L6 was higher than in L5, and the selection pressures differed between both groups. Regulatory proteins appeared to evolve neutrally in L5 but under purifying selection in L6. Conversely, human T cell epitopes were under purifying selection in L5, but under positive selection in L6. Although only 10% of the T cell epitopes were variable, mutations were mostly lineage-specific. Our findings indicate thatMafL5 and L6 are genomically distinct, possibly reflecting different ecological niches.

Antimicrobial resistance (AMR) is one of the major threats to human and animal health worldwide, yet few high-throughput tools exist to analyse and predict the resistance of a bacterial isolate from sequencing data. Here we present a new tool, ARIBA, that identifies AMR-associated genes and single nucleotide polymorphisms directly from short reads, and generates detailed and customisable output. The accuracy and advantages of ARIBA over other tools are demonstrated on three datasets from Gram-positive and Gram-negative bacteria, with ARIBA outperforming existing methods. ARIBA is available at https://github.com/sanger-pathogens/ariba.

Background Horizontal transfer of mobile genetic elements (MGEs) that carry virulence and antimicrobial resistance genes mediates the evolution of methicillin-resistant Staphylococcus aureus, and the emergence of new MRSA clones. Most MRSA lineages show an association with specific MGEs and the evolution of MGE composition following clonal expansion has not been widely studied. Results We investigated the genomes of 1193 S. aureus bloodstream isolates, 1169 of which were MRSA, collected in the UK and the Republic of Ireland between 2001 and 2010. The majority of isolates belonged to clonal complex (CC)22 (n = 923), which contained diverse MGEs including elements that were found in other MRSA lineages. Several MGEs showed variable distribution across the CC22 phylogeny, including two antimicrobial resistance plasmids (pWBG751-like and SAP078A-like, carrying erythromycin and heavy metal resistance genes, respectively), a pathogenicity island carrying the enterotoxin C gene and two phage types Sa1int and Sa6int. Multiple gains and losses of these five MGEs were identified in the CC22 phylogeny using ancestral state reconstruction. Analysis of the temporal distribution of the five MGEs between 2001 and 2010 revealed an unexpected reduction in prevalence of the two plasmids and the pathogenicity island, and an increase in the two phage types. This occurred across the lineage and was not correlated with changes in the relative prevalence of CC22, or of any sub-lineages within in. Conclusions Ancestral state reconstruction coupled with temporal trend analysis demonstrated that epidemic MRSA CC22 has an evolving MGE composition, and indicates that this important MRSA lineage has continued to adapt to changing selective pressure since its emergence. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4065-z) contains supplementary material, which is available to authorized users.

This work was supported by the Animal Health Trust. The zoonotic bacterium Streptococcus equi subsp. zooepidemicus (S. zooepidemicus) is a diverse, opportunistic pathogen that can cause mastitis in dairy sheep and goats. We used multilocus sequence typing (MLST) to define the genetic diversity of 60 isolates of S. zooepidemicus, which were recovered from sheep and goats in Spain between 2003 and 2010. We identify a novel clonal complex based on sequence type (ST), ST-236, which accounted for 39 of the 60 isolates. A representative ST-236 strain, S. zooepidemicus strain C7 (SzC7), was sequenced and interrogated for the presence of novel nutritional uptake or utilisation systems, the acquisition of which have previously been shown to be important for environmental adaptation in other streptococcal pathogens. A novel phosphoenolpyruvate sugar phosphotransferase system (PTS), which enabled the utilisation of arbutin, was identified. Functionality of the PTS was confirmed following deletion of the PTS from SzC7. Arbutin is found in multiple animal foodstuffs and we propose that the ability to utilise arbutin may have conferred a selective advantage to strains infecting animals, the diet of which contains this sugar. Postprint

Countries of the African 'meningitis belt' are susceptible to meningococcal meningitis outbreaks. While in the past major epidemics have been primarily caused by serogroup A meningococci, W strains are currently responsible for most of the cases. After an epidemic in Mecca in 2000, W:ST-11 strains have caused many outbreaks worldwide. An unrelated W:ST-2881 clone was described for the first time in 2002, with the first meningitis cases caused by these bacteria reported in 2003. Here we describe results of a comparative whole-genome analysis of 74 W:ST-2881 strains isolated within the framework of two longitudinal colonization and disease studies conducted in Ghana and Burkina Faso. Genomic data indicate that the W:ST-2881 clone has emerged from Y:ST-175(CC175) bacteria by capsule switching. The circulating W:ST-2881 populations were composed of a variety of closely related but distinct genomic variants with no systematic differences between colonization and disease isolates. Two distinct and geographically clustered phylogenetic clonal variants were identified in Burkina Faso and a third in Ghana. On the basis of the presence or absence of 17 recombination fragments, the Ghanaian variant could be differentiated into five clusters. All 25 Ghanaian disease isolates clustered together with 23 out of 40 Ghanaian isolates associated with carriage within one cluster, indicating that W:ST-2881 clusters differ in virulence. More than half of the genes affected by horizontal gene transfer encoded proteins of the 'cell envelope' and the 'transport/binding protein' categories, which indicates that exchange of non-capsular antigens plays an important role in immune evasion.

Introduction People living in Britain who have sex abroad are more likely to report sexual behaviour that puts them at greater risk of acquiring STIs, including Neisseria gonorrhoeae (NG). Antimicrobial resistant (AMR) NG is a global public health concern, which may emerge de novo or be imported to the UK when individuals infected abroad have subsequent sexual partners at home. We investigated whether patients who reported sex outside the UK (‘sex abroad’) were more or less likely to be diagnosed with AMR NG. Methods Logistic regression was used to model the association between reported recent sex abroad and decreased susceptibility (DS) to ceftriaxone (MIC (mg/L)>0.015) and cefixime (0.125) and azithromycin AMR (>1) stratifying by sexual orientation (men who have sex with men (MSM) and heterosexual men and women) from isolates in England and Wales collected within the Gonococcal Resistance to Antimicrobials Surveillance Programme, 2004–2015. Results Over 10% of MSM and heterosexuals reported sex abroad. Among heterosexuals, infection with a strain of NG with DS to ceftriaxone was associated with sex abroad after adjusting for potential confounders: ceftriaxone (DS prevalence, adjusted odds ratio (95% confidence interval)): 14%, 1.8 (1.3–2.3). Infection with NG DS/AMR to cefixime or azithromycin was not associated with reported sex abroad after adjusting for potential confounders: cefixime 4%, 1.6 (0.9–2.7); azithromycin 2%, 1.5 (0.7–3.3). For MSM, no association was found between infections with DS/AMR NG and sex abroad. Conclusion In the UK, heterosexuals with NG infection who report sex abroad are at a higher risk of DS to ceftriaxone, suggesting that sex abroad might be the source of some AMR NG within heterosexual networks and highlighting the importance of condom use for travellers. In contrast, DS/AMR NG was not associated with sex abroad among MSM, which might reflect established AMR within MSM networks in the UK. Genetic comparison of these isolates using whole genome sequencing might further elucidate how AMR NG is imported and disseminated in the UK.

Introduction In England, men who have sex with men (MSM) who are HIV-positive are disproportionately affected by STIs, in part probably due to HIV sero-adaptive behaviours. Neisseria gonorrhoeae (NG) is of particular concern because treatment is threatened by antimicrobial resistance (AMR). In England, AMR NG has typically spread rapidly within sexual networks of MSM. We investigated whether the emergence and/or spread of AMR NG was associated with HIV-positive status. Methods The prevalence of NG decreased susceptibility (DS) to ceftriaxone (MIC (mg/L)≥0.015), cefixime (≥0.125), and azithromycin AMR (≥1) from 2004–2015 was plotted by HIV status to investigate the emergence of DS/AMR using data from England and Wales collected within the Gonococcal Resistance to Antimicrobials Surveillance Programme. Differences were assessed using the Kolmogorov-Smirnov (KS) test. Logistic regression was used to model the association between HIV status and susceptibility to these antimicrobials in separate models adjusting for year. Results Among all 5,630 MSM with NG, 25% of samples had DS/AMR to ceftriaxone, 8% to cefixime and 3% to azithromycin. A third (2024/5630) of MSM were HIV-positive. The distribution of prevalence of NG DS/AMR to ceftriaxone, cefixime and azithromycin was similar in HIV-positive and HIV-negative MSM across 2004–2015 (p>0.05 for each antimicrobial). In the logistic regression models, HIV-positive MSM were as likely as HIV-negative MSM to be infected with NG DS to ceftriaxone (DS/AMR prevalence in HIV-positive MSM vs HIV-negative MSM, adjusted odds ratio [95% confidence interval]) 25% vs 25%, 1.0 [0.9–1.1], cefixime 7% vs 8%, 1.1 [0.9–1.4] or azithromycin: 3% vs 3%, 0.9 [0.6–1.2]. Conclusion From these epidemiological data there is no evidence that MSM with HIV are at greater risk of DS/AMR NG compared to those without HIV. Whole genome sequencing will assist further investigations to explore relatedness of isolates and understand whether distinct populations of NG are spread more efficiently within sexual networks of HIV-positive MSM.

We study a dyadic branching Brownian motion on the real line with absorption at 0, drift $\mu \in \mathbb{R}$ and started from a single particle at position $x>0.$ When $\mu$ is large enough so that the process has a positive probability of survival, we consider $K(t),$ the number of individuals absorbed at 0 by time $t$ and for $s\ge 0$ the functions $\omega_s(x):= \mathbb{E}^x[s^{K(\infty)}].$ We show that $\omega_s1$ and we study the properties of these functions. Furthermore, for $s=0, \omega(x) := \omega_0(x) =\mathbb{P}^x(K(\infty)=0)$ is the cumulative distribution function of the all time minimum of the branching Brownian motion with drift started at 0 without absorption. We give three descriptions of the family $\omega_s, s\in [0,s_0]$ through a single pair of functions, as the two extremal solutions of the Kolmogorov-Petrovskii-Piskunov (KPP) traveling wave equation on the half-line, through a martingale representation and as an explicit series expansion. We also obtain a precise result concerning the tail behavior of $K(\infty)$. In addition, in the regime where $K(\infty)>0$ almost surely, we show that $u(x,t) := \mathbb{P}^x(K(t)=0)$ suitably centered converges to the KPP critical travelling wave on the whole real line., Comment: Grant information added

Legionella pneumophila is an environmental bacterium and the causative agent of Legionnaires’ disease. Previous genomic studies have shown that recombination accounts for a high proportion (>96%) of diversity within several major disease-associated sequence types (STs) of L. pneumophila. This suggests that recombination represents a potentially important force shaping adaptation and virulence. Despite this, little is known about the biological effects of recombination in L. pneumophila, particularly with regards to homologous recombination (whereby genes are replaced with alternative allelic variants). Using newly available population genomic data, we have disentangled events arising from homologous and non-homologous recombination in six major disease-associated STs of L. pneumophila (subsp. pneumophila), and subsequently performed a detailed characterisation of the dynamics and impact of homologous recombination. We identified genomic “hotspots” of homologous recombination that include regions containing outer membrane proteins, the lipopolysaccharide (LPS) region and Dot/Icm effectors, which provide interesting clues to the selection pressures faced by L. pneumophila. Inference of the origin of the recombined regions showed that isolates have most frequently imported DNA from isolates belonging to their own clade, but also occasionally from other major clades of the same subspecies. This supports the hypothesis that the possibility for horizontal exchange of new adaptations between major clades of the subspecies may have been a critical factor in the recent emergence of several clinically important STs from diverse genomic backgrounds. However, acquisition of recombined regions from another subspecies, L. pneumophila subsp. fraseri, was rarely observed, suggesting the existence of a recombination barrier and/or the possibility of ongoing speciation between the two subspecies. Finally, we suggest that multi-fragment recombination may occur in L. pneumophila, whereby multiple non-contiguous segments that originate from the same molecule of donor DNA are imported into a recipient genome during a single episode of recombination., Author summary Legionella pneumophila is an environmental bacterium that causes Legionnaires’ disease, a serious and potentially fatal pneumonia. Previous studies have shown that members of this species undergo a process called recombination, whereby DNA is imported from another bacterial cell into the recipient genome. The imported DNA can either replace an equivalent segment of the genome (homologous recombination) or can comprise novel genes that are new to the recipient genome (non-homologous recombination). Whilst recombination plays an undoubtedly important role in L. pneumophila evolution, accounting for more than 96% of the diversity observed within some lineages, little is known about its biological impact. In this study, we performed a detailed characterisation of the dynamics and effect of homologous recombination on L. pneumophila evolution in six clinically important lineages of L. pneumophila. We identified “hotspot” regions of the genome in which an excess of homologous recombination events was observed, which provided important clues to the selection pressures faced by L. pneumophila. By determining the donors of the recombined regions, we also revealed that recombination has occurred most frequently between isolates from the same clade, but also occurred between isolates from different major clades. This demonstrates the possibility of new adaptations arising in one lineage and being transferred to another distantly related lineage, which we predict has been an important factor in the emergence of several major disease-causing strains from diverse genomic backgrounds.

Background Chlamydia abortus (formerly Chlamydophila abortus) is an economically important livestock pathogen, causing ovine enzootic abortion (OEA), and can also cause zoonotic infections in humans affecting pregnancy outcome. Large-scale genomic studies on other chlamydial species are giving insights into the biology of these organisms but have not yet been performed on C. abortus. Our aim was to investigate a broad collection of European isolates of C. abortus, using next generation sequencing methods, looking at diversity, geographic distribution and genome dynamics. Results Whole genome sequencing was performed on our collection of 57 C. abortus isolates originating primarily from the UK, Germany, France and Greece, but also from Tunisia, Namibia and the USA. Phylogenetic analysis of a total of 64 genomes shows a deep structural division within the C. abortus species with a major clade displaying limited diversity, in addition to a branch carrying two more distantly related Greek isolates, LLG and POS. Within the major clade, seven further phylogenetic groups can be identified, demonstrating geographical associations. The number of variable nucleotide positions across the sampled isolates is significantly lower than those published for C. trachomatis and C. psittaci. No recombination was identified within C. abortus, and no plasmid was found. Analysis of pseudogenes showed lineage specific loss of some functions, notably with several Pmp and TMH/Inc proteins predicted to be inactivated in many of the isolates studied. Conclusions The diversity within C. abortus appears to be much lower compared to other species within the genus. There are strong geographical signatures within the phylogeny, indicating clonal expansion within areas of limited livestock transport. No recombination has been identified within this species, showing that different species of Chlamydia may demonstrate different evolutionary dynamics, and that the genome of C. abortus is highly stable. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3657-y) contains supplementary material, which is available to authorized users.

An epidemic of post-surgical wound infections, caused by a non-tuberculous mycobacterium, has been on-going in Brazil. It has been unclear whether one or multiple lineages are responsible and whether their wide geographical distribution across Brazil is due to spread from a single point source or is the result of human-mediated transmission. 188 isolates, collected from nine Brazilian states, were whole genome sequenced and analysed using phylogenetic and comparative genomic approaches. The isolates from Brazil formed a single clade, which was estimated to have emerged in 2003. We observed temporal and geographic structure within the lineage that enabled us to infer the movement of sub-lineages across Brazil. The genome size of the Brazilian lineage was reduced relative to most strains in the three subspecies of Mycobacterium abscessus and contained a novel plasmid, pMAB02, in addition to the previously described pMAB01 plasmid. One lineage, which emerged just prior to the initial outbreak, is responsible for the epidemic of post-surgical wound infections in Brazil. Phylogenetic analysis indicates that multiple transmission events led to its spread. The presence of a novel plasmid and the reduced genome size suggest that the lineage has undergone adaptation to the surgical niche.