Prevalence of avian-pathogenic Escherichia coli strain O1 genomic islands among extraintestinal and commensal E. coli isolates

Escherichia coli strains that cause disease outside the intestine are known as extraintestinal pathogenic E. coli (ExPEC) and include pathogens of humans and animals. Previously, the genome of avian-pathogenic E. coli (APEC) O1:K1:H7 strain O1, from ST95, was sequenced and compared to those of several other E. coli strains, identifying 43 genomic islands. Here, the genomic islands of APEC O1 were compared to those of other sequenced E. coli strains, and the distribution of 81 genes belonging to 12 APEC O1 genomic islands among 828 human and avian ExPEC and commensal E. coli isolates was determined. Multiple islands were highly prevalent among isolates belonging to the O1 and O18 serogroups within phylogenetic group B2, which are implicated in human neonatal meningitis. Because of the extensive genomic similarities between APEC O1 and other human ExPEC strains belonging to the ST95 phylogenetic lineage, its ability to cause disease in a rat model of sepsis and meningitis was assessed. Unlike other ST95 lineage strains, APEC O1 was unable to cause bacteremia or meningitis in the neonatal rat model and was significantly less virulent than uropathogenic E. coli (UPEC) CFT073 in a mouse sepsis model, despite carrying multiple neonatal meningitis E. coli (NMEC) virulence factors and belonging to the ST95 phylogenetic lineage. These results suggest that host adaptation or genome modifications have occurred either in APEC O1 or in highly virulent ExPEC isolates, resulting in differences in pathogenicity. Overall, the genomic islands examined provide targets for further discrimination of the different ExPEC subpathotypes, serogroups, phylogenetic types, and sequence types.     

Comprehensive molecular,genomic and phenotypic analysis of a major clone of Enterococcus faecalis MLST ST40

Background

Enterococcus faecalis is a multifaceted microorganism known to act as a beneficial intestinal commensal bacterium. It is also a dreaded nosocomial pathogen causing life-threatening infections in hospitalised patients. Isolates of a distinct MLST type ST40 represent the most frequent strain type of this species, distributed worldwide and originating from various sources (animal, human, environmental) and different conditions (colonisation/infection). Since enterococci are known to be highly recombinogenic we determined to analyse the microevolution and niche adaptation of this highly distributed clonal type.

Results

We compared a set of 42 ST40 isolates by assessing key molecular determinants, performing whole genome sequencing (WGS) and a number of phenotypic assays including resistance profiling, formation of biofilm and utilisation of carbon sources. We generated the first circular closed reference genome of an E. faecalis isolate D32 of animal origin and compared it with the genomes of other reference strains. D32 was used as a template for detailed WGS comparisons of high-quality draft genomes of 14 ST40 isolates. Genomic and phylogenetic analyses suggest a high level of similarity regarding the core genome, also demonstrated by similar carbon utilisation patterns. Distribution of known and putative virulence-associated genes did not differentiate between ST40 strains from a commensal and clinical background or an animal or human source. Further analyses of mobile genetic elements (MGE) revealed genomic diversity owed to: (1) a modularly structured pathogenicity island; (2) a site-specifically integrated and previously unknown genomic island of 138 kb in two strains putatively involved in exopolysaccharide synthesis; and (3) isolate-specific plasmid and phage patterns. Moreover, we used different cell-biological and animal experiments to compare the isolate D32 with a closely related ST40 endocarditis isolate whose draft genome sequence was also generated. D32 generally showed a greater capacity of adherence to human cell lines and an increased pathogenic potential in various animal models in combination with an even faster growth in vivo (not in vitro).

Conclusion

Molecular, genomic and phenotypic analysis of representative isolates of a major clone of E. faecalis MLST ST40 revealed new insights into the microbiology of a commensal bacterium which can turn into a conditional pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1367-x) contains supplementary material, which is available to authorized users.     

Evolution of Streptococcus pneumoniae and its close commensal relatives

Streptococcus pneumoniae is a member of the Mitis group of streptococci which, according to 16S rRNA-sequence based phylogenetic reconstruction, includes 12 species. While other species of this group are considered prototypes of commensal bacteria, S. pneumoniae is among the most frequent microbial killers worldwide. Population genetic analysis of 118 strains, supported by demonstration of a distinct cell wall carbohydrate structure and competence pheromone sequence signature, shows that S. pneumoniae is one of several hundred evolutionary lineages forming a cluster separate from Streptococcus oralis and Streptococcus infantis. The remaining lineages of this distinct cluster are commensals previously collectively referred to as Streptococcus mitis and each represent separate species by traditional taxonomic standard. Virulence genes including the operon for capsule polysaccharide synthesis and genes encoding IgA1 protease, pneumolysin, and autolysin were randomly distributed among S. mitis lineages. Estimates of the evolutionary age of the lineages, the identical location of remnants of virulence genes in the genomes of commensal strains, the pattern of genome reductions, and the proportion of unique genes and their origin support the model that the entire cluster of S. pneumoniae, S. pseudopneumoniae, and S. mitis lineages evolved from pneumococcus-like bacteria presumably pathogenic to the common immediate ancestor of hominoids. During their adaptation to a commensal life style, most of the lineages gradually lost the majority of genes determining virulence and became genetically distinct due to sexual isolation in their respective hosts.     

Comparison of 61 Sequenced Escherichia coli Genomes

Escherichia coli is an important component of the biosphere and is an ideal model for studies of processes involved in bacterial genome evolution. Sixty-one publically available E. coli and Shigella spp. sequenced genomes are compared, using basic methods to produce phylogenetic and proteomics trees, and to identify the pan- and core genomes of this set of sequenced strains. A hierarchical clustering of variable genes allowed clear separation of the strains into clusters, including known pathotypes; clinically relevant serotypes can also be resolved in this way. In contrast, when in silico MLST was performed, many of the various strains appear jumbled and less well resolved. The predicted pan-genome comprises 15,741 gene families, and only 993 (6%) of the families are represented in every genome, comprising the core genome. The variable or ‘accessory’ genes thus make up more than 90% of the pan-genome and about 80% of a typical genome; some of these variable genes tend to be co-localized on genomic islands. The diversity within the species E. coli, and the overlap in gene content between this and related species, suggests a continuum rather than sharp species borders in this group of Enterobacteriaceae.     

Genome-based analysis of virulence genes in a non-biofilm-forming Staphylococcus epidermidis strain (ATCC 12228)

Staphylococcus epidermidis strains are diverse in their pathogenicity; some are invasive and cause serious nosocomial infections, whereas others are non-pathogenic commensal organisms. To analyse the implications of different virulence factors in Staphylococcus epidermidis infections, the complete genome of Staphylococcus epidermidis strain ATCC 12228, a non-biofilm forming, non-infection associated strain used for detection of residual antibiotics in food products, was sequenced. This strain showed low virulence by mouse and rat experimental infections. The genome consists of a single 2499 279 bp chromosome and six plasmids. The chromosomal G + C content is 32.1% and 2419 protein coding sequences (CDS) are predicted, among which 230 are putative novel genes. Compared to the virulence factors in Staphylococcus aureus, aside from delta-haemolysin and beta-haemolysin, other toxin genes were not found. In contrast, the majority of adhesin genes are intact in ATCC 12228. Most strikingly, the ica operon coding for the enzymes synthesizing interbacterial cellular polysaccharide is missing in ATCC 12228 and rearrangements of adjacent genes are shown. No mec genes, IS256, IS257, were found in ATCC 12228. It is suggested that the absence of the ica operon is a genetic marker in commensal Staphylococcus epidermidis strains which are less likely to become invasive.     

Comparative genomics inferred two distinct populations of piscine pathogenic Streptococcus agalactiae,serotype Ia ST7 and serotype III ST283, in Thailand and Vietnam

The genomes of Streptococcus agalactiae (group B streptococcus; GBS) collected from diseased fish in Thailand and Vietnam over a nine-year period (2008–2016) were sequenced and compared (n = 21). Based on capsular serotype and multilocus sequence typing (MLST), GBS isolates are divided into 2 groups comprised of i) serotype Ia; sequence type (ST)7 and ii) serotype III; ST283. Population structure inferred by core genome (cg)MLST and Bayesian clustering analysis also strongly indicated distribution of two GBS populations in both Thailand and Vietnam. Deep phylogenetic analysis implied by CRISPR array's spacer diversity was able to cluster GBS isolates according to their temporal and geographic origins, though ST7 has varying CRISPR1-spacer profiles when compared to ST283 strains. Based on overall genotypic features, Thai ST283 strains were closely related to the Singaporean ST283 strain causing foodborne illness in humans in 2015, thus, signifying zoonotic potential of this GBS population in the country.     

A Year of Infection in the Intensive Care Unit: Prospective Whole Genome Sequencing of Bacterial Clinical Isolates Reveals Cryptic Transmissions and Novel Microbiota

Bacterial whole genome sequencing holds promise as a disruptive technology in clinical microbiology, but it has not yet been applied systematically or comprehensively within a clinical context. Here, over the course of one year, we performed prospective collection and whole genome sequencing of nearly all bacterial isolates obtained from a tertiary care hospital’s intensive care units (ICUs). This unbiased collection of 1,229 bacterial genomes from 391 patients enables detailed exploration of several features of clinical pathogens. A sizable fraction of isolates identified as clinically relevant corresponded to previously undescribed species: 12% of isolates assigned a species-level classification by conventional methods actually qualified as distinct, novel genomospecies on the basis of genomic similarity. Pan-genome analysis of the most frequently encountered pathogens in the collection revealed substantial variation in pan-genome size (1,420 to 20,432 genes) and the rate of gene discovery (1 to 152 genes per isolate sequenced). Surprisingly, although potential nosocomial transmission of actively surveilled pathogens was rare, 8.7% of isolates belonged to genomically related clonal lineages that were present among multiple patients, usually with overlapping hospital admissions, and were associated with clinically significant infection in 62% of patients from which they were recovered. Multi-patient clonal lineages were particularly evident in the neonatal care unit, where seven separate Staphylococcus epidermidis clonal lineages were identified, including one lineage associated with bacteremia in 5/9 neonates. Our study highlights key differences in the information made available by conventional microbiological practices versus whole genome sequencing, and motivates the further integration of microbial genome sequencing into routine clinical care.     

HBV Genotypic Variability in Cuba

The genetic diversity of HBV in human population is often a reflection of its genetic admixture. The aim of this study was to explore the genotypic diversity of HBV in Cuba. The S genomic region of Cuban HBV isolates was sequenced and for selected isolates the complete genome or precore-core sequence was analyzed. The most frequent genotype was A (167/250, 67%), mainly A2 (149, 60%) but also A1 and one A4. A total of 77 isolates were classified as genotype D (31%), with co-circulation of several subgenotypes (56 D4, 2 D1, 5 D2, 7 D3/6 and 7 D7). Three isolates belonged to genotype E, two to H and one to B3. Complete genome sequence analysis of selected isolates confirmed the phylogenetic analysis performed with the S region. Mutations or polymorphisms in precore region were more common among genotype D compared to genotype A isolates. The HBV genotypic distribution in this Caribbean island correlates with the Y lineage genetic background of the population, where a European and African origin prevails. HBV genotypes E, B3 and H isolates might represent more recent introductions.     

Genomic characterization of two Staphylococcus epidermidis bacteriophages with anti-biofilm potential

ABSTRACT: BACKGROUND: Staphylococcus epidermidis is a commensal bacterium but can colonize the hospital environment due to its ability to form biofilms favouring adhesion to host tissues, medical devices and increasing resistance to antibiotics. In this context, the use of phages to destroy biofilms is an interesting alternative. RESULTS: The complete genomes of two Staphylococcus epidermidis bacteriophages, vB_SepiSphiIPLA5 and vB_SepiS-phiIPLA7, have been analyzed. Their genomes are 43,581 bp and 42,123 bp, and contain 67 and 59 orfs. Bioinformatic analyses enabled the assignment of putative functions to 36 and 29 gene products, respectively, including DNA packaging and morphogenetic proteins, lysis components, and proteins necessary for DNA recombination, regulation, modification and replication. A point mutation in vB_SepiS-phiIPLA5 lysogeny control-associated genes explained its strictly lytic behaviour. Comparative analysis of phi- IPLA5 and phi-IPLA7 genome structure resembled those of S. epidermidis phiPH15 and phiCNPH82 phages. A mosaic structure of S. epidermidis prophage genomes was revealed by PCR analysis of three marker genes (integrase, major head protein and holin). Using these genes, high prevalence (73%) of phage DNA in a representative S. epidermidis strain collection consisting of 60 isolates from women with mastitis and healthy women was determined. Putative pectin lyase-like domains detected in virion-associated proteins of both phages could be involved in exopolysaccharide (EPS) depolymerization, as evidenced by both the presence of a clear halo surrounding the phage lysis zone and the phage-mediated biofilm degradation. CONCLUSIONS: Staphylococcus epidermidis bacteriophages, vB_SepiS-phiIPLA5 and vB_SepiS-phiIPLA7, have a mosaic structure similar to other widespread S. epidermidis prophages. Virions of these phages are provided of pectin lyase-like domains, which may be regarded as promising anti-biofilm tools.     

Phylogenomics of Enterococcus faecalis from wild birds: new insights into host-associated differences in core and accessory genomes of the species

Wild birds have been suggested to be reservoirs of antimicrobial resistant and/or pathogenic Enterococcus faecalis (Efs) strains, but the scarcity of studies and available sequences limit our understanding of the population structure of the species in these hosts. Here, we analysed the clonal and plasmid diversity of 97 Efs isolates from wild migratory birds. We found a high diversity, with most sequence types (STs) being firstly described here, while others were found in other hosts including some predominant in poultry. We found that pheromone-responsive plasmids predominate in wild bird Efs while 35% of the isolates entirely lack plasmids. Then, to better understand the ecology of the species, the whole genome of fivestrains with known STs (ST82, ST170, ST16 and ST55) were sequenced and compared with all the Efs genomes available in public databases. Using several methods to analyse core and accessory genomes (AccNET, PLACNET, hierBAPS and PANINI), we detected differences in the accessory genome of some lineages (e.g. ST82) demonstrating specific associations with birds. Conversely, the genomes of other Efs lineages exhibited divergence in core and accessory genomes, reflecting different adaptive trajectories in various hosts. This pangenome divergence, horizontal gene transfer events and occasional epidemic peaks could explain the population structure of the species.     

Rapid and high-throughput genotyping of Staphylococcus epidermidis isolates by automated multilocus variable-number of tandem repeats: a tool for real-time epidemiology

Fast and reliable genotyping methods allowing real-time epidemiology would be instrumental to discriminate Staphylococcus epidermidis isolates, in order to evaluate potential cross-infections or to follow genome content of infecting strains of this important opportunistic pathogen. We describe an automated multilocus variable-number tandem repeat-based assay (MLVA) for the rapid genotyping of S. epidermidis. Multiplex PCR amplifications using 6 primer pairs targeting gene-regions containing variable numbers of tandem repeats and the mecA gene are resolved by micro-capillary electrophoresis and automatically assessed by cluster analysis. This genotyping technique was evaluated for discriminatory power and reproducibility on 2 sequenced strains, on a collection of 21 strains previously characterized using genotyping reference methods and finally on 65 clinical isolates identified in two different institutions. All steps of this new procedure were developed to ensure rapid turn-around time and moderate costs. Our results suggest that this rapid approach is a valuable epidemiological tool to genotype S. epidermidis isolates in real-time. The rapid analysis of a limited number of evolutionary markers showed a power of discrimination similar to that of pulse-field gel electrophoresis (PFGE) or multilocus sequence type (MLST). This type of rapid and high-throughput methodology opens the possibility to rapidly assess long-term nosocomial transmission or to characterize infecting strains in the general procedure of routine laboratories, in real-time.     

Short term evolution of a highly transmissible methicillin-resistant Staphylococcus aureus clone (ST228) in a tertiary care hospital

Staphylococcus aureus is recognized as one of the major human pathogens and is by far one of the most common nosocomial organisms. The genetic basis for the emergence of highly epidemic strains remains mysterious. Studying the microevolution of the different clones of S. aureus is essential for identifying the forces driving pathogen emergence and spread. The aim of the present study was to determine the genetic changes characterizing a lineage belonging to the South German clone (ST228) that spread over ten years in a tertiary care hospital in Switzerland. For this reason, we compared the whole genome of eight isolates recovered between 2001 and 2008 at the Lausanne hospital. The genetic comparison of these isolates revealed that their genomes are extremely closely related. Yet, a few more important genetic changes, such as the replacement of a plasmid, the loss of large fragments of DNA, or the insertion of transposases, were observed. These transfers of mobile genetic elements shaped the evolution of the ST228 lineage that spread within the Lausanne hospital. Nevertheless, although the strains analyzed differed in their dynamics, we have not been able to link a particular genetic element with spreading success. Finally, the present study showed that new sequencing technologies improve considerably the quality and quantity of information obtained for a single strain; but this information is still difficult to interpret and important investments are required for the technology to become accessible for routine investigations.     

Genomic and Genetic Diversity within the Pseudomonas fluorescens Complex

The Pseudomonas fluorescens complex includes Pseudomonas strains that have been taxonomically assigned to more than fifty different species, many of which have been described as plant growth-promoting rhizobacteria (PGPR) with potential applications in biocontrol and biofertilization. So far the phylogeny of this complex has been analyzed according to phenotypic traits, 16S rDNA, MLSA and inferred by whole-genome analysis. However, since most of the type strains have not been fully sequenced and new species are frequently described, correlation between taxonomy and phylogenomic analysis is missing. In recent years, the genomes of a large number of strains have been sequenced, showing important genomic heterogeneity and providing information suitable for genomic studies that are important to understand the genomic and genetic diversity shown by strains of this complex. Based on MLSA and several whole-genome sequence-based analyses of 93 sequenced strains, we have divided the P. fluorescens complex into eight phylogenomic groups that agree with previous works based on type strains. Digital DDH (dDDH) identified 69 species and 75 subspecies within the 93 genomes. The eight groups corresponded to clustering with a threshold of 31.8% dDDH, in full agreement with our MLSA. The Average Nucleotide Identity (ANI) approach showed inconsistencies regarding the assignment to species and to the eight groups. The small core genome of 1,334 CDSs and the large pan-genome of 30,848 CDSs, show the large diversity and genetic heterogeneity of the P. fluorescens complex. However, a low number of strains were enough to explain most of the CDSs diversity at core and strain-specific genomic fractions. Finally, the identification and analysis of group-specific genome and the screening for distinctive characters revealed a phylogenomic distribution of traits among the groups that provided insights into biocontrol and bioremediation applications as well as their role as PGPR.     

Archaeology and evolution of transfer RNA genes in the Escherichia coli genome

Transfer RNA genes tend to be presented in multiple copies in the genomes of most organisms, from bacteria to eukaryotes. The evolution and genomic structure of tRNA genes has been a somewhat neglected area of molecular evolution. Escherichia coli, the first phylogenetic species for which more than two different strains have been sequenced, provides an invaluable framework to study the evolution of tRNA genes. In this work, a detailed analysis of the tRNA structure of the genomes of Escherichia coli strains K12, CFT073, and O157:H7, Shigella flexneri 2a 301, and Salmonella typhimurium LT2 was carried out. A phylogenetic analysis of these organisms was completed, and an archaeological map depicting the main events in the evolution of tRNA genes was drawn. It is shown that duplications, deletions, and horizontal gene transfers are the main factors driving tRNA evolution in these genomes. On average, 0.64 tRNA insertions/duplications occur every million years (Myr) per genome per lineage, while deletions occur at the slower rate of 0.30 per million years per genome per lineage. This work provides a first genomic glance at the problem of tRNA evolution as a repetitive process, and the relationship of this mechanism to genome evolution and codon usage is discussed.     

The prevalence of tetracycline resistance genes in clinical and commensal strains of Enterococcus faecalis isolated in Gdańsk

The common use of antibiotics is responsible for selecting of drug resistance not only in pathogenic, clinical bacteria but also in commensal, not pathogenic strains which could cause the rapid dissemination of the resistance to these antibacterial agents. However, information regarding the antibiotic resistance of commensal bacteria is very scarce, and the data is based mostly on phenotypical research. Therefore the use of genotyping methods for detection of tetracycline resistance genes, in commensal and medical isolates of bacteria, is essential, for understanding the spread of antibiotic resistance. In this study 24 commensal and 27 clinical isolates of Enterococcus faecalis has been screened by PCR methods for tet(M), tet(S) genes and Tn916 and Tn5397 transpozons. Subsequently, the tet(M) gene amplicones were sequenced and phylogenetic analysis was performed. We have found that the prevalence of tet(S) gene varied significantly between commensal and clinical strains. Moreover, the frequency of transpozons in clinical isolates was much higher comparing to strains isolated from healthy individuals. The phylogenetic analysis did not show significant differences between clinical and commensal strains but it could suggest that the genetic similarity between these two groups could be favourable factor for broad range spread of tet(M) gene.     

Pan-Genomic Analysis Provides Insights into the Genomic Variation and Evolution of Salmonella Paratyphi A

Salmonella Paratyphi A (S. Paratyphi A) is a highly adapted, human-specific pathogen that causes paratyphoid fever. Cases of paratyphoid fever have recently been increasing, and the disease is becoming a major public health concern, especially in Eastern and Southern Asia. To investigate the genomic variation and evolution of S. Paratyphi A, a pan-genomic analysis was performed on five newly sequenced S. Paratyphi A strains and two other reference strains. A whole genome comparison revealed that the seven genomes are collinear and that their organization is highly conserved. The high rate of substitutions in part of the core genome indicates that there are frequent homologous recombination events. Based on the changes in the pan-genome size and cluster number (both in the core functional genes and core pseudogenes), it can be inferred that the sharply increasing number of pseudogene clusters may have strong correlation with the inactivation of functional genes, and indicates that the S. Paratyphi A genome is being degraded.     

The dominant Australian community-acquired methicillin-resistant Staphylococcus aureus clone ST93-IV [2B] is highly virulent and genetically distinct

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 has spread rapidly across North America, and CA-MRSA is also increasing in Australia. However, the dominant Australian CA-MRSA strain, ST93-IV [2B] appears distantly related to USA300 despite strikingly similar clinical and epidemiological profiles. Here, we compared the virulence of a recent Australian ST93 isolate (JKD6159) to other MRSA, including USA300, and found that JKD6159 was the most virulent in a mouse skin infection model. We fully sequenced the genome of JKD6159 and confirmed that JKD6159 is a distinct clone with 7616 single nucleotide polymorphisms (SNPs) distinguishing this strain from all other S. aureus genomes. Despite its high virulence there were surprisingly few virulence determinants. However, genes encoding α-hemolysin, Panton-Valentine leukocidin (PVL) and α-type phenol soluble modulins were present. Genome comparisons revealed 32 additional CDS in JKD6159 but none appeared to encode new virulence factors, suggesting that this clone's enhanced pathogenicity could lie within subtler genome changes, such as SNPs within regulatory genes. To investigate the role of accessory genome elements in CA-MRSA epidemiology, we next sequenced three additional Australian non-ST93 CA-MRSA strains and compared them with JKD6159, 19 completed S. aureus genomes and 59 additional S. aureus genomes for which unassembled genome sequence data was publicly available (82 genomes in total). These comparisons showed that despite its distinctive genotype, JKD6159 and other CA-MRSA clones (including USA300) share a conserved repertoire of three notable accessory elements (SSCmecIV, PVL prophage, and pMW2). This study demonstrates that the genetically distinct ST93 CA-MRSA from Australia is highly virulent. Our comparisons of geographically and genetically diverse CA-MRSA genomes suggest that apparent convergent evolution in CA-MRSA may be better explained by the rapid dissemination of a highly conserved accessory genome from a common source.     

Evolutionary dynamics of methicillin-resistant Staphylococcus aureus within a healthcare system

BackgroundIn the past decade, several countries have seen gradual replacement of endemic multi-resistant healthcare-associated methicillin-resistant Staphylococcus aureus (MRSA) with clones that are more susceptible to antibiotic treatment. One example is Singapore, where MRSA ST239, the dominant clone since molecular profiling of MRSA began in the mid-1980s, has been replaced by ST22 isolates belonging to EMRSA-15, a recently emerged pandemic lineage originating from Europe.ResultsWe investigated the population structure of MRSA in Singaporean hospitals spanning three decades, using whole genome sequencing. Applying Bayesian phylogenetic methods we report that prior to the introduction of ST22, the ST239 MRSA population in Singapore originated from multiple introductions from the surrounding region; it was frequently transferred within the healthcare system resulting in a heterogeneous hospital population. Following the introduction of ST22 around the beginning of the millennium, this clone spread rapidly through Singaporean hospitals, supplanting the endemic ST239 population. Coalescent analysis revealed that although the genetic diversity of ST239 initially decreased as ST22 became more dominant, from 2007 onwards the genetic diversity of ST239 began to increase once more, which was not associated with the emergence of a sub-clone of ST239. Comparative genomic analysis of the accessory genome of the extant ST239 population identified that the Arginine Catabolic Mobile Element arose multiple times, thereby introducing genes associated with enhanced skin colonization into this population.ConclusionsOur results clearly demonstrate that, alongside clinical practice and antibiotic usage, competition between clones also has an important role in driving the evolution of nosocomial pathogen populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0643-z) contains supplementary material, which is available to authorized users.     

Phylogenetic incongruence in E. coli O104: understanding the evolutionary relationships of emerging pathogens in the face of homologous recombination

Escherichia coli O104:H4 was identified as an emerging pathogen during the spring and summer of 2011 and was responsible for a widespread outbreak that resulted in the deaths of 50 people and sickened over 4075. Traditional phenotypic and genotypic assays, such as serotyping, pulsed field gel electrophoresis (PFGE), and multilocus sequence typing (MLST), permit identification and classification of bacterial pathogens, but cannot accurately resolve relationships among genotypically similar but pathotypically different isolates. To understand the evolutionary origins of E. coli O104:H4, we sequenced two strains isolated in Ontario, Canada. One was epidemiologically linked to the 2011 outbreak, and the second, unrelated isolate, was obtained in 2010. MLST analysis indicated that both isolates are of the same sequence type (ST678), but whole-genome sequencing revealed differences in chromosomal and plasmid content. Through comprehensive phylogenetic analysis of five O104:H4 ST678 genomes, we identified 167 genes in three gene clusters that have undergone homologous recombination with distantly related E. coli strains. These recombination events have resulted in unexpectedly high sequence diversity within the same sequence type. Failure to recognize or adjust for homologous recombination can result in phylogenetic incongruence. Understanding the extent of homologous recombination among different strains of the same sequence type may explain the pathotypic differences between the ON2010 and ON2011 strains and help shed new light on the emergence of this new pathogen.     

Detection and transformation of genome segments that differ within a coastal population of Vibrio cholerae strains

Vibrio cholerae is an autochthonous member of diverse aquatic ecosystems around the globe. Collectively, the genomes of environmental V. cholerae strains comprise a large repository of encoded functions which can be acquired by individual V. cholerae lineages through uptake and recombination. To characterize the genomic diversity of environmental V. cholerae, we used comparative genome hybridization to study 41 environmental strains isolated from diverse habitats along the central California coast, a region free of endemic cholera. These data were used to classify genes of the epidemic V. cholerae O1 sequenced strain N16961 as conserved, variably present, or absent from the isolates. For the most part, absent genes were restricted to large mobile elements and have known functions in pathogenesis. Conversely, genes present in some, but not all, California isolates were in smaller contiguous clusters and were less likely to be near genes with functions in DNA mobility. Two such clusters of variable genes encoding different selectable metabolic phenotypes (mannose and diglucosamine utilization) were transformed into the genomes of environmental isolates by chitin-dependent competence, indicating that this mechanism of general genetic exchange is conserved among V. cholerae. The transformed DNA had an average size of 22.7 kbp, demonstrating that natural competence can mediate the movement of large chromosome fragments. Thus, whether variable genes arise through the acquisition of new sequences by horizontal gene transfer or by the loss of preexisting DNA though deletion, natural transformation provides a mechanism by which V. cholerae clones can gain access to the V. cholerae pan-genome.