- Volume 4, Issue 3, 2018

In December 2013, Public Health England (PHE) observed an increase in the number of cases of Shigella sonnei linked to the Orthodox Jewish Community (OJC). Ultimately, 52 cases of S. sonnei phage type (PT) P and PT7 were notified between November 2013 and July 2014. Whole-genome sequencing (WGS) was performed on a HiSeq 2500 platform (Illumina) on isolates of S. sonnei submitted to PHE during the investigation. Quality trimmed sequence reads were mapped to a reference genome using BWA-MEM, and single-nucleotide polymorphisms (SNPs) were identified using GATK2. Analysis of the core genome SNP positions (>90 % consensus, minimum depth 10×, MQ≥30) revealed that isolates linked to the outbreak could be categorized as members of distinct monophyletic clusters (MPCs) representing concurrent regional outbreaks occurring in the OJCs across the United Kingdom. A dated phylogeny predicted the date of the most recent common ancestor of the MPCs to be approximately 3.1 years previously [95 % highest posterior density (HPD), 2.4–3.4]. Isolates of S. sonnei from cases from the OJCs in Israel included in the phylogeny, branched from nodes basal to the UK OJC outbreak clusters, indicating they were ancestral to the UK OJC isolates, and that the UK isolates represented multiple importations of S. sonnei into the UK population from Israel. The level of discrimination exhibited by WGS facilitated the identification of clusters of isolates within the closely related bacterial populations circulating in the OJC that may be linked to a unique point sources or transmission routes, thus enabling a more appropriate public health response and targeted interventions.

There is growing evidence that patients with Clostridiumdifficile-associated diarrhoea often acquire their infecting strain before hospital admission. Wastewater is known to be a potential source of surface water that is contaminated with C. difficile spores. Here, we describe a study that used genome sequencing to compare C. difficile isolated from multiple wastewater treatment plants across the East of England and from patients with clinical disease at a major hospital in the same region. We confirmed that C. difficile from 65 patients were highly diverse and that most cases were not linked to other active cases in the hospital. In total, 186 C. difficile isolates were isolated from effluent water obtained from 18 municipal treatment plants at the point of release into the environment. Whole genome comparisons of clinical and environmental isolates demonstrated highly related populations, and confirmed extensive release of toxigenic C. difficile into surface waters. An analysis based on multilocus sequence types (STs) identified 19 distinct STs in the clinical collection and 38 STs in the wastewater collection, with 13 of 44 STs common to both clinical and wastewater collections. Furthermore, we identified five pairs of highly similar isolates (≤2 SNPs different in the core genome) in clinical and wastewater collections. Strategies to control community acquisition should consider the need for bacterial control of treated wastewater.

Salmonella enterica serovar Enteritidis is a major agent of foodborne diseases worldwide. In Uruguay, this serovar was almost negligible until the mid 1990s but since then it has become the most prevalent. Previously, we characterized a collection of strains isolated from 1988 to 2005 and found that the two oldest strains were the most genetically divergent. In order to further characterize these strains, we sequenced and annotated eight genomes including those of the two oldest isolates. We report on the identification and characterization of a novel 44 kbp Salmonella prophage found exclusively in these two genomes. Sequence analysis reveals that the prophage is a mosaic, with homologous regions in different Salmonella prophages. It contains 60 coding sequences, including two genes, gogB and sseK3, involved in virulence and modulation of host immune response. Analysis of serovar Enteritidis genomes available in public databases confirmed that this prophage is absent in most of them, with the exception of a group of 154 genomes. All 154 strains carrying this prophage belong to the same sequence type (ST-1974), suggesting that its acquisition occurred in a common ancestor. We tested this by phylogenetic analysis of 203 genomes representative of the intraserovar diversity. The ST-1974 forms a distinctive monophyletic lineage, and the newly described prophage is a phylogenetic signature of this lineage that could be used as a molecular marker. The phylogenetic analysis also shows that the major ST (ST-11) is polyphyletic and might have given rise to almost all other STs, including ST-1974.

Extensively drug-resistant Klebsiella pneumoniae (XDR-KP) infections cause high mortality and are disseminating globally. Identifying the genetic basis underpinning resistance allows for rapid diagnosis and treatment. XDR isolates sourced from Greece and Brazil, including 19 polymyxin-resistant and five polymyxin-susceptible strains, were subjected to whole genome sequencing. Seventeen of the 19 polymyxin-resistant isolates harboured variations upstream or within mgrB. The most common mutation identified was an insertion at nucleotide position 75 in mgrB via an ISKpn26-like element in the ST258 lineage and ISKpn13 in one ST11 isolate. Three strains acquired an IS1 element upstream of mgrB and another strain had an ISKpn25 insertion at 133 bp. Other isolates had truncations (C28STOP, Q30STOP) or a missense mutation (D29E) affecting mgrB. Complementation assays revealed all mgrB perturbations contributed to resistance. Missense mutations in phoQ (T281M, G385C) were also found to facilitate resistance. Several variants in phoPQ co-segregating with the ISKpn26-like insertion were identified as potential partial suppressor mutations. Three ST258 samples were found to contain subpopulations with different resistance-conferring mutations, including the ISKpn26-like insertion colonizing with a novel mutation in pmrB (P158R), both confirmed via complementation assays. These findings highlight the broad spectrum of chromosomal modifications which can facilitate and regulate resistance against polymyxins in K. pneumoniae.

In Bacteria, a working consensus of species circumscription may have been reached and one of the most prominent criteria is high average nucleotide identity (ANI). ANI in effect groups strains that may recombine more or less frequently, depending on their biology, as opposed to rare interspecies gene transfer. For bacteriophages, which show various lifestyles, the nature of the fundamental natural unit, if it exists in a biological sense, is not well understood and defined. The approaches based on dot-plots are useful to group similar bacteriophages, yet are not quantitative and use arbitrarily set cut-offs. Here, we focus on lytic Myoviridae and test the ANI metric for group delineation. We show that ANI-based groups are in agreement with the International Committee on Taxonomy of Viruses (ICTV) classification and already established dot-plot groups, which are occasionally further refined owing to higher resolution of ANI analysis. Furthermore, these groups are separated among themselves by clear ANI discontinuities. Their members readily exchange core genes with each other while they do not with bacteriophages of other ANI groups, not even with the most similar. Thus, ANI-delineated groups may represent the natural units in lytic Myoviridae evolution with features that resemble those encountered in bacterial species.

The Arctic is being disproportionally affected by climate change compared with other geographic locations, and is currently experiencing unprecedented melt rates. The Greenland Ice Sheet (GrIS) can be regarded as the largest supraglacial ecosystem on Earth, and ice algae are the dominant primary producers on bare ice surfaces throughout the course of a melt season. Ice-algal-derived pigments cause a darkening of the ice surface, which in turn decreases albedo and increases melt rates. The important role of ice algae in changing melt rates has only recently been recognized, and we currently know little about their community compositions and functions. Here, we present the first analysis of ice algal communities across a 100 km transect on the GrIS by high-throughput sequencing and subsequent oligotyping of the most abundant taxa. Our data reveal an extremely low algal diversity with Ancylonema nordenskiöldii and a Mesotaenium species being by far the dominant taxa at all sites. We employed an oligotyping approach and revealed a hidden diversity not detectable by conventional clustering of operational taxonomic units and taxonomic classification. Oligotypes of the dominant taxa exhibit a site-specific distribution, which may be linked to differences in temperatures and subsequently the extent of the melting. Our results help to better understand the distribution patterns of ice algal communities that play a crucial role in the GrIS ecosystem.

Members of the Mycobacterium tuberculosis complex (MTBC) are the causative agents of tuberculosis in a range of mammals, including humans. A key feature of MTBC pathogens is their high degree of genetic identity yet distinct host tropism. Notably, while Mycobacterium bovis is highly virulent and pathogenic for cattle, the human pathogen M. tuberculosis is attenuated in cattle. Previous research also suggests that host preference amongst MTBC members has a basis in host innate immune responses. To explore MTBC host tropism, we present in-depth profiling of the MTBC reference strains M. bovis AF2122/97 and M. tuberculosis H37Rv at both the global transcriptional and the translational level via RNA-sequencing and SWATH MS. Furthermore, a bovine alveolar macrophage infection time course model was used to investigate the shared and divergent host transcriptomic response to infection with M. tuberculosis H37Rv or M. bovis AF2122/97. Significant differential expression of virulence-associated pathways between the two bacilli was revealed, including the ESX-1 secretion system. A divergent transcriptional response was observed between M. tuberculosis H37Rv and M. bovis AF2122/97 infection of bovine alveolar macrophages, in particular cytosolic DNA-sensing pathways at 48 h post-infection, and highlights a distinct engagement of M. bovis with the bovine innate immune system. The work presented here therefore provides a basis for the identification of host innate immune mechanisms subverted by virulent host-adapted mycobacteria to promote their survival during the early stages of infection.

Acinetobacter baumannii is a common causative agent of hospital-acquired infections and a leading cause of infection in burns patients. Carbapenem-resistant A. baumannii is considered a major public-health threat and has been identified by the World Health Organization as the top priority organism requiring new antimicrobials. The most common mechanism for carbapenem resistance in A. baumannii is via horizontal acquisition of carbapenemase genes. In this study, we sampled 20 A. baumannii isolates from a patient with extensive burns, and characterized the evolution of carbapenem resistance over a 45 day period via Illumina and Oxford Nanopore sequencing. All isolates were multidrug resistant, carrying two genomic islands that harboured several antibiotic-resistance genes. Most isolates were genetically identical and represented a single founder genotype. We identified three novel non-synonymous substitutions associated with meropenem resistance: F136L and G288S in AdeB (part of the AdeABC efflux pump) associated with an increase in meropenem MIC to ≥8 µg ml−1; and A515V in FtsI (PBP3, a penicillin-binding protein) associated with a further increase in MIC to 32 µg ml−1. Structural modelling of AdeB and FtsI showed that these mutations affected their drug-binding sites and revealed mechanisms for meropenem resistance. Notably, one of the adeB mutations arose prior to meropenem therapy but following ciprofloxacin therapy, suggesting exposure to one drug whose resistance is mediated by the efflux pump can induce collateral resistance to other drugs to which the bacterium has not yet been exposed.

Increasingly rich metadata are now being linked to samples that have been whole-genome sequenced. However, much of this information is ignored. This is because linking this metadata to genes, or regions of the genome, usually relies on knowing the gene sequence(s) responsible for the particular trait being measured and looking for its presence or absence in that genome. Examples of this would be the spread of antimicrobial resistance genes carried on mobile genetic elements (MGEs). However, although it is possible to routinely identify the resistance gene, identifying the unknown MGE upon which it is carried can be much more difficult if the starting point is short-read whole-genome sequence data. The reason for this is that MGEs are often full of repeats and so assemble poorly, leading to fragmented consensus sequences. Since mobile DNA, which can carry many clinically and ecologically important genes, has a different evolutionary history from the host, its distribution across the host population will, by definition, be independent of the host phylogeny. It is possible to use this phenomenon in a genome-wide association study to identify both the genes associated with the specific trait and also the DNA linked to that gene, for example the flanking sequence of the plasmid vector on which it is encoded, which follows the same patterns of distribution as the marker gene/sequence itself. We present PlasmidTron, which utilizes the phenotypic data normally available in bacterial population studies, such as antibiograms, virulence factors, or geographical information, to identify traits that are likely to be present on DNA that can randomly reassort across defined bacterial populations. It is also possible to use this methodology to associate unknown genes/sequences (e.g. plasmid backbones) with a specific molecular signature or marker (e.g. resistance gene presence or absence) using PlasmidTron. PlasmidTron uses a k-mer-based approach to identify reads associated with a phylogenetically unlinked phenotype. These reads are then assembled de novo to produce contigs in a fast and scalable-to-large manner. PlasmidTron is written in Python 3 and is available under the open source licence GNU GPL3 from https://github.com/sanger-pathogens/plasmidtron.