Genomic epidemiology of meticillin-resistant Staphylococcus aureus ST22 widespread in communities of the Gaza Strip, 2009.

BackgroundRemarkably high carriage prevalence of a community-associated meticillin-resistant Staphylococcus aureus (MRSA) strain of sequence type (ST) 22 in the Gaza strip was reported in 2012. This strain is linked to the pandemic hospital-associated EMRSA-15. The origin and evolutionary history of ST22 in Gaza communities and the genomic elements contributing to its widespread predominance are unknown. Methods: We generated high-quality draft genomes of 61 ST22 isolates from Gaza communities and, along with 175 ST22 genomes from global sources, reconstructed the ST22 phylogeny and examined genotypes unique to the Gaza isolates. Results: The Gaza isolates do not exhibit a close relationship with hospital-associated ST22 isolates, but rather with a basal population from which EMRSA-15 emerged. There were two separate resistance acquisitions by the same MSSA lineage, followed by diversification of other genetic determinants. Nearly all isolates in the two distinct clades, one characterised by staphylococcal cassette chromosome mec (SCCmec) IVa and the other by SCCmec V and MSSA isolates, contain the toxic shock syndrome toxin-1 gene. Discussion: The genomic diversity of Gaza ST22 isolates is not consistent with recent emergence in the region. The results indicate that two divergent Gaza clones evolved separately from susceptible isolates. Researchers should not assume that isolates identified as ST22 in the community are examples of EMRSA-15 that have escaped their healthcare roots. Future surveillance of MRSA is essential to the understanding of ST22 evolutionary dynamics and to aid efforts to slow the further spread of this lineage.

Genomic Epidemiology of Penicillin-Nonsusceptible Pneumococci with Nonvaccine Serotypes Causing Invasive Disease in the United States.

Conjugate vaccination against seven pneumococcal serotypes (PCV7) reduced disease prevalence due to antibiotic-resistant strains throughout the 2000s. However, diseases caused by resistant nonvaccine type (NVT) strains increased. Some of these emerging strains were derived from vaccine types (VT) that had changed their capsule by recombination. The introduction of a vaccine targeting 13 serotypes (PCV13) in 2010 has led to concern that this scenario will repeat itself. We generated high-quality draft genomes from 265 isolates of NVT pneumococci not susceptible to penicillin (PNSP) in 2009 and compared them with the genomes of 581 isolates from 2012 to 2013 collected by the Active Bacterial Core surveillance (ABCs) of the Centers for Disease Control and Prevention (CDC). Of the seven sequence clusters (SCs) identified, three SCs fell into a single lineage associated with serogroup 23, which had an origin in 1908 as dated by coalescent analysis and included isolates with a divergent 23B capsule locus. Three other SCs represented relatively deep-branching lineages associated with serotypes 35B, 15A, and 15BC. In all cases, the resistant clones originated prior to 2010, indicating that PNSP are at present dominated by descendants of NVT clones present before vaccination. With one exception (15BC/ST3280), these SCs were related to clones identified by the Pneumococcal Molecular Epidemiology Network (PMEN). We conclude that postvaccine diversity in NVT PNSP between 2009 and 2013 was driven mainly by the persistence of preexisting strains rather than through de novo adaptation, with few cases of serotype switching. Future surveillance is essential for documenting the long-term dynamics and resistance of NVT PNSP.

Evaluation of calibration curve-based approaches to predict clinical inducers and noninducers of CYP3A4 with plated human hepatocytes.

Cytochrome P450 (P450) induction is often considered a liability in drug development. Using calibration curve-based approaches, we assessed the induction parameters R3 (a term indicating the amount of P450 induction in the liver, expressed as a ratio between 0 and 1), relative induction score, Cmax/EC50, and area under the curve (AUC)/F2 (the concentration causing 2-fold increase from baseline of the dose-response curve), derived from concentration-response curves of CYP3A4 mRNA and enzyme activity data in vitro, as predictors of CYP3A4 induction potential in vivo. Plated cryopreserved human hepatocytes from three donors were treated with 20 test compounds, including several clinical inducers and noninducers of CYP3A4. After the 2-day treatment, CYP3A4 mRNA levels and testosterone 6ß-hydroxylase activity were determined by real-time reverse transcription polymerase chain reaction and liquid chromatography-tandem mass spectrometry analysis, respectively. Our results demonstrated a strong and predictive relationship between the extent of midazolam AUC change in humans and the various parameters calculated from both CYP3A4 mRNA and enzyme activity. The relationships exhibited with non-midazolam in vivo probes, in aggregate, were unsatisfactory. In general, the models yielded better fits when unbound rather than total plasma Cmax was used to calculate the induction parameters, as evidenced by higher R(2) and lower root mean square error (RMSE) and geometric mean fold error. With midazolam, the R3 cut-off value of 0.9, as suggested by US Food and Drug Administration guidance, effectively categorized strong inducers but was less effective in classifying midrange or weak inducers. This study supports the use of calibration curves generated from in vitro mRNA induction response curves to predict CYP3A4 induction potential in human. With the caveat that most compounds evaluated here were not strong inhibitors of enzyme activity, testosterone 6ß-hydroxylase activity was also demonstrated to be a strong predictor of CYP3A4 induction potential in this assay model.

Inferring bacterial transmission dynamics using deep sequencing genomic surveillance data.

Identifying and interrupting transmission chains is important for controlling infectious diseases. One way to identify transmission pairs - two hosts in which infection was transmitted from one to the other - is using the variation of the pathogen within each single host (within-host variation). However, the role of such variation in transmission is understudied due to a lack of experimental and clinical datasets that capture pathogen diversity in both donor and recipient hosts. In this work, we assess the utility of deep-sequenced genomic surveillance (where genomic regions are sequenced hundreds to thousands of times) using a mouse transmission model involving controlled spread of the pathogenic bacterium Citrobacter rodentium from infected to naïve female animals. We observe that within-host single nucleotide variants (iSNVs) are maintained over multiple transmission steps and present a model for inferring the likelihood that a given pair of sequenced samples are linked by transmission. In this work we show that, beyond the presence and absence of within-host variants, differences arising in the relative abundance of iSNVs (allelic frequency) can infer transmission pairs more precisely. Our approach further highlights the critical role bottlenecks play in reserving the within-host diversity during transmission.

Rapid inference of antibiotic resistance and susceptibility by genomic neighbour typing.

Surveillance of drug-resistant bacteria is essential for healthcare providers to deliver effective empirical antibiotic therapy. However, traditional molecular epidemiology does not typically occur on a timescale that could affect patient treatment and outcomes. Here, we present a method called 'genomic neighbour typing' for inferring the phenotype of a bacterial sample by identifying its closest relatives in a database of genomes with metadata. We show that this technique can infer antibiotic susceptibility and resistance for both Streptococcus pneumoniae and Neisseria gonorrhoeae. We implemented this with rapid k-mer matching, which, when used on Oxford Nanopore MinION data, can run in real time. This resulted in the determination of resistance within 10 min (91% sensitivity and 100% specificity for S. pneumoniae and 81% sensitivity and 100% specificity for N. gonorrhoeae from isolates with a representative database) of starting sequencing, and within 4 h of sample collection (75% sensitivity and 100% specificity for S. pneumoniae) for clinical metagenomic sputum samples. This flexible approach has wide application for pathogen surveillance and may be used to greatly accelerate appropriate empirical antibiotic treatment.

Population genomics of pneumococcal carriage in Massachusetts children following introduction of PCV-13.

The 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large paediatric carriage sample collected from shortly after the introduction of PCV-7 to several years after the introduction of PCV-13, we investigate alterations in the composition of the pneumococcal population following the introduction of PCV-13, evaluating the extent to which the post-vaccination non-vaccine type (NVT) population mirrors that from prior to vaccine introduction and the effect of PCV-13 on vaccine type lineages. Draft genome assemblies from 736 newly sequenced and 616 previously published pneumococcal carriage isolates from children in Massachusetts between 2001 and 2014 were analysed. Isolates were classified into one of 22 sequence clusters (SCs) on the basis of their core genome sequence. We calculated the SC diversity for each sampling period as the probability that any two randomly drawn isolates from that period belong to different SCs. The sampling period immediately after the introduction of PCV-13 (2011) was found to have higher diversity than preceding (2007) or subsequent (2014) sampling periods {Simpson's D 2007: 0.915 [95 % confidence interval (CI) 0.901, 0.929]; 2011:  0.935 [0.927, 0.942]; 2014 :  0.912 [0.901, 0.923]}. Amongst NVT isolates, we found the distribution of SCs in 2011 to be significantly different from that in 2007 or 2014 (Fisher's exact test P=0.018, 0.0078), but did not find a difference comparing 2007 to 2014 (Fisher's exact test P=0.24), indicating greater similarity between samples separated by a longer time period than between samples from closer time periods. We also found changes in the accessory gene content of the NVT population between 2007 and 2011 to have been reduced by 2014. Amongst the new serotypes targeted by PCV-13, four were present in our sample. The proportion of our sample composed of PCV-13-only vaccine serotypes 19A, 6C and 7F decreased between 2007 and 2014, but no such reduction was seen for serotype 3. We did, however, observe differences in the genetic composition of the pre- and post-PCV-13 serotype 3 population. Our isolates were collected during discrete sampling periods from a small geographical area, which may limit the generalizability of our findings. Pneumococcal diversity increased immediately following the introduction of PCV-13, but subsequently returned to pre-vaccination levels. This is reflected in the distribution of NVT lineages, and, to a lesser extent, their accessory gene frequencies. As such, there may be a period during which the population is particularly disrupted by vaccination before returning to a more stable distribution. The persistence and shifting genetic composition of serotype 3 is a concern and warrants further investigation.

Diversification of bacterial genome content through distinct mechanisms over different timescales.

Bacterial populations often consist of multiple co-circulating lineages. Determining how such population structures arise requires understanding what drives bacterial diversification. Using 616 systematically sampled genomes, we show that Streptococcus pneumoniae lineages are typically characterized by combinations of infrequently transferred stable genomic islands: those moving primarily through transformation, along with integrative and conjugative elements and phage-related chromosomal islands. The only lineage containing extensive unique sequence corresponds to a set of atypical unencapsulated isolates that may represent a distinct species. However, prophage content is highly variable even within lineages, suggesting frequent horizontal transmission that would necessitate rapidly diversifying anti-phage mechanisms to prevent these viruses sweeping through populations. Correspondingly, two loci encoding Type I restriction-modification systems able to change their specificity over short timescales through intragenomic recombination are ubiquitous across the collection. Hence short-term pneumococcal variation is characterized by movement of phage and intragenomic rearrangements, with the slower transfer of stable loci distinguishing lineages.

A multi-endpoint evaluation of cytochrome P450 1A2, 2B6 and 3A4 induction response in human hepatocyte cultures after treatment with ß-naphthoflavone, phenobarbital and rifampicin.

U.S. FDA and EMEA guidance recommend that the preferred in vitro model for cytochrome P450 induction testing is human hepatocytes coupled with acceptable inducers as controls. However, there are surprisingly few published studies characterizing this model system for dose and time-dependence response to model inducing compounds. The concentration-dependent response and time-course for the induction of CYP1A2, CYP2B6 and CYP3A4 by inducing agents ß-naphthoflavone, phenobarbital and rifampicin, respectively were examined in two or more donors using multiple end-points (mRNA, enzyme activity and Western blot analysis). Depending on the endpoint, exposure time for maximal response of CYP induction potential for the three enzymes ranged from 24 to 72 hours. Of the concentrations of BNF, PB and RIF tested, those which gave the maximal response were found to be 33 µM, > 2 mM and 10 µM, respectively.