Microbial Genomics

Visualization is frequently used to aid our interpretation of complex datasets. Within microbial genomics, visualizing the relationships between multiple genomes as a tree provides a framework onto which associated data (geographical, temporal, phenotypic and epidemiological) are added to generate hypotheses and to explore the dynamics of the system under investigation. Selected static images are then used within publications to highlight the key findings to a wider audience. However, these images are a very inadequate way of exploring and interpreting the richness of the data. There is, therefore, a need for flexible, interactive software that presents the population genomic outputs and associated data in a user-friendly manner for a wide range of end users, from trained bioinformaticians to front-line epidemiologists and health workers. Here, we present Microreact, a web application for the easy visualization of datasets consisting of any combination of trees, geographical, tempor...

Background
Escherichia coli is a model prokaryote, an important pathogen, and a key organism for industrial biotechnology. E. coli W (ATCC 9637), one of four strains designated as safe for laboratory purposes, has not been sequenced. E. coli W is a fast-growing strain and is the only safe strain that can utilize sucrose as a carbon source. Lifecycle analysis has demonstrated that sucrose from sugarcane is a preferred carbon source for industrial bioprocesses.

Results
We have sequenced and annotated the genome of E. coli W. The chromosome is 4,900,968 bp and encodes 4,764 ORFs. Two plasmids, pRK1 (102,536 bp) and pRK2 (5,360 bp), are also present. W has unique features relative to other sequenced laboratory strains (K-12, B and Crooks): it has a larger genome and belongs to phylogroup B1 rather than A. W also grows on a much broader range of carbon sources than does K-12. A genome-scale reconstruction was developed and validated in order to interrogate metabolic properties.

Conclusions
The genome of W is more similar to commensal and pathogenic B1 strains than phylogroup A strains, and therefore has greater utility for comparative analyses with these strains. W should therefore be the strain of choice, or 'type strain' for group B1 comparative analyses. The genome annotation and tools created here are expected to allow further utilization and development of E. coli W as an industrial organism for sucrose-based bioprocesses. Refinements in our E. coli metabolic reconstruction allow it to more accurately define E. coli metabolism relative to previous models.

Grass pea (Lathyrus sativus L.) is widely cultivated for food and feed in some developing countries including Ethiopia. However, due to its overexaggerated neuro-lathyrism alkaloid causing paralysis of limbs, it failed to attract attention of the research community and is one of the most neglected orphan crops in the world. But, the crop is considered an insurance crop by resource-poor farmers due to its strong abiotic stress tolerance and ability to produce high yields when all other crops fail due to unfavorable environmental conditions. This study was aimed at screening rhizobial isolates of grass pea and evaluating their symbiotic nitrogen fixation efficiency and tolerance to abiotic stresses. Fifty rhizobial isolates collected from grass pea nodules were isolated, screened, and characterized based on standard microbiological methods. The rhizobial isolates showed diversity in nodulation, symbiotic nitrogen fixation, and nutrient utilization. The 16S rRNA gene sequencing of 14 rhizobial isolates showed that two of them were identified as Rhizobium leguminosarum and the remaining twelve as Rhizobium species. Based on their overall performance, strains AAUGR-9, AAUGR-11, and AAUGR-14 that performed top and identified as Rhizobium species were recommended for field trials. This study screened and identified effective and competitive rhizobial isolates enriched with high nitrogen-fixing and abiotic stress tolerant traits, which contributes much to the application of microbial inoculants as alternative to chemical fertilizers.

The enteric, pathogenic spirochaete Brachyspira pilosicoli colonizes and infects a variety of birds and mammals, including humans. However, there is a paucity of genomic data available for this organism. This study introduces 12 newly sequenced draft genome assemblies, boosting the cohort of examined isolates by fourfold and cataloguing the intraspecific genomic diversity of the organism more comprehensively. We used several in silico techniques to define a core genome of 1751 genes and qualitatively and quantitatively examined the intraspecific species boundary using phylogenetic analysis and average nucleotide identity, before contextualizing this diversity against other members of the genus Brachyspira . Our study revealed that an additional isolate that was unable to be species typed against any other Brachyspira lacked putative virulence factors present in all other isolates. Finally, we quantified that homologous recombination has as great an effect on the evolution of the cor...

The identification of glycoside hydrolases (GHs) for efficient polysaccharide deconstruction is essential for the development of biofuels. Here, we investigate the potential of sequential HMM-profile identification for the rapid and precise identification of the multi-domain architecture of GHs from various datasets. First, as a validation, we successfully reannotated >98% of the biochemically characterized enzymes listed on the CAZy database. Next, we analyzed the 43 million non-redundant sequences from the M5nr data and identified 322,068 unique GHs. Finally, we searched 129 assembled metagenomes retrieved from MG-RAST for environmental GHs and identified 160,790 additional enzymes. Although most identified sequences corresponded to single domain enzymes, many contained several domains, including known accessory domains and some domains never identified in association with GH. several sequences displayed multiple catalytic domains and few of these potential multi-activity proteins combined potentially synergistic domains. Finally, we produced and confirmed the biochemical activities of a GH5-GH10 cellulase-xylanase and a GH11-CE4 xylanase-esterase. Globally, this "gene to enzyme pipeline" provides a rationale for mining large datasets in order to identify new catalysts combining unique properties for the efficient deconstruction of polysaccharides. Glycoside Hydrolases (GHs) are Carbohydrate-Active Enzymes (CAZy) that catalyze the hydrolysis of the glyco-sidic linkage in polysaccharides (e.g., cellulose, chitin) and oligosaccharides (e.g., cellobiose, chitobiose) 1. GHs are found as single domain proteins (SDGHs) or associated with accessory domains such as carbohydrate binding modules (CBMs) within multi-domain GHs (MDGHs) 2. In MDGHs, CBMs enhance enzyme-substrate interaction by anchoring the catalytic domain to the substrate 3. The anchoring reduces diffusion from the substrate and locally increases the concentration of catalytic domains 4 , thus improving the overall polysaccharide degradation 3. GHs support essential processes for ecosystem function and for biotechnology. Among others, in land ecosystems , the deconstruction of plant biomass by microbial GHs is essential 5,6 , whereas the breakdown of chitin, from arthropods and fungi, is important in both marine 7,8 and terrestrial ecosystems 9-12. Next, in the gut of animals, microbial GHs target polysaccharides, supplement the lack of endogenous enzymes 13,14 and thus contribute to the processing of complex carbohydrates during digestion 15-18. Finally, GHs are essential for the biofuel industry, as plant based polysaccharides constitute a major source of sustainable and renewable material capable of providing liquid transportation fuel 19-22. Many GH-genes and proteins have been identified in a growing number of sequenced genomes and environmental samples thanks to the use of activity-driven screening 23,24 and bioinformatic annotation systems 12,16,18,25. The precise identification of GH-genes and proteins is essential in order to understand how microbes support key functions across ecosystems 17,25,26 and to identify new enzymes for biotechnological application 18,21,27. In order to identify new catalysts for biomass degradation, we examined the performance of sequential Hidden Markov Model (HMM) identifications 28 combined with publicly accessible HMM-profiles from the PFam database 29 , here referred to as the GeneHunt approach 2,30 , to detect GH-sequences and investigate their detailed architecture (i.e., the precise domain organization of MDGHs) 2. More precisely, we first validated the GeneHunt approach by re-annotating the biochemically characterized GHs listed on the CAZy database (as of June 2018) 1. As described for cellulases, xylanases, and chitinases 30 , we expected the PFam-based annotation to correctly identify most of the proteins from the major GH families, although rare and recently identified GH-families would display inconsistencies. Next, we identified GHs in the M5nr database (version 13.12.15) containing 43,098,145 non-redundant, mostly microbial, protein sequences 31. This collection of sequences derived from the major sequence database serves as the reference database for the MG-RAST annotation pipeline 31,32. Finally, we identified the detailed multi-domain architecture of GH proteins in assembled, publicly accessible, metagenomes from

Background: The deep-sea bacterium Photobacterium profundum is an established model for studying high pressure adaptation. In this paper we analyse the parental strain DB110 and the toxR mutant TW30 by massively parallel cDNA sequencing (RNA-seq). ToxR is a transmembrane DNA-binding protein first discovered in Vibrio cholerae, where it regulates a considerable number of genes involved in environmental adaptation and virulence. In P. profundum the abundance and activity of this protein is influenced by hydrostatic pressure and its role is related to the regulation of genes in a pressure-dependent manner. Results: To better characterize the ToxR regulon, we compared the expression profiles of wt and toxR strains in response to pressure changes. Our results revealed a complex expression pattern with a group of 22 genes having expression profiles similar to OmpH that is an outer membrane protein transcribed in response to high hydrostatic pressure. Moreover, RNA-seq allowed a deep characterization of the transcriptional landscape that led to the identification of 460 putative small RNA genes and the detection of 298 protein-coding genes previously unknown. We were also able to perform a genome-wide prediction of operon structure, transcription start and termination sites, revealing an unexpected high number of genes (992) with large 5 0 -UTRs, long enough to harbour cis-regulatory RNA structures, suggesting a correlation between intergenic region size and UTR length. Conclusion: This work led to a better understanding of high-pressure response in P. profundum. Furthermore, the high-resolution RNA-seq analysis revealed several unexpected features about transcriptional landscape and general mechanisms of controlling bacterial gene expression.

High-temperature environments (>70°C) contain diverse and abundant members of the crenarchaeal order Thermoproteales. However, a comprehensive study of the distribution and function of diverse members of this group across different habitat types has not been conducted. Consequently, the goals of this study were to determine the distribution of different Thermoproteales genera across geochemically distinct geothermal habitats of Yellowstone National Park, and to identify key functional attributes of major genera that correlate with environmental parameters. Curated sequence assemblies belonging to five genera were characterized in replicate samples of 11 high-temperature communities ranging in pH from 3 to 9. Thermocladium, Vulcanisaeta and Caldivirga spp. were the primary Thermoproteales populations present in low pH (pH < 5) habitats, whereas Thermoproteus populations were found in mildly-acidic (pH 5–6) sulfur sediments, and Pyrobaculum populations were confined to higher pH (pH>6) sulfur sediments and/or filamentous ‘streamer’ communities. Metabolic reconstruction and comparative genomics among assemblies show that these populations are primarily chemoorganotrophs that utilize different electron acceptors depending on geochemical conditions. The presence of potential CO2 fixation pathways in some Thermoproteales populations appears to be linked with NiFe hydrogenases, which combined with high levels
of H2 in many sulfidic systems, may provide the energy required to fix inorganic C.

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.

An outbreak of multi-drug resistant (MDR) tuberculosis (TB) has been reported on Daru Island, Papua New Guinea. Mycobacterium tuberculosis strains driving this outbreak and the temporal accrual of drug resistance mutations have not been described. Whole genome sequencing of 100 of 165 clinical isolates referred from Daru General Hospital to the Supranational reference laboratory, Brisbane, during 2012-2015 revealed that 95 belonged to a single modern Beijing sub-lineage strain. Molecular dating suggested acquisition of streptomycin and isoniazid resistance in the 1960s, with potentially enhanced virulence mediated by an mycP1 mutation. The Beijing sub-lineage strain demonstrated a high degree of co-resistance between isoniazid and ethionamide (80/95; 84.2 %) attributed to an inhA promoter mutation combined with inhA and ndh coding mutations. Multi-drug resistance, observed in 78/95 samples, emerged with the acquisition of a typical rpoB mutation together with a compensatory rpoC mut...

Oceans cover approximately 70% of the Earth's surface with an average depth of 3800 m and a pressure of 38 MPa, thus a large part of the biosphere is occupied by high pressure environments. Piezophilic (pressure-loving) organisms are adapted to deep-sea life and grow optimally at pressures higher than 0.1 MPa. To better understand high pressure adaptation from a genomic point of view three different Photobacterium profundum strains were compared. Using the sequenced piezophile P. profundum strain SS9 as a reference, microarray technology was used to identify the genomic regions missing in two other strains: a pressure adapted strain (named DSJ4) and a pressure-sensitive strain (named 3TCK). Finally, the transcriptome of SS9 grown under different pressure (28 MPa; 45 MPa) and temperature (4°C; 16°C) conditions was analyzed taking into consideration the differentially expressed genes belonging to the flexible gene pool.

Human tuberculosis is caused by members of the Mycobacterium tuberculosis Complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto as well as two lineages (L5 and L6) traditionally referred to as M. africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known on their genomic diversity, phylogeography and evolution. Here, we analyzed the genomes of 365 L5 and 326 L6 strains, plus five related genomes that had not been classified into any of the known MTBC lineages, isolated from patients from 21 African countries.Our population genomic and phylogeographical analyses show that the unclassified genomes belonged to a new group that we propose to name MTBC Lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences betwe...

Vehicular data collection applications are emerging as an appealing technology to monitor urban areas, where a high concentration of connected vehicles with onboard sensors is a near future scenario. In this context, smartphones are, on one side, effective enablers of Vehicle-to-Vehicle (V2V) and Vehicle-to-Infrastructure (V2I) applications and, on the other side, highly sophisticated sensing platforms. In this paper, we introduce an effective and efficient system, denoted as D4V, to disseminate vehicle-related information and sensed data using smartphones as V2I devices. D4V relies on a Peer-to-Peer (P2P) overlay scheme, denoted as Distributed Geographic Table (DGT), which unifies the concepts of physical and virtual neighborhoods in a scalable and robust infrastructure for application-level services. First, we investigate the discovery procedure of the DGT overlay network, through analytical and simulation results. Then, we present and discuss an extensive simulation-based performance evaluation (considering relevant performance indicators) of the D4V system, in a 4G wireless communication scenario. The simulation methodology combines DEUS (an application-level simulation tool for the study of large-scale systems) with ns-3 (a well-known network simulator, which takes into account lower layers), in order to provide a D4V proof-of-concept. The observed results show that D4V-based information sharing among vehicles allows to significantly reduce risks and nuisances (e.g., due to road defects and congestions).

Pseudomonas aeruginosa is a ubiquitous opportunistic pathogen with several clones being frequently associated with outbreaks in hospital settings. ST395 is among these so-called 'international' clones. We aimed here to define the biological features that could have helped the implantation and spread of the clone ST395 in hospital settings. The complete genome of a multidrug resistant index isolate (DHS01) of a large hospital outbreak was analysed. We identified DHS01-specific genetic elements, among which were identified those shared with a panel of six independent ST395 isolates responsible for outbreaks in other hospitals. DHS01 has the fifth largest chromosome of the species (7.1 Mbp), with most of its 1555 accessory genes borne by either genomic islands (GIs, n=48) or integrative and conjugative elements (ICEs, n=5). DHS01 is multidrug resistant mostly due to chromosomal mutations. It displayed signatures of adaptation to chronic infection in part due to the loss of a 131 kbp chromosomal fragment. Four GIs were specific to the clone ST395 and contained genes involved in metabolism (GI-4), in virulence (GI-6) and in resistance to copper (GI-7). GI-7 harboured an array of six copper transporters and was shared with nonpathogenic Pseudomonas sp. retrieved from copper-contaminated environments. Copper resistance was confirmed phenotypically in all other ST395 isolates and possibly accounted for the spreading capability of the clone in hospital outbreaks, where water networks have been incriminated. This suggests that genes transferred from copper-polluted environments may have favoured the implantation and spread of the international clone P. aeruginosa ST395 in hospital settings. DATA SUMMARY 1. This whole-genome sequencing project has been deposited at GenBank/EMBL/DDBJ; accession number: CP013993.1 (urlhttps://www.ncbi.nlm.nih.gov/nuccore/CP013993.1). 2. The collections of ST395 and non-ST395 isolates have been deposited at GenBank/EMBL/DDBJ; Bio-Projects: PRJNA379554 (urlhttps://www.ncbi.nlm.nih. gov/bioproject/?term=PRJNA379554) and PRJNA380885 (urlhttps://www.ncbi.nlm.nih.gov/bioproject/?term= PRJNA380885) (Table S2, available in the online Supplementary Material). 3. The 82 Pseudomonas aeruginosa isolates for which complete genomes were available from the National Center for Biotechnology Information at the time of the study (June 2017) were included in the study.

There have been two anthrax cases affecting people that played and/or made animal-skin drums in the UK during the last 10 years, with single fatal occurrences in Scotland in 2006 and London in 2008. Investigations by the Health Protection Agency (now Public Health England) employing multi-locus-variable number tandem repeat analysis had previously linked the clinical cases to spores associated with animal skins and drums the patients had been in contact with. In this study, whole-genome sequencing of 23 Bacillus anthracis isolates harvested during the investigations was performed. High-quality draft assemblies of these genomes provided greater characterization of the B. anthracis strains present and placed them all upon a new branch of the global phylogeny. Although closely related, the clinical isolates from the two events, and another isolated from a drum-skin-associated case in New York in 2006, were distinct from each other. Multiple distinct genotypes were found during both investigations, implying either multiple contamination events or a single heterogeneous contamination. One environmental isolate from the Scottish incident was more closely related to London isolates than to the other Scottish isolates. As B. anthracis of this subgroup was present at both geographically and temporally distinct events, it may be more widespread at the source of contamination. All isolates were distinct from currently characterized West African strains, despite this being the likely origin of the drums and hides, therefore adding to our knowledge of B. anthracis diversity in the region.

The electrically conductive pili (e-pili) of Geobactersulfurreducens have environmental and practical significance because they can facilitate electron transfer to insoluble Fe(III) oxides; to other microbial species; and through electrically conductive biofilms. E-pili conductivity has been attributed to the truncated PilA monomer, which permits tight packing of aromatic amino acids to form a conductive path along the length of e-pili. In order to better understand the evolution and distribution of e-pili in the microbial world, type IVa PilA proteins from various Gram-negative and Gram-positive bacteria were examined with a particular emphasis on Fe(III)-respiring bacteria. E-pilin genes are primarily restricted to a tight phylogenetic group in the order Desulfuromonadales. The downstream gene in all but one of the Desulfuromonadales that possess an e-pilin gene is a gene previously annotated as 'pilA-C' that has characteristics suggesting that it may encode an outermembrane protein. Other genes associated with pilin function are clustered with e-pilin and 'pilA-C' genes in the Desulfuromonadales. In contrast, in the few bacteria outside the Desulfuromonadales that contain e-pilin genes, the other genes required for pilin function may have been acquired through horizontal gene transfer. Of the 95 known Fe(III)-reducing microorganisms for which genomes are available, 80 % lack e-pilin genes, suggesting that e-pili are just one of several mechanisms involved in extracellular electron transport. These studies provide insight into where and when e-pili are likely to contribute to extracellular electron transport processes that are biogeochemically important and involved in bioenergy conversions.

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

Carbapenem resistance is a rapidly growing threat to our ability to treat refractory bacterial infections. To understand how carbapenem resistance is mobilized and spread between pathogens, it is important to study the genetic context of the underlying resistance mechanisms. In this study, the resistomes of six clinical carbapenem-resistant isolates of five different species-Acinetobacter baumannii, Escherichia coli, two Klebsiella pneumoniae, Proteus mirabilis and Pseudomonas aeruginosa were characterized using whole genome sequencing. All Enterobacteriaceae isolates and the A. baumannii isolate had acquired a large number of antimicrobial resistance genes (7-18 different genes per isolate), including the following encoding carbapenemases: bla KPC-2 , bla OXA-48 , bla OXA-72 , bla NDM-1 , bla NDM-7 and bla VIM-1. In addition, a novel version of bla SHV was discovered. Four new resistance plasmids were identified and their fully assembled sequences were verified using optical DNA mapping. Most of the resistance genes were co-localized on these and other plasmids, suggesting a risk for coselection. In contrast, five out of six carbapenemase genes were present on plasmids with no or few other resistance genes. The expected level of resistancebased on acquired resistance determinantswas concordant with measured levels in most cases. There were, however, several important discrepancies for four of the six isolates concerning multiple classes of antibiotics. In conclusion, our results further elucidate the diversity of carbapenemases, their mechanisms of horizontal transfer and possible patterns of co-selection. The study also emphasizes the difficulty of using whole genome sequencing for antimicrobial susceptibility testing of pathogens with complex genotypes.

Shiga-toxin-producing Escherichia coli (STEC) O157:H7 is a recently emerged zoonotic pathogen with considerable morbidity. Since the emergence of this serotype in the 1980s, research has focussed on unravelling the evolutionary events from the E. coli O55:H7 ancestor to the contemporaneous globally dispersed strains observed today. In this study, the genomes of over 1000 isolates from both human clinical cases and cattle, spanning the history of STEC O157:H7 in the UK, were sequenced. Phylogenetic analysis revealed the ancestry, key acquisition events and global context of the strains. Dated phylogenies estimated the time to evolution of the most recent common ancestor of the current circulating global clone to be 175 years ago. This event was followed by rapid diversification. We show the acquisition of specific virulence determinates has occurred relatively recently and coincides with its recent detection in the human population. We used clinical outcome data from 493 cases of STEC O157:H7 to assess the relative risk of severe disease including haemolytic uraemic syndrome from each of the defined clades in the population and show the dramatic effect Shiga toxin repertoire has on virulence. We describe two strain replacement events that have occurred in the cattle population in the UK over the last 30 years, one resulting in a highly virulent strain that has accounted for the majority of clinical cases in the UK over the last decade. There is a need to understand the selection pressures maintaining Shiga-toxin-encoding bacteriophages in the ruminant reservoir and the study affirms the requirement for close surveillance of this pathogen in both ruminant and human populations.

We report two cases of severe pneumonia due to clone ST93 methicillin-resistant Staphylococcus aureus (MRSA) presenting from a remote Australian Indigenous community within a 2-week period, and the utilization of whole genome sequences to determine whether these were part of an outbreak. S. aureus was isolated from 12 of 92 nasal swabs collected from 25 community households (including the two index households); one isolate was ST93. Three of five skin lesion S. aureus isolates obtained at the community were ST93. Whole genome sequencing of the ST93 isolates from this study and a further 20 ST93 isolates from the same region suggested that recent transmission and progression to disease had not taken place. The proximity in time and space of the two severe pneumonia cases is probably a reflection of the high burden of disease due to ST93 MRSA in this population where skin infections and household crowding are common.

Genome wide microsatellite maps shall support Phytophthora systematics through the development of reliable markers, enabling species discrimination and variability analyses. Whole genome sequences of 17 Phytophthora accessions belonging to 14 species were retrieved from GenBank and the genome-wide microsatellites in each species were mined. A total of 51,200 microsatellites, including dinucleotide to decanucleotide motifs, have been identified across all the species and each one was characterized for uniqueness and repeat number. The P. infestans T30-4 genome had the highest (6873) and P. multivora 3378 had the lowest number of microsatellites (1802). Dinucleotide motifs (63.6%) followed by trinucleotide motifs (30.1%) were most abundant in all the genome. From 14 species, 250 microsatellites which are unique for the respective genomes are detailed along with their primer combinations and product sizes. P. sojae had the highest number of unique microsatellite motifs. Genome wide microsatellite maps for all the 14 Phytophthora species including the chromosome, position, motif, repeat number, forward and reverse primer sequences and expected PCR product size, for every microsatellite are presented. Markers based on the unique microsatellites could be used to identify each species, whereas the ones common to all species could be used to identify the genetic variability. Furthermore, to confirm the results, pure cultures of P. capsici, P. nicotianae and P. palmivora were procured from the Phytophthora Repository, DNA was isolated and the unique markers were screened across the species. The characteristic markers developed have confirmed the genome analysis results.

Background: The deep-sea bacterium Photobacterium profundum is an established model for studying high pressure adaptation. In this paper we analyse the parental strain DB110 and the toxR mutant TW30 by massively parallel cDNA sequencing (RNA-seq). ToxR is a transmembrane DNA-binding protein first discovered in Vibrio cholerae, where it regulates a considerable number of genes involved in environmental adaptation and virulence. In P. profundum the abundance and activity of this protein is influenced by hydrostatic pressure and its role is related to the regulation of genes in a pressure-dependent manner. Results: To better characterize the ToxR regulon, we compared the expression profiles of wt and toxR strains in response to pressure changes. Our results revealed a complex expression pattern with a group of 22 genes having expression profiles similar to OmpH that is an outer membrane protein transcribed in response to high hydrostatic pressure. Moreover, RNA-seq allowed a deep characterization of the transcriptional landscape that led to the identification of 460 putative small RNA genes and the detection of 298 protein-coding genes previously unknown. We were also able to perform a genome-wide prediction of operon structure, transcription start and termination sites, revealing an unexpected high number of genes (992) with large 5 0 -UTRs, long enough to harbour cis-regulatory RNA structures, suggesting a correlation between intergenic region size and UTR length. Conclusion: This work led to a better understanding of high-pressure response in P. profundum. Furthermore, the high-resolution RNA-seq analysis revealed several unexpected features about transcriptional landscape and general mechanisms of controlling bacterial gene expression.

Quorum sensing (QS) regulates the onset of bacterial social responses in function to cell density having an important impact in virulence. Autoinducer-2 (AI-2) is a signal that has the peculiarity of mediating both intra-and interspecies bacterial QS. We analyzed the diversity of all components of AI-2 QS across 44 complete genomes of Escherichia coli and Shigella strains. We used phylogenetic tools to study its evolution and determined the phenotypes of single-deletion mutants to predict phenotypes of natural strains. Our analysis revealed many likely adaptive polymorphisms both in gene content and in nucleotide sequence. We show that all natural strains possess the signal emitter (the luxS gene), but many lack a functional signal receptor (complete lsr operon) and the ability to regulate extracellular signal concentrations. This result is in striking contrast with the canonical species-specific QS systems where one often finds orphan receptors, without a cognate synthase, but not orphan emitters. Our analysis indicates that selection actively maintains a balanced polymorphism for the presence/absence of a functional lsr operon suggesting diversifying selection on the regulation of signal accumulation and recognition. These results can be explained either by niche-specific adaptation or by selection for a coercive behavior where signal-blind emitters benefit from forcing other individuals in the population to haste in cooperative behaviors.

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organizations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak an...

Pseudomonas aeruginosa is an extremely successful pathogen able to cause both acute and chronic infections in a range of hosts, utilizing a diverse arsenal of cell-associated and secreted virulence factors. A major cell-associated virulence factor, the Type IV pilus (T4P), is required for epithelial cell adherence and mediates a form of surface translocation termed twitching motility, which is necessary to establish a mature biofilm and actively expand these biofilms. P. aeruginosa twitching motility-mediated biofilm expansion is a coordinated, multicellular behaviour, allowing cells to rapidly colonize surfaces, including implanted medical devices. Although at least 44 proteins are known to be involved in the biogenesis, assembly and regulation of the T4P, with additional regulatory components and pathways implicated, it is unclear how these components and pathways interact to control these processes. In the current study, we used a global genomics-based random-mutagenesis technique, transposon directed insertion-site sequencing (TraDIS), coupled with a physical segregation approach, to identify all genes implicated in twitching motility-mediated biofilm expansion in P. aeruginosa. Our approach allowed identification of both known and novel genes, providing new insight into the complex molecular network that regulates this process in P. aeruginosa. Additionally, our data suggest that the flagellum-associated gene products have a differential effect on twitching motility, based on whether components are intra-or extracellular. Overall the success of our TraDIS approach supports the use of this global genomic technique for investigating virulence genes in bacterial pathogens.

Shigella sonnei is a major contributor to the global burden of diarrhoeal disease, generally associated with dysenteric diarrhoea in developed countries but also emerging in developing countries. The reason for the recent success of S. sonnei is unknown, but is likely catalysed by its ability to acquire resistance against multiple antimicrobials. Between 2011 and 2013, S. sonnei exhibiting resistance to fluoroquinolones, the first-line treatment recommended for shigellosis, emerged in Bhutan. Aiming to reconstruct the introduction and establishment of fluoroquinolone-resistant S. sonnei populations in Bhutan, we performed whole-genome sequencing on 71 S. sonnei samples isolated in Bhutan between 2011 and 2013.We found that these strains represented an expansion of a clade within the previously described lineage III, found specifically in Central Asia. Temporal phylogenetic reconstruction demonstrated that all of the sequenced Bhutanese S. sonnei diverged from a single ancestor that was introduced into Bhutan around 2006. Our data additionally predicted that fluoroquinolone resistance, conferred by mutations in gyrA and parC, arose prior to the introduction of the founder strain into Bhutan. Once established in Bhutan, these S. sonnei had access to a broad gene pool, as indicated by the acquisition of extendedspectrum b-lactamase-encoding plasmids and genes encoding type IV pili. The data presented here outline a model for the introduction and maintenance of fluoroquinolone-resistant S. sonnei in a new setting. Given the current circulation of fluoroquinolone-resistant S. sonnei in Asia, we speculate that this pattern of introduction is being recapitulated across the region and beyond.

Keywords: Salmonella enterica; Serotype typhi; Cell Invasion Protein SipB; Natural Product Library; Molecular Docking Background: The most dangerous issue in the healthcare arises is outbreak of the antibiotic-resistant bacteria worldwide. The randomness of Salmonella infections rely on its serotypes and the immune response of the host. Children up to age of 5 years and patients of older age are prone to the salmonella infections. The virulence genes of Salmonella encodes five different groups of Cell invasion proteins (CIPs), including Sip A, B, C, D and E. Upon contact with the target cell, SipB undergo ‘type III’ export from the bacterial cytoplasm and translocate into the cell membrane. Caspase-1 activity is essential for the cytotoxicity, and it has been proven that binding of SipB to caspase-1 induces macrophage apoptosis, so functional inhibition of caspase-1 blocks macrophage cytotoxicity. Methods: This study is to determine the components of the natural origin compounds (NO...

Most Staphylococcus aureus isolates can cause invasive disease given the right circumstances, but it is unknown if some isolates are more likely to cause severe infections than others. S. aureus bloodstream isolates from 120 patients with definite infective endocarditis and 121 with S. aureus bacteraemia without infective endocarditis underwent whole-genome sequencing. Genome-wide association analysis was performed using a variety of bioinformatics approaches including SNP analysis, accessory genome analysis and k-mer based analysis. Core and accessory genome analyses found no association with either of the two clinical groups. In this study, the genome sequences of S. aureus bloodstream isolates did not discriminate between bacteraemia and infective endocarditis. Based on our study and the current literature, it is not convincing that a specific S. aureus genotype is clearly associated to infective endocarditis in patients with S. aureus bacteraemia.

 Users may download and print one copy of any publication from the public portal for the purpose of private study or research.  You may not further distribute the material or use it for any profit-making activity or commercial gain  You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

Background: Despite the frequent isolation of Salmonella enterica sub. enterica serovars Derby and Mbandaka from livestock in the UK and USA little is known about the biological processes maintaining their prevalence. Statistics for Salmonella isolations from livestock production in the UK show that S. Derby is most commonly associated with pigs and turkeys and S. Mbandaka with cattle and chickens. Here we compare the first sequenced genomes of S. Derby and S. Mbandaka as a basis for further analysis of the potential host adaptations that contribute to their distinct host species distributions. Results: Comparative functional genomics using the RAST annotation system showed that predominantly mechanisms that relate to metabolite utilisation, in vivo and ex vivo persistence and pathogenesis distinguish S. Derby from S. Mbandaka. Alignment of the genome nucleotide sequences of S. Derby D1 and D2 and S. Mbandaka M1 and M2 with Salmonella pathogenicity islands (SPI) identified unique complements of genes associated with host adaptation. We also describe a new genomic island with a putative role in pathogenesis, SPI-23. SPI-23 is present in several S. enterica serovars, including S. Agona, S. Dublin and S. Gallinarum, it is absent in its entirety from S. Mbandaka. Conclusions: We discovered a new 37 Kb genomic island, SPI-23, in the chromosome sequence of S. Derby, encoding 42 ORFS, ten of which are putative TTSS effector proteins. We infer from full-genome synonymous SNP analysis that these two serovars diverged, between 182kya and 625kya coinciding with the divergence of domestic pigs. The differences between the genomes of these serovars suggest they have been exposed to different stresses including, phage, transposons and prolonged externalisation. The two serovars possess distinct complements of metabolic genes; many of which cluster into pathways for catabolism of carbon sources.

The lipopolysaccharide (O) and flagellar (H) surface antigens of Escherichia coli are targets for serotyping that have traditionally been used to identify pathogenic lineages. These surface antigens are important for the survival of E. coli within mammalian hosts. However, traditional serotyping has several limitations, and public health reference laboratories are increasingly moving towards whole genome sequencing (WGS) to characterize bacterial isolates. Here we present a method to rapidly and accurately serotype E. coli isolates from raw, short read WGS data. Our approach bypasses the need for de novo genome assembly by directly screening WGS reads against a curated database of alleles linked to known and novel E. coli O-groups and H-types (the EcOH database) using the software package SRST2. We validated the approach by comparing in silico results for 197 enteropathogenic E. coli isolates with those obtained by serological phenotyping in an independent laboratory. We then demonstrated the utility of our method to characterize isolates in public health and clinical settings, and to explore the genetic diversity of >1500 E. coli genomes from multiple sources. Importantly, we showed that transfer of O-and H-antigen loci between E. coli chromosomal backbones is common, with little evidence of constraints by host or pathotype, suggesting that E. coli 'strain space' may be virtually unlimited, even within specific pathotypes. Our findings show that serotyping is most useful when used in combination with strain genotyping to characterize microevolution events within an inferred population structure.

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

Food-borne outbreak investigation currently relies on the time-consuming and challenging bacterial isolation from food, to be able to link food-derived strains to more easily obtained isolates from infected people. When no food isolate can be obtained, the source of the outbreak cannot be unambiguously determined. Shotgun metagenomics approaches applied to the food samples could circumvent this need for isolation from the suspected source, but require downstream strain-level data analysis to be able to accurately link to the human isolate. Until now, this approach has not yet been applied outside research settings to analyse real food-borne outbreak samples. In September 2019, a Salmonella outbreak occurred in a hotel school in Bruges, Belgium, affecting over 200 students and teachers. Following standard procedures, the Belgian National Reference Center for human salmonellosis and the National Reference Laboratory for Salmonella in food and feed used conventional analysis based on i...

Sporulation is a survival strategy, adapted by bacterial cells in response to harsh environmental adversities. The adaptation potential differs between strains and the variations may arise from differences in gene regulation. Gene networks are a valuable way of studying such regulation processes and establishing associations between genes. We reconstructed and compared sporulation gene co-expression networks (GCNs) of the model laboratory strain Bacillus subtilis 168 and the food-borne industrial isolate Bacillus amyloliquefaciens. Transcriptome data obtained from samples of six stages during the sporulation process were used for network inference. Subsequently, a gene set enrichment analysis was performed to compare the reconstructed GCNs of B. subtilis 168 and B. amyloliquefaciens with respect to biological functions, which showed the enriched modules with coherent functional groups associated with sporulation. On basis of the GCNs and time-evolution of differentially expressed ge...

Background: The most dangerous issue in the healthcare arises is outbreak of the antibiotic resistant bacteria worldwide. The randomness of Salmonella infections rely on its serotypes and the immune response of the host. Children up to age of 5 years and patients of older age are prone to the salmonella infections. The virulence genes of Salmonella encodes five different groups of Cell invasion proteins (CIPs), including Sip A, B, C, D and E. Upon contact with the target cell, SipB undergo ‘type III’ export from the bacterial cytoplasm and translocate into the cellmembrane. Caspase-1 activity is essential for the cytotoxicity, and it has been proven that binding of SipB to caspase-1 induces macrophage apoptosis, so functional inhibition of caspase-1 blocks macrophage cytotoxicity. Methods: This study is to determine the components of the natural origin compounds (NOCs) have an anti-bacterial effect (ABE) and capable to prevent the humans from bacterial infection Salmonella enterica ...

Background: Drug-resistant tuberculosis is a high priority threat to global public health. There are still critical gaps in understanding how novel drug-resistant M. tuberculosis strains emerge and, once emergent, what drives the differential propagation of certain epidemiologically-successful strains over others. This study sought to describe the joint evolutionary and epidemiological histories of a novel multidrug-resistant M. tuberculosis strain recently identified in the capital city of the Republic of Moldova (MDR Ural/4.2). Methods: Using whole genome sequence data and Bayesian phylogenomic methods, we reconstruct the stepwise acquisition of drug-resistance mutations in the MDR Ural/4.2 strain, estimate its historical bacterial population size over time, and infer the migration history of this strain between Eastern European countries. Results: We infer that MDR Ural/4.2 likely evolved (via acquisition of rpoB S450L, which confers resistance to rifampin) in the early 1990s, du...

The global spread of Klebsiella pneumoniae producing Klebsiella pneumoniae carbapenemase (KPC) has been mainly associated with the dissemination of high-risk clones. In the last decade, hospital outbreaks involving KPC-producing K. pneumoniae have been predominantly attributed to isolates belonging to clonal group (CG) 258. However, results of recent epidemiological analysis indicate that KPC-producing sequence type (ST) 307, is emerging in different parts of the world and is a candidate to become a prevalent high-risk clone in the near future. Here we show that the ST307 genome encodes genetic features that may provide an advantage in adaptation to the hospital environment and the human host. Sequence analysis revealed novel plasmid-located virulence factors, including a cluster for glycogen synthesis. Glycogen production is considered to be one of the possible adaptive responses to long-term survival and growth in environments outside the host. Chromosomally-encoded virulence traits in the clone comprised fimbriae, an integrative conjugative element carrying the yersiniabactin siderophore, and two different capsular loci. Compared with the ST258 clone, capsulated ST307 isolates showed higher resistance to complement-mediated killing. The acquired genetic features identified in the genome of this new emerging clone may contribute to increased persistence of ST307 in the hospital environment and shed light on its potential epidemiological success.

Avian pathogenic Escherichia coli (APEC) cause widespread economic losses in poultry production and are potential zoonotic pathogens. Genome sequences of 95 APEC from commercial poultry operations in four Australian states that carried the class 1 integrase gene intI1, a proxy for multiple drug resistance (MDR), were characterized. Sequence types ST117 (22/95), ST350 (10/95), ST429 and ST57 (each 9/95), ST95 (8/95) and ST973 (7/95) dominated, while 24 STs were represented by one or two strains. FII and FIB repA genes were the predominant (each 93/95, 98 %) plasmid incompatibility groups identified, but those of B/O/K/Z (25/95, 26 %) and I1 (24/95, 25 %) were also identified frequently. Virulence-associated genes (VAGs) carried by ColV and ColBM virulence plasmids, including those encoding protectins [iss (91/95, 96 %), ompT (91/95, 96 %) and traT (90/95, 95 %)], iron-acquisition systems [sitA (88/95, 93 %), etsA (87/95, 92 %), iroN (84/95, 89 %) and iucD/iutA (84/95, 89 %)] and the putative avian haemolysin hylF (91/95, 96 %), featured prominently. Notably, mobile resistance genes conferring resistance to fluoroquinolones, colistin, extended-spectrum b-lactams and carbapenems were not detected in the genomes of these 95 APEC but carriage of the sulphonamide resistance gene, sul1 (59/95, 63 %), the trimethoprim resistance gene cassettes dfrA5 (48/95, 50 %) and dfrA1 (25/95, 27 %), the tetracycline resistance determinant tet(A) (51/95, 55 %) and the ampicillin resistance genes bla TEM-1A/B/C (48/95, 52 %) was common. IS26 (77/95, 81 %), an insertion element known to capture and mobilize a wide spectrum of antimicrobial resistance genes, was also frequently identified. These studies provide a baseline snapshot of drug-resistant APEC in Australia and their role in the carriage of ColV-like virulence plasmids.

The spread of antibiotic resistance within and between different bacterial populations is a major health problem on a global scale. The identification of genetic transformation in genomic data from Neisseria meningitidis, the meningococcus (Mc), and other bacteria is problematic, since similar or even identical alleles may be involved. A particular challenge in naturally transformable bacteria generally is to distinguish between common ancestry and true recombined sites in sampled genome sequences. Furthermore, the identification of recombination following experimental transformation of homologous alleles requires identifiable differences between donor and recipient, which in itself influences the propensity for homologous recombination (HR). This study identifies the distribution of HR events following intraspecies and interspecies Mc transformations of rpoB alleles encoding rifampicin resistance by whole-genome DNA sequencing and single nucleotide variant analysis. The HR events analysed were confined to the genomic region surrounding the single nucleotide genetic marker used for selection. An exponential length distribution of these recombined events was found, ranging from a few nucleotides to about 72 kb stretches. The lengths of imported sequences were on average found to be longer following experimental transformation of the recipient with genomic DNA from an intraspecies versus an interspecies donor (P<0.001). The recombination events were generally observed to be mosaic, with donor sequences interspersed with recipient sequence. Here, we present four models to explain these observations, by fragmentation of the transformed DNA, by interruptions of the recombination mechanism, by secondary recombination of endogenous self-DNA, or by repair/replication mechanisms.

The decreased costs of genome sequencing have increased capability to apply wholegenome sequence on epidemiological surveillance of zoonotic Campylobacter jejuni. However, knowledge about how genetically similar epidemiologically linked isolates can be is vital for correct application of this methodology. To address this issue in C. jejuni we investigated the spatial and temporal signals in the genomes of a major clonal complex and generalist lineage, ST-45 CC, by exploiting the population structure and genealogy and applying genome-wide association analysis of 340 isolates from across Europe collected over a wide time-range. The occurrence and strength of the geographical signal varied between sublineages and followed the clonal frame when present, while no evidence of a temporal signal was found. Certain sublineages of ST-45 CC formed discrete and genetically isolated clades to which geography and time had left only negligible traces in the genomes. We hypothesize that these ST-45 CC clades form globally expanded monomorphic clones possibly spread across Europe by migratory birds. In addition, we observed an incongruence between the genealogy of the strains and MLST typing, thereby challenging the existing clonal complex definition and use of a common MLST-based nomenclature for the ST-45 CC of C. jejuni.

Wolbachia are alpha-proteobacteria symbionts infecting a large range of arthropod species and two different families of nematodes. Interestingly, these endosymbionts are able to induce diverse phenotypes in their hosts: they are reproductive parasites within many arthropods, nutritional mutualists within some insects and obligate mutualists within their filarial nematode hosts. Defining Wolbachia “species” is controversial and so they are commonly classified into 16 different phylogenetic lineages, termed supergroups, named A to S. However, available genomic data remains limited and not representative of the full Wolbachia diversity; indeed, of the 24 complete genomes and 55 draft genomes of Wolbachia available to date, 84% belong to supergroups A and B, exclusively composed of Wolbachia from arthropods.For the current study, we took advantage of a recently developed DNA enrichment method to produce four complete genomes and two draft genomes of Wolbachia from filarial nematodes. Tw...

Cholera is a severe, waterborne diarrheal disease caused by toxin-producing strains of the bacterium Vibrio cholerae. Comparative genomics has revealed &quot;waves&quot; of cholera transmission and evolution, in which clones are successively replaced over decades and centuries. However, the extent of V. cholerae genetic diversity within an epidemic or even within an individual patient is poorly understood. Here, we characterized V. cholerae genomic diversity at a micro-epidemiological level within and between individual patients from Bangladesh and Haiti. To capture within-patient diversity, we isolated multiple (8 to 20) V. cholerae colonies from each of eight patients, sequenced their genomes and identified point mutations and gene gain/loss events. We found limited but detectable diversity at the level of point mutations within hosts (zero to three single nucleotide variants within each patient), and comparatively higher gene content variation within hosts (at least one gain/loss...

Carbapenem-resistant Klebsiella pneumoniae (CRKP) increasingly cause high-mortality outbreaks in hospital settings globally. Following a patient fatality at a hospital in Beijing due to a bla KPC-2-positive CRKP infection, close monitoring was put in place over the course of 14 months to characterize all bla KPC-2-positive CRKP in circulation in the hospital. Whole genome sequences were generated for 100 isolates from bla KPC-2-positive isolates from infected patients, carriers and the hospital environment. Phylogenetic analyses identified a closely related cluster of 82 sequence type 11 (ST11) isolates circulating in the hospital for at least a year prior to admission of the index patient. The majority of inferred transmissions for these isolates involved patients in intensive care units. Whilst the 82 ST11 isolates collected during the surveillance effort all had closely related chromosomes, we observed extensive diversity in their antimicrobial resistance (AMR) phenotypes. We were able to reconstruct the major genomic changes underpinning this variation in AMR profiles, including multiple gains and losses of entire plasmids and recombination events between plasmids, including transposition of bla KPC-2. We also identified specific cases where variation in plasmid copy number correlated with the level of phenotypic resistance to drugs, suggesting that the number of resistance elements carried by a strain may play a role in determining the level of AMR. Our findings highlight the epidemiological value of whole genome sequencing for investigating multi-drug-resistant hospital infections and illustrate that standard typing schemes cannot capture the extraordinarily fast genome evolution of CRKP isolates.

Shigella sonnei is a major contributor to the global burden of diarrhoeal disease, generally associated with dysenteric diarrhoea in developed countries but also emerging in developing countries. The reason for the recent success of S. sonnei is unknown, but is likely catalysed by its ability to acquire resistance against multiple antimicrobials. Between 2011 and 2013, S. sonnei exhibiting resistance to fluoroquinolones, the first-line treatment recommended for shigellosis, emerged in Bhutan. Aiming to reconstruct the introduction and establishment of fluoroquinolone-resistant S. sonnei populations in Bhutan, we performed whole-genome sequencing on 71 S. sonnei samples isolated in Bhutan between 2011 and 2013.We found that these strains represented an expansion of a clade within the previously described lineage III, found specifically in Central Asia. Temporal phylogenetic reconstruction demonstrated that all of the sequenced Bhutanese S. sonnei diverged from a single ancestor that ...

Stenotrophomonas maltophilia has evolved as one of the leading multidrug-resistant pathogens responsible for a variety of nosocomial infections especially in highly debilitated patients. As information on the genomic and intraspecies diversity of this clinically important pathogen is limited, we sequenced the whole genome of 27 clinical isolates from hospitalized patients. Phylogenomic analysis along with the genomes of type strains suggested that the clinical isolates are distributed over the Stenotrophomonas maltophilia complex (Smc) within the genus Stenotrophomonas. Further genome-based taxonomy coupled with the genomes of type strains of the genus Stenotrophomonas allowed us to identify five cryptic genomospecies, which are associated with the clinical isolates of S. maltophilia and are potentially novel species. These isolates share a very small core genome that implies a high level of genetic diversity within the isolates. Recombination analysis of core genomes revealed that ...

Serotype 1 Streptococcus pneumoniae is a leading cause of invasive pneumococcal disease (IPD) worldwide, with the highest burden in developing countries. We report the whole-genome sequencing analysis of 448 serotype 1 isolates from 27 countries worldwide (including 11 in Africa). The global serotype 1 population shows a strong phylogeographic structure at the continental level, and within Africa there is further region-specific structure. Our results demonstrate that regionspecific diversification within Africa has been driven by limited cross-region transfer events, genetic recombination and antimicrobial selective pressure. Clonal replacement of the dominant serotype 1 clones circulating within regions is uncommon; however, here we report on the accessory gene content that has contributed to a rare clonal replacement event of ST3081 with ST618 as the dominant cause of IPD in the Gambia.

Acinetobacter baumannii is nowadays a relevant nosocomial pathogen characterized by multidrug resistance (MDR) and concomitant difficulties to treat infections. OmpA is the most abundant A. baumannii outer membrane (OM) protein, and is involved in virulence, host cell recognition, biofilm formation, regulation of OM stability, permeability, and antibiotic resistance. OmpA members are two-domain proteins with an N-terminal eight-stranded β-barrel domain with four external loops (ELs) interacting with the environment, and a C-terminal periplasmic domain binding non-covalently to the peptidoglycan. Here, we combined data from genome sequencing, phylogenetic, and multilocus sequence analyses from 242 strains of the Acinetobacter calcoaceticus/Acinetobacter baumannii complex (ACB), 222 from A. baumannii, to explore ompA microevolutionary divergence. Five major ompA variant groups were identified (V1 to V5) comprising 50 different alleles coding for 29 different proteins. Polymorphisms we...