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.

561. Genomic Epidemiology of Methicillin-Resistant Staphylococcus aureus in Two Cohorts of High-Risk Military Trainees

Background

Methicillin-resistant Staphylococcus aureus (MRSA) skin and soft-tissue infections (SSTIs) are common among military recruits. Identifying which strains are responsible for SSTI and understanding the underlying transmission dynamics is critical to developing appropriate interventions for this high-risk population.

Methods

A cohort study of US Army Infantry trainees at Fort Benning, GA (June and September 2015). Participants from two training Companies were screened for colonization on multiple anatomic sites throughout the 14-week cycle as well as the time of clinical infection. MRSA+ samples were sequenced with Illumina HiSeq. Multi-locus sequence type (MLST) and virulence genes were identified in silico. Single nucleotide polymorphism (SNP) distances between soldiers’ bacteria were compared with assessing for potential transmission.

Results

Of 383 soldiers enrolled, 84 (22%) were colonized with MRSA during the study. Forty-two of 84 had a single positive colonization sample, of which 76% were from anatomical sites other than the nares (36% oropharyngeal, 26% perianal, 14% inguinal). Twelve trainees had MRSA SSTI during training (50% had colonization detected prior to or at infection). All were PFGE-type US300 (ST8) and were lukS/lukF-positive. SNP-based phylogenetic analyses and epidemiologic data indicate that most MRSA positives at baseline were due to unique importations from various community origins, suggesting that the ongoing MRSA epidemic is not due to a single endemic strain circulating on base. Following importation, extensive transmission then occurred, with multiple STs implicated. Transmission appeared restricted to within Companies, and predominantly within platoons.

Conclusion

Frequent colonization at baseline suggests a need for extensive MRSA screening and decolonization upon arrival to base, followed by ongoing infection control measures throughout training to prevent recolonization/infection. As multiple anatomical sites appear to play a role in transmission of MRSA, this may have important implications for screening protocols and control both in community and hospital-based settings.

Disclosures

All authors: No reported disclosures.

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

Remarkably 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.

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.

Populations of the pathogenic bacterium Streptococcus pneumoniae consist of distinct co-circulating lineages. Here, the authors show lineages are characterized by particular combinations of stable genomic islands, whereas prophage and restriction-modification systems vary over short timescales.