Santiago Castillo-Ramírez

A genomic portrait of the emergence, evolution and global spread of a methicillin resistant Staphylococcus aureus pandemic

The widespread use of antibiotics in association with high-density clinical care has driven the emergence of drug-resistant bacteria that are adapted to thrive in hospitalized patients. Of particular concern are globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clones that cause outbreaks and epidemics associated with health care. The most rapidly spreading and tenacious health-care-associated clone in Europe currently is EMRSA-15, which was first detected in the UK in the early 1990s and subsequently spread throughout Europe and beyond. Using phylogenomic methods to analyze the genome sequences for 193 S. aureus isolates, we were able to show that the current pandemic population of EMRSA-15 descends from a health-care-associated MRSA epidemic that spread throughout England in the 1980s, which had itself previously emerged from a primarily community-associated methicillin-sensitive population. The emergence of fluoroquinolone resistance in this EMRSA-15 subclone in the English Midlands during the mid-1980s appears to have played a key role in triggering pandemic spread, and occurred shortly after the first clinical trials of this drug. Genome-based coalescence analysis estimated that the population of this subclone over the last 20 yr has grown four times faster than its progenitor. Using comparative genomic analysis we identified the molecular genetic basis of 99.8% of the antimicrobial resistance phenotypes of the isolates, highlighting the potential of pathogen genome sequencing as a diagnostic tool. We document the genetic changes associated with adaptation to the hospital environment and with increasing drug resistance over time, and how MRSA evolution likely has been influenced by country-specific drug use regimens.

A very early-branching Staphylococcus aureus lineage lacking the carotenoid pigment staphyloxanthin

Here we discuss the evolution of the northern Australian Staphylococcus aureus isolate MSHR1132 genome. MSHR1132 belongs to the divergent clonal complex 75 lineage. The average nucleotide divergence between orthologous genes in MSHR1132 and typical S. aureus is approximately sevenfold greater than the maximum divergence observed in this species to date. MSHR1132 has a small accessory genome, which includes the well-characterized genomic islands, νSAα and νSaβ, suggesting that these elements were acquired well before the expansion of the typical S. aureus population. Other mobile elements show mosaic structure (the prophage φSa3) or evidence of recent acquisition from a typical S. aureus lineage (SCCmec, ICE6013 and plasmid pMSHR1132). There are two differences in gene repertoire compared with typical S. aureus that may be significant clues as to the genetic basis underlying the successful emergence of S. aureus as a pathogen. First, MSHR1132 lacks the genes for production of staphyloxanthin, the carotenoid pigment that confers upon S. aureus its characteristic golden color and protects against oxidative stress. The lack of pigment was demonstrated in 126 of 126 CC75 isolates. Second, a mobile clustered regularly interspaced short palindromic repeat (CRISPR) element is inserted into orfX of MSHR1132. Although common in other staphylococcal species, these elements are very rare within S. aureus and may impact accessory genome acquisition. The CRISPR spacer sequences reveal a history of attempted invasion by known S. aureus mobile elements. There is a case for the creation of a new taxon to accommodate this and related isolates.

The Impact of Recombination on dN/dS within Recently Emerged Bacterial Clones

The development of next-generation sequencing platforms is set to reveal an unprecedented level of detail on short-term molecular evolutionary processes in bacteria. Here we re-analyse genome-wide single nucleotide polymorphism (SNP) datasets for recently emerged clones of methicillin resistant Staphylococcus aureus (MRSA) and Clostridium difficile. We note a highly significant enrichment of synonymous SNPs in those genes which have been affected by recombination, i.e. those genes on mobile elements designated “non-core” (in the case of S. aureus), or those core genes which have been affected by homologous replacements (S. aureus and C. difficile). This observation suggests that the previously documented decrease in dN/dS over time in bacteria applies not only to genomes of differing levels of divergence overall, but also to horizontally acquired genes of differing levels of divergence within a single genome. We also consider the role of increased drift acting on recently emerged, highly specialised clones, and the impact of recombination on selection at linked sites. This work has implications for a wide range of genomic analyses.

Phylogeographic variation in recombination rates within a global clone of methicillin-resistant Staphylococcus aureus

BackgroundNext-generation sequencing (NGS) is a powerful tool for understanding both patterns of descent over time and space (phylogeography) and the molecular processes underpinning genome divergence in pathogenic bacteria. Here, we describe a synthesis between these perspectives by employing a recently developed Bayesian approach, BRATNextGen, for detecting recombination on an expanded NGS dataset of the globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clone ST239.ResultsThe data confirm strong geographical clustering at continental, national and city scales and demonstrate that the rate of recombination varies significantly between phylogeographic sub-groups representing independent introductions from Europe. These differences are most striking when mobile non-core genes are included, but remain apparent even when only considering the stable core genome. The monophyletic ST239 sub-group corresponding to isolates from South America shows heightened recombination, the sub-group predominantly from Asia shows an intermediate level, and a very low level of recombination is noted in a third sub-group representing a large collection from Turkey.ConclusionsWe show that the rapid global dissemination of a single pathogenic bacterial clone results in local variation in measured recombination rates. Possible explanatory variables include the size and time since emergence of each defined sub-population (as determined by the sampling frame), variation in transmission dynamics due to host movement, and changes in the bacterial genome affecting the propensity for recombination.

Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe

ABSTRACT The implementation of routine whole-genome sequencing (WGS) promises to transform our ability to monitor the emergence and spread of bacterial pathogens. Here we combined WGS data from 308 invasive Staphylococcus aureus isolates corresponding to a pan-European population snapshot, with epidemiological and resistance data. Geospatial visualization of the data is made possible by a generic software tool designed for public health purposes that is available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc). Our analysis demonstrates that high-risk clones can be identified on the basis of population level properties such as clonal relatedness, abundance, and spatial structuring and by inferring virulence and resistance properties on the basis of gene content. We also show that in silico predictions of antibiotic resistance profiles are at least as reliable as phenotypic testing. We argue that this work provides a comprehensive road map illustrating the three vital components for future molecular epidemiological surveillance: (i) large-scale structured surveys, (ii) WGS, and (iii) community-oriented database infrastructure and analysis tools. IMPORTANCE The spread of antibiotic-resistant bacteria is a public health emergency of global concern, threatening medical intervention at every level of health care delivery. Several recent studies have demonstrated the promise of routine whole-genome sequencing (WGS) of bacterial pathogens for epidemiological surveillance, outbreak detection, and infection control. However, as this technology becomes more widely adopted, the key challenges of generating representative national and international data sets and the development of bioinformatic tools to manage and interpret the data become increasingly pertinent. This study provides a road map for the integration of WGS data into routine pathogen surveillance. We emphasize the importance of large-scale routine surveys to provide the population context for more targeted or localized investigation and the development of open-access bioinformatic tools to provide the means to combine and compare independently generated data with publicly available data sets. The spread of antibiotic-resistant bacteria is a public health emergency of global concern, threatening medical intervention at every level of health care delivery. Several recent studies have demonstrated the promise of routine whole-genome sequencing (WGS) of bacterial pathogens for epidemiological surveillance, outbreak detection, and infection control. However, as this technology becomes more widely adopted, the key challenges of generating representative national and international data sets and the development of bioinformatic tools to manage and interpret the data become increasingly pertinent. This study provides a road map for the integration of WGS data into routine pathogen surveillance. We emphasize the importance of large-scale routine surveys to provide the population context for more targeted or localized investigation and the development of open-access bioinformatic tools to provide the means to combine and compare independently generated data with publicly available data sets.

A Common Genomic Framework for a Diverse Assembly of Plasmids in the Symbiotic Nitrogen Fixing Bacteria

This work centres on the genomic comparisons of two closely-related nitrogen-fixing symbiotic bacteria, Rhizobium leguminosarum biovar viciae 3841 and Rhizobium etli CFN42. These strains maintain a stable genomic core that is also common to other rhizobia species plus a very variable and significant accessory component. The chromosomes are highly syntenic, whereas plasmids are related by fewer syntenic blocks and have mosaic structures. The pairs of plasmids p42f-pRL12, p42e-pRL11 and p42b-pRL9 as well large parts of p42c with pRL10 are shown to be similar, whereas the symbiotic plasmids (p42d and pRL10) are structurally unrelated and seem to follow distinct evolutionary paths. Even though purifying selection is acting on the whole genome, the accessory component is evolving more rapidly. This component is constituted largely for proteins for transport of diverse metabolites and elements of external origin. The present analysis allows us to conclude that a heterogeneous and quickly diversifying group of plasmids co-exists in a common genomic framework.

Rapid evolutionary change of common bean (Phaseolus vulgaris L) plastome, and the genomic diversification of legume chloroplasts

BackgroundFabaceae (legumes) is one of the largest families of flowering plants, and some members are important crops. In contrast to what we know about their great diversity or economic importance, our knowledge at the genomic level of chloroplast genomes (cpDNAs or plastomes) for these crops is limited.ResultsWe sequenced the complete genome of the common bean (Phaseolus vulgari s cv. Negro Jamapa) chloroplast. The plastome of P. vulgaris is a 150,285 bp circular molecule. It has gene content similar to that of other legume plastomes, but contains two pseudogenes, rpl 33 and rps 16. A distinct inversion occurred at the junction points of trn H-GUG/rpl 14 and rps 19/rps 8, as in adzuki bean [1]. These two pseudogenes and the inversion were confirmed in 10 varieties representing the two domestication centers of the bean. Genomic comparative analysis indicated that inversions generally occur in legume plastomes and the magnitude and localization of insertions/deletions (indels) also vary. The analysis of repeat sequences demonstrated that patterns and sequences of tandem repeats had an important impact on sequence diversification between legume plastomes and tandem repeats did not belong to dispersed repeats. Interestingly, P. vulgaris plastome had higher evolutionary rates of change on both genomic and gene levels than G. max, which could be the consequence of pressure from both mutation and natural selection.ConclusionLegume chloroplast genomes are widely diversified in gene content, gene order, indel structure, abundance and localization of repetitive sequences, intracellular sequence exchange and evolutionary rates. The P. vulgaris plastome is a rapidly evolving genome.

The Rhizobium etli RpoH1 and RpoH2 sigma factors are involved in different stress responses.

The physiological role and transcriptional expression of Rhizobium etli sigma factors rpoH1 and rpoH2 are reported in this work. Both rpoH1 and rpoH2 were able to complement the temperature-sensitive phenotype of an Escherichia coli rpoH mutant. The R. etli rpoH1 mutant was sensitive to heat shock, sodium hypochlorite and hydrogen peroxide, whereas the rpoH2 mutant was sensitive to NaCl and sucrose. The rpoH2 rpoH1 double mutant had increased sensitivity to heat shock and oxidative stress when compared with the rpoH1 single mutant. This suggests that in R. etli, RpoH1 is the main heat-shock sigma factor, but a more complete protective response could be achieved with the participation of RpoH2. Conversely, RpoH2 is involved in osmotic tolerance. In symbiosis with bean plants, the R. etli rpoH1 and rpoH2 rpoH1 mutants still elicited nodule formation, but exhibited reduced nitrogenase activity and bacterial viability in early and late symbiosis compared with nodules produced by rpoH2 mutants and wild-type strains. In addition, nodules formed by R. etli rpoH1 and rpoH2 rpoH1 mutants showed premature senescence. It was also determined that fixNf and fixKf expression was affected in rpoH1 mutants. Both rpoH genes were induced under microaerobic conditions and in the stationary growth phase, but not in response to heat shock. Analysis of the upstream region of rpoH1 revealed a sigma70 and a probable sigmaE promoter, whereas in rpoH2, one probable sigmaE-dependent promoter was detected. In conclusion, the two RpoH proteins operate under different stress conditions, RpoH1 in heat-shock and oxidative responses, and RpoH2 in osmotic tolerance.

Genetic Relationships of Vibrio parahaemolyticus Isolates from Clinical, Human Carrier, and Environmental Sources in Thailand, Determined by Multilocus Sequence Analysis

ABSTRACT Vibrio parahaemolyticus is a seafood-borne pathogenic bacterium that is a major cause of gastroenteritis worldwide. We investigated the genetic and evolutionary relationships of 101 V. parahaemolyticus isolates originating from clinical, human carrier, and various environmental and seafood production sources in Thailand using multilocus sequence analysis. The isolates were recovered from clinical samples (n = 15), healthy human carriers (n = 18), various types of fresh seafood (n = 18), frozen shrimp (n = 16), fresh-farmed shrimp tissue (n = 18), and shrimp farm water (n = 16). Phylogenetic analysis revealed a high degree of genetic diversity within the V. parahaemolyticus population, although isolates recovered from clinical samples and from farmed shrimp and water samples represented distinct clusters. The tight clustering of the clinical isolates suggests that disease-causing isolates are not a random sample of the environmental reservoir, although the source of infection remains unclear. Extensive serotypic diversity occurred among isolates representing the same sequence types and recovered from the same source at the same time. These findings suggest that the O- and K-antigen-encoding loci are subject to exceptionally high rates of recombination. There was also strong evidence of interspecies horizontal gene transfer and intragenic recombination involving the recA locus in a large proportion of isolates. As the majority of the intragenic recombinational exchanges involving recA occurred among clinical and carrier isolates, it is possible that the human intestinal tract serves as a potential reservoir of donor and recipient strains that is promoting horizontal DNA transfer, driving evolutionary change, and leading to the emergence of new, potentially pathogenic strains.

Horizontal gene transfer and diverse functional constrains within a common replication-partitioning system in Alphaproteobacteria: the repABC operon

BackgroundThe repABC plasmid family, which is extensively present within Alphaproteobacteria, and some secondary chromosomes of the Rhizobiales have the particular feature that all the elements involved in replication and partitioning reside within one transcriptional unit, the repABC operon. Given the functional interactions among the elements of the repABC operon, and the fact that they all reside in the same operon, a common evolutionary history would be expected if the entire operon had been horizontally transferred. Here, we tested whether there is a common evolutionary history within the repABC operon. We further examined different incompatibility groups in terms of their differentiation and degree of adaptation to their host.ResultsWe did not find a single evolutionary history within the repABC operon. Each protein had a particular phylogeny, horizontal gene transfer events of the individual genes within the operon were detected, and different functional constraints were found within and between the Rep proteins. When different repABC operons coexisted in the same genome, they were well differentiated from one another. Finally, we found different levels of adaptation to the host genome within and between repABC operons coexisting in the same species.ConclusionHorizontal gene transfer with conservation of the repABC operon structure provides a highly dynamic operon in which each member of this operon has its own evolutionary dynamics. In addition, it seems that different incompatibility groups present in the same species have different degrees of adaptation to their host genomes, in proportion to the amount of time the incompatibility group has coexisted with the host genome.