Editor's Pick Research Article

Genomic Analysis and Comparison of Two Gonorrhea Outbreaks

Xavier Didelot, Janina Dordel, Lilith K. Whittles, Caitlin Collins, Nicole Bilek, Cynthia J. Bishop, Peter J. White, David M. Aanensen, Julian Parkhill, Stephen D. Bentley, Brian G. Spratt, Simon R. Harris
Paul Stephen Keim, Editor
DOI: 10.1128/mBio.00525-16

ABSTRACT

Gonorrhea is a sexually transmitted disease causing growing concern, with a substantial increase in reported incidence over the past few years in the United Kingdom and rising levels of resistance to a wide range of antibiotics. Understanding its epidemiology is therefore of major biomedical importance, not only on a population scale but also at the level of direct transmission. However, the molecular typing techniques traditionally used for gonorrhea infections do not provide sufficient resolution to investigate such fine-scale patterns. Here we sequenced the genomes of 237 isolates from two local collections of isolates from Sheffield and London, each of which was resolved into a single type using traditional methods. The two data sets were selected to have different epidemiological properties: the Sheffield data were collected over 6 years from a predominantly heterosexual population, whereas the London data were gathered within half a year and strongly associated with men who have sex with men. Based on contact tracing information between individuals in Sheffield, we found that transmission is associated with a median time to most recent common ancestor of 3.4 months, with an upper bound of 8 months, which we used as a criterion to identify likely transmission links in both data sets. In London, we found that transmission happened predominantly between individuals of similar age, sexual orientation, and location and also with the same HIV serostatus, which may reflect serosorting and associated risk behaviors. Comparison of the two data sets suggests that the London epidemic involved about ten times more cases than the Sheffield outbreak.

IMPORTANCE The recent increases in gonorrhea incidence and antibiotic resistance are cause for public health concern. Successful intervention requires a better understanding of transmission patterns, which is not uncovered by traditional molecular epidemiology techniques. Here we studied two outbreaks that took place in Sheffield and London, United Kingdom. We show that whole-genome sequencing provides the resolution to investigate direct gonorrhea transmission between infected individuals. Combining genome sequencing with rich epidemiological information about infected individuals reveals the importance of several transmission routes and risk factors, which can be used to design better control measures.

INTRODUCTION

Gonorrhea is a sexually transmitted disease (STD) caused by the bacterium Neisseria gonorrhoeae. In the United Kingdom, gonorrhea is one of the most common bacterial STDs, and its reported incidence has markedly increased since 2008 in both men and women, reaching a total of 35,000 diagnosed cases in 2014 (1). Treatment with antimicrobials is usually successful, but increased resistance to many frontline antibiotics has recently been observed (2, 3). Implementation of effective control measures to mitigate the spread of gonorrhea is difficult due to a lack of understanding of the importance of the complex transmission routes and reservoirs (4). Traditional epidemiological studies are complicated by the facts that gonorrhea can be carried asymptomatically for months in about 10% of men and 50% of women (5) and in the United Kingdom infects disproportionately men who have sex with men (MSM) and young heterosexuals of black ethnicity and features significant geographical and temporal variations (6).

Molecular epidemiology approaches—for example, opa typing (7), N. gonorrhoeae multiantigen sequence typing (NG-MAST) (8), or multilocus sequence typing (MLST) (9)—have proved helpful to demonstrate that multiple strains often circulate simultaneously within a host population, to detect emerging, often resistant, clones, and to identify clusters of individuals infected with the same strain. Primarily, these typing schemes are useful to rule out transmission links between individuals who carry different types and do not have sufficient resolution to shed light on fine patterns of transmission between individuals carrying the same type.

In recent years, bacterial epidemiology has started to be transformed by the availability of fast, affordable whole-genome sequencing, which can help identify transmission links sometimes even at the level of direct transmission between individuals (10–12). The value of genomic data to investigate local outbreaks has been demonstrated for several bacterial pathogens, including Staphylococcus aureus (13, 14), Clostridium difficile (15, 16), and Mycobacterium tuberculosis (17, 18). In a study of gonococcal epidemiology, whole-genome sequencing has recently been applied to a collection of isolates from across the United States, showing important structuring of the pathogen population with both geography and sexual orientation of the hosts (19).

To test the usefulness of gonococcal genome sequencing to track transmission at the finer scale of outbreaks occurring within a city, we selected two local isolate collections with different epidemiological properties, both of which have been previously described without genomic data. The Sheffield isolate collection was assembled between 1995 and 2000 at the single genitourinary clinic in the city, from a mostly heterosexual population in which contact tracing was performed. Multiple strains of the pathogen have been found to circulate using opa typing (20, 21) and NG-MAST typing (22), and in both cases molecular typing was correlated with known sexual contact links. We applied whole-genome sequencing to 132 isolates with the most prevalent NG-MAST type, ST12.

The London collection consists of 2,045 isolates sampled between June and November 2004 from 13 major sexual health clinics throughout London, representing 54% of the 3,754 cases reported in London at that time (23, 24). Contact tracing information was not available, but rich metadata were recorded for each infected individual, including ethnic background, HIV status where known, postcode, and reporting clinic. All isolates were previously typed using NG-MAST, revealing the coexistence of a large number of strains, some of which were mostly MSM associated and others mostly heterosexual associated (23). Forty-five percent of infected individuals had one of 21 major strains, among which NG-MAST type ST225 was the most strongly MSM associated (with 92% of cases being MSM) and not geographically clustered within London (24). We applied whole-genome sequencing to 105 isolates of ST225.

Our aim was to assess the potential benefits that whole-genome sequencing data could provide for epidemiology within single NG-MAST types. A relatively recent history of transmission would be expected due to the shared type and the localized sampling frames. The epidemiological backgrounds of the two studies are very different, with the Sheffield data set coming from a single STD clinic in a city more than 10 times smaller than the London data set, which was collected at 13 clinics. Furthermore, the Sheffield isolates are heterosexual associated, whereas the London isolates are MSM associated. Consequently, this will also allow us to investigate how differences in variation within the two data sets reflect their underlying context.

RESULTS AND DISCUSSION

Genomic analysis of the Sheffield data set.We sequenced the genomes of 132 isolates from Sheffield dated from 1995 to 2000, all of which had the same genotype, ST12, as defined by NG-MAST (see Table S1 in the supplemental material). Few recombination events were detected using Gubbins (25), with a relative effect of recombination compared to mutation of r/m = 0.04, which was much lower than those of previous reports based on species-wide diversity (19, 26). This may reflect a difference in the effect of recombination when measured at different scales or diversity or between different lineages, as previously reported in other bacterial pathogens, such as Staphylococcus aureus (27), Clostridium difficile (15), or Streptococcus pneumoniae (28). The few detected recombination events occurred mainly in genes coding for outer membrane proteins that undergo antigenic variation and have previously been described as highly recombinant—for example, opa and pil genes (29). After both repetitive and recombinant regions were removed, a total of 156 variable sites were found to distinguish the genomes from each other. A maximum likelihood phylogeny was constructed using these data and showed a strong temporal signal based on the correlation between root-to-tip distances and isolation dates for all leaves in the phylogeny (R² = 0.49) (see Fig. S1 in the supplemental material).

Figure S1 

(A) Maximum likelihood tree reconstructed for the Sheffield genomes. (B) Temporal signal in the tree. The signal was strong enough to estimate the evolutionary rate shown by the solid line. Download Figure S1, PDF file, 0.01 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Table S1 

List of 132 genomes in the Sheffield data set. Table S1, PDF file, 0.03 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

We therefore applied the Bayesian evolutionary analysis software BEAST (30) to these data in order to reconstruct a timed phylogeny (Fig. 1). The molecular clock rate was estimated to be 1.41 × 10−6 single nucleotide polymorphisms (SNPs) per site per year with a 95% confidence interval (CI) of 1.15 × 10−6 to 1.70 × 10−6. This rate was equivalent to 3.05 (95% CI, 2.47 to 3.67) mutations per year across the genome, is in good agreement with previous estimates in N. gonorrhoeae (19, 26) and the related species Neisseria meningitidis (31), and was around the middle of the range of reported values for other bacterial pathogens (32, 33).

FIG 1 
FIG 1 

Timed phylogeny for the Sheffield data set reconstructed using BEAST. Known sexual contacts are indicated by uniquely colored circles and squares. (Inset) Intervals for the time to the last common ancestor (TMRCA) of each pair within a group of known sexual contacts.

Known sexual contacts were significantly clustered on the timed phylogeny (Fig. 1 [permutation test]; see Fig. S2 in the supplemental material [P < 10−4]). For the 25 pairs of known sexual contacts (see Table S1 in the supplemental material), the median time to the most recent common ancestor (TMRCA) of the isolated bacteria was 3.4 months, with an interquartile range from 2.3 to 5.1 months (Fig. 1, inset). When considering direct transmission between individuals sampled at roughly the same time, with no diversity being transmitted (due to a strong transmission bottleneck), the most recent common ancestor of the two sampled genomes would have existed within the pathogen population of the infector (34–36). The TMRCA is therefore a lower bound for the time from infection to sampling of the donor and an upper bound for the time from infection to sampling of the recipient, so that the mean TMRCA for pairs of contacts is expected to be approximately equal to the average duration of infection. Our result is in good agreement with previous estimates ranging from 2 to 6 months for the average duration of gonorrhea infection in several modeling studies (37–39). This duration is significantly longer than the few days reported for the incubation period in male experimental challenges (40), even accounting for the few additional days taken from symptom onset to care seeking (41). However, since the Sheffield population is predominantly heterosexual, the mean duration of infection is an average between men and women and takes into account the fact that a large fraction of women as well as a smaller proportion of men can remain asymptomatic for extended periods of time (42).

Figure S2 

Permutation test comparing the mean TMRCA for known sexual contacts (in red) with values that would be obtained if the sexual contact labels were permuted at random (histogram). Download Figure S2, PDF file, 0 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Since all pairs of known sexual contacts have a TMRCA of less than 8 months, we decided to use this value as the maximum threshold for the TMRCA between two individuals who have directly infected each other. This threshold should conservatively rule out transmission for pairs of genomes with a higher TMRCA, because all pairs of known sexual contacts fulfill this criterion, and yet not all of them are direct transmission pairs since the data include four triplets and one quadruplet of sexual contacts (Fig. 1, inset). Pairs of genomes sampled more than 8 months apart are unlikely to be transmission links. We compared all pairs of genomes sampled within 3 months of each other, which applies to all 25 pairs of known sexual contacts, and found that 24% (398/1,632) had a TMRCA of less than 8 months, including the known sexual contacts (see Fig. S3 in the supplemental material). Since all of these pairs were sampled within a short period of time in the same city and had the same NG-MAST type, a nongenomic analysis could not rule out transmission for any of them, but with the help of genomic data, we can confidently rule out transmission for the majority.

Figure S3 

Links between the Sheffield genomes, based on the 8-month maximum TMRCA criterion. Download Figure S3, PDF file, 0.1 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

To complement the phylogenetic analysis of the Sheffield data, we applied Outbreaker, which allows the direct reconstruction of a transmission tree representing transmission pathways within a sample, including the possibility of unsampled missing links in the transmission chains (43). Known sexual contacts clustered together on the reconstructed transmission tree, with at most three intermediates in the transmission chain between contacts (Fig. 2). When looking at all pairs of individuals in the Sheffield data set sampled within 3 months of each other, a strong correlation was found between their TMRCA in the BEAST tree and whether or not Outbreaker inferred them to be linked (see Fig. S4 in the supplemental material) (Kruskal-Wallis test [KWT], P < 10−15; mean TMRCA of 0.28 and 1.54 years for pairs linked and unlinked by Outbreaker, respectively), indicating a good agreement between the two approaches.

Figure S4 

Distribution of TMRCA between pairs of genomes in the Sheffield data set that have been assessed to be directly linked by Outbreaker (red) or not (blue). Download Figure S4, PDF file, 0 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 2 
FIG 2 

Transmission tree for the Sheffield data set reconstructed using Outbreaker. Cases are indicated by black dots, except for known sexual contacts, who are indicated using the same markers as in Fig. 1. Each case is aligned on the x axis with its reporting date, and the y axis is arbitrary. Black links between cases indicate inferred direct transmission, and gray links indicate indirect transmission through at least one unsampled case.

Genomic analysis of the London data set.We sequenced the genomes of 105 isolates sampled in London between June and November 2004, all of which were representatives of NG-MAST type ST225 (see Table S2 in the supplemental material). Significant recombination was detected using Gubbins (25), with r/m = 1.17, which was higher and in better agreement with previous reports (19, 26) than the value for the Sheffield data. Most of the recombination events occurred on deep branches and mainly affected outer membrane and pilus genes, repeat regions, and a prophage. A total of 167 variable sites were found after repetitive and recombinant regions were removed.

Table S2 

List of 105 genomes in the London data set. Table S2, PDF file, 0.04 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

The temporal signal was weak in this data set due to the short sampling span of only 6 months, so that the molecular clock rate could not be directly estimated with confidence from this data set. The distribution of root-to-tip distance versus sampling dates was, however, compatible with that expected under the rate estimated for the Sheffield data set (see Fig. S5 in the supplemental material). A timed tree was therefore produced using BEAST (30), but forcing the molecular clock rate to be equal to that estimated for the Sheffield data (Fig. 3). We also tried to apply Outbreaker to this data set as we did for the Sheffield data set, but the results were not biologically meaningful due to the shorter sampling interval for the London data, which is incompatible with Outbreaker’s assumption of a constant sampling density from beginning to end of an outbreak (43).

Figure S5 

(A) Maximum likelihood tree reconstructed for the London genomes. (B) Temporal signal in the tree. The signal was not strong enough to estimate the evolutionary rate, and instead the solid line represents the rate estimated in Fig. S1 based on the Sheffield data. Download Figure S5, PDF file, 0.01 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 3 
FIG 3 

Timed phylogeny for the London data set reconstructed using BEAST. Each isolate is annotated on the right-hand side as follows. First column: black for MSM, blue for heterosexual men, and red for heterosexual women. Second column: black for HIV negative and red for HIV positive. Third column: number of reported United Kingdom partners in the last 3 months, with black for zero, gray for one, dark blue for two to five, and light blue for six or more.

The analysis of the Sheffield data suggested a threshold of at most 8 months for the TMRCA compatible with direct transmission from one individual to another. The London data come from a predominantly MSM population, but the same threshold should be applicable since the asymptomatic frequency of anorectal cases among MSM is ~80% (44), similar to the ~75% in heterosexual females (42). We therefore applied the same criterion to all pairs of individuals in the London data who were sampled within 3 months of each other. We found that transmission was a possibility for 4% (165/4,251) of such pairs (see Fig. S6 in the supplemental material). Likely transmission links were significantly associated with shorter geographical distances between postcodes of residence in London (KWT, P = 6 × 10−3; mean, 10.3 versus 19.1 km), as would be expected if sexual partnerships tend to be geographically clustered. The age difference was significantly lower for linked individuals than for unlinked individuals (KWT, P = 1.3 × 10−4; mean, 7.5 versus 9.8 years), likely reflecting a tendency for sexual partners to be about the same age, as found by behavioral surveys such as Natsal-3 (the third of the National Surveys of Sexual Attitudes and Lifestyles) (45). Linked individuals reported the same sexual orientation more often than expected by chance for this sample of individuals (Fisher’s exact test [FET], P = 6 × 10−3). This significant similarity in location, age, and sexuality of linked individuals suggests that the transmission analysis, which was based only on pathogen genome similarity, successfully captured the correct transmission links. Among the 105 individuals, 24 were HIV positive, 47 were HIV negative, and 33 had unknown HIV status. Relative to the number of transmission links involving an HIV-positive individual and an HIV-negative individual, there was an excess of transmission links found between pairs of individuals who were both HIV positive (FET, P = 2.6 × 10−3), and simultaneously there was a dearth of transmission links when both individuals were HIV negative (FET, P = 0.03).

Figure S6 

Links between the London genomes, based on the 8-month maximum TMRCA criterion. Download Figure S6, PDF file, 0.01 MB.
Copyright © 2016 Didelot et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Among the 105 individuals in the London sample, the genomic analysis above suggested that 29 did not have any likely transmission link, 19 had one, 8 had two, 14 had three, and the remaining 35 had between four and eleven links. This number of likely transmission links for a given individual was not significantly associated with sexual orientation (KWT, P = 0.13; mean, 3.19 versus 1.28 links), reported sex abroad (KWT, P = 0.41; mean, 2.37 versus 3.41 links), or reported previous gonorrhea (KWT, P = 0.33; mean, 3.52 versus 2.89 links) but was significantly increased for individuals who were HIV positive (KWT, P = 1.6 × 10−4; mean, 5.25 versus 2.08 links), in accordance with the observation above of a higher number of links between HIV-positive individuals. This observation may indicate a role of HIV infection in susceptibility and transmissibility of gonorrhea (46, 47), may be caused by both HIV and gonorrhea infection being linked with the same high-risk behaviors, or both. It is also probable that HIV-positive MSM have on average more partners and engage in serosorting and riskier behaviors (48, 49). A significant correlation was found between the number of transmission links and the reported number of sexual partners in the United Kingdom in the last 3 months (Spearman’s rank correlation test, rho = 0.4; P = 2 × 10−5). This correlation was partly explained by a higher number of reported partners for HIV-positive individuals (KWT, P = 2.2 × 10−3; mean, 8.03 versus 2.61), but a correlation was also suggested between number of links and number of partners when analyzing HIV-positive and HIV-negative individuals separately (rho = 0.5, P = 0.01, and rho = 0.3, P = 0.04, respectively).

Comparison between the Sheffield and London data sets.Having analyzed the Sheffield and London data sets separately, we now turn to comparative analysis between them. The temporal sampling frames are different in the two data sets since sampling in Sheffield happened over 6 years, whereas the London data cover only 6 months. It is therefore not possible, for example, to compare the frequencies with which a putative transmission donor is found for cases in both data sets, because this frequency would be expected to be smaller in London, since donors would more often have been reported before sampling started. However, direct comparisons become possible by focusing on pairs of isolates sampled at approximately the same time, since the relationships between such pairs should not be affected by differences in sampling frame duration. Having reconstructed dated phylogenies for both data sets (Fig. 1 and 3), we can compare the distributions of TMRCA for all pairs of isolates sampled within 3 months of each other (Fig. 4). As previously noted, a larger proportion of pairs have a TMRCA under 8 months in Sheffield relative to London (24% versus 4%), and overall the London pairs have an average TMRCA of 4.4 years compared to 1.5 years for the Sheffield pairs, with the difference being highly significant (KWT, P < 1 × 10−15). More pairs of individuals in Sheffield are therefore likely to have infected each other compared to in London, and this result is robust to the choice of other TMRCA thresholds than 8 months for when to rule out transmission (Fig. 4).

FIG 4 
FIG 4 

Density histograms of the time to the most recent common ancestor for all pairs of cases sampled within 3 months of each other. The Sheffield data are shown in red when a sexual contact was reported and in green otherwise, whereas the London data are shown in blue.

This difference between the Sheffield and London data sets can be explained from a population genetics viewpoint by a pathogen effective population size (Ne) roughly three times higher in London relative to Sheffield (50). Using a metapopulation analogy, the effective population size of a pathogen can be shown to be proportional to the number of infected individuals and inversely proportional to the transmission rate (51–53). A recent study estimated that gonorrhea spreads two to three times faster in MSM than heterosexual networks (54), which is consistent with a lower proportion of asymptomatic cases in men compared to women and symptoms that are more likely to result in rapid care seeking (5, 42). This difference should contribute to a lower Ne, suggesting that the number of infected individuals in London at a given time may have been six to nine times higher than in Sheffield for the two specific lineages of gonorrhea studied here, namely, NG-MAST ST12 in Sheffield and ST225 in London. This estimation is in broad agreement with the similar number of reported cases in Sheffield and London, despite the sampling interval being about ten times longer in the former.

Conclusions.We have presented the first genomic epidemiology investigation of two large, localized outbreaks of gonorrhea—one including 132 isolates from Sheffield collected over 6 years among a mostly heterosexual population, and the other including 105 isolates from London gathered over half a year from a mainly MSM population. We showed that whole-genome sequencing can be used to predict person-to-person transmission events, which we combined with epidemiological information about the infected individuals to reveal patterns of transmission and infection risk factors. Importantly, both isolate collections had previously been assigned to single NG-MAST types, highlighting again the superior resolution provided by genomic data for molecular epidemiology. Pairs of cases detected in the same town, within 3 months of each other and carrying the same molecular type, would be considered to be “linked” using traditional molecular epidemiology, and yet we showed that genome sequencing ruled out transmission for the majority of such pairs in both outbreaks. Furthermore, the proportions of pairs of cases for which transmission were found to be likely using genomic comparisons were very different in Sheffield (24%) and in London (4%), and we estimated that the number of cases was six to nine times higher in London compared to Sheffield. Molecular typing using NG-MAST remains useful to identify cases that are part of the same local outbreak, but our results clearly show that a type does not represent a uniform epidemiological entity. As cost effectiveness and turnaround time of whole-genome sequencing continue to improve, it will become an increasingly important tool in the investigation and control of gonococcal outbreaks, but unlocking its full potential requires the simultaneous application of traditional epidemiological techniques such as the use of questionnaires or contact tracing.

MATERIALS AND METHODS

Bacterial isolates.From the previously described Sheffield collection (20–22), which was sampled between 1995 and 2000, we included 132 out of 140 isolates of the most prevalent NG-MAST type, ST12. Table S1 in the supplemental material contains the list of these isolates and associated metadata. From the previously described London collection (23, 24), which consists of 2,045 isolates sampled between June and November 2004, we included 105 out of 124 isolates from the most MSM-associated NG-MAST type, ST225. Table S2 contains the list of these isolates and associated metadata. Exclusion of 8 and 19 isolates from the Sheffield and London collections, respectively, was due to loss of samples, loss of associated metadata, or failure to grow or sequence.

Whole-genome sequencing.For each of the two data sets, the oldest available isolate was used for the production of a reference genome using 454. A 3-kb library was prepared and sequenced on a 1/4 plate run. For the London reference isolate, this produced 139.2 Mbp. De novo assembly was carried out using Newbler, resulting in 12 scaffolds with a total length of 2.14 Mbp. The sequencing of the Sheffield isolate produced 91.1 Mbp and 13 scaffolds with a total length of 2.14 Mbp. The length of these two reference genomes is in good agreement with previous reports of the lengths of N. gonorrhoeae genomes (19, 26, 55). Both assemblies were annotated using the RAST annotation pipeline (56) and manually curated. Prophages were detected using the online tool PHAST (57). Further mobile genetic elements were detected through genome comparisons with already published N. gonorrhoeae genomes and BLAST analysis. Repeats were detected using the EMBOSS tools einverted and equicktandem (58).

The genomes of all Sheffield isolates and all London isolates were sequenced on an Illumina HiSeq2000 with 100-bp paired-end reads. Paired-end reads were mapped against the corresponding reference genome using SMALTv0.7.5 (http://www.sanger.ac.uk/science/tools/smalt-0) with subsequent realignment around indels using GATKv1.5.9 (59). Single nucleotide polymorphisms (SNPs) were called as previously described (60). Recombination was detected using Gubbins (25), and recombinant SNPs were excluded. Furthermore, SNPs within mobile genetic elements as well as repetitive regions were also excluded. This resulted in alignments containing 156 and 167 SNPs for Sheffield and London, respectively.

Phylogenetic analysis.Maximum likelihood phylogenies were constructed using phyml (61). In the Sheffield tree, we included a single genome from London, and in the London tree, we included a single genome from Sheffield, to be used as outgroups in order to root the trees. Comparison of root-to-tip distances with isolation dates for each genome revealed a strong temporal signal in the Sheffield data (see Fig. S1 in the supplemental material) but not in London (see Fig. S5 in the supplemental material) due to a short sampling interval. The Bayesian evolutionary analysis software BEAST (30) was applied to both data sets, with the molecular clock of the London data set being forced equal to the rate estimated for the Sheffield data set. The resulting timed phylogenies are shown in Fig. 1 and 3. We used BEAST version 1.8.2 with the default HKY substitution model, the default coalescent model with constant population size, and a strict clock model with rate prior distribution Exponential(1) for the Sheffield data and fixed rate for the London data. For each of the two data sets, we performed four runs of 107 iterations, which were compared in Tracer (http://beast.bio.ed.ac.uk/tracer) to confirm convergence.

Transmission analysis.We sought a value for the threshold on the time to the most recent common ancestry of two genomes beyond which direct transmission between the two corresponding individuals can be discounted (15). Based on known sexual contacts in the Sheffield data set, we estimated that 8 months was an appropriate value for this threshold. We applied this criterion to all pairs of genomes sampled within 3 months of each other in both Sheffield and London to determine for each pair whether transmission was likely or not. Pairs of genomes for which transmission was likely are shown in Fig. S3 and S6 in the supplemental material for the Sheffield and London data sets, respectively. We also used Outbreaker, a Bayesian approach for reconstructing a transmission tree from dated genetic data (43). For the generation time, we used a discretized gamma distribution with a mean of 90 days and a standard deviation of 40 days, which has high variance reflecting our lack of knowledge of exact values (62). The Outbreaker output for the Sheffield data is shown in Fig. 2 and compared with the BEAST-based approach in Fig. S4 in the supplemental material. Outbreaker did not produce meaningful results for the London data due to the short sampling frame, which implies that the donor of many cases would have occurred and been reported before the sampling frame.

Nucleotide sequence accession numbers.Sequence data have been deposited in the European Nucleotide Archive. Accession numbers for the Illumina data from each isolate can be found in Tables S1 and S2 in the supplemental material.

FOOTNOTES

    • Received 22 March 2016
    • Accepted 17 May 2016
    • Published 28 June 2016

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Public Health England
    . 2015. Sexually transmitted infections and chlamydia screening in England: 2014. Public Health Engl Infect Rep 9:1–20.
    OpenUrl
  2. 2.↵
    1. Goire N,
    2. Lahra MM,
    3. Chen M,
    4. Donovan B,
    5. Fairley CK,
    6. Guy R,
    7. Kaldor J,
    8. Regan D,
    9. Ward J,
    10. Nissen MD,
    11. Sloots TP,
    12. Whiley DM
    . 2014. Molecular approaches to enhance surveillance of gonococcal antimicrobial resistance. Nat Rev Microbiol 12:223–229. doi:10.1038/nrmicro3217.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Chisholm SA,
    2. Wilson J,
    3. Alexander S,
    4. Tripodo F,
    5. Al-Shahib A,
    6. Schaefer U,
    7. Lythgow K,
    8. Fifer H
    . 24 November 2015. An outbreak of high-level azithromycin resistant Neisseria gonorrhoeae in England. Sex Transm Infect doi:10.1136/sextrans-2015-052312.
    OpenUrlCrossRef
  4. 4.↵
    1. Grad YH,
    2. Goldstein E,
    3. Lipsitch M,
    4. White PJ
    . 2016. Improving control of antibiotic-resistant gonorrhea by integrating research agendas across disciplines: key questions arising from mathematical modeling. J Infect Dis 213:883–890. doi:10.1093/infdis/jiv517.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Bignell C,
    2. FitzGerald M, Guideline Development Group, British Association for Sexual Health and HIV UK
    . 2011. UK national guideline for the management of gonorrhoea in adults, 2011. Int J STD AIDS 22:541–547. doi:10.1258/ijsa.2011.011267.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Le Polain De Waroux O,
    2. Harris RJ,
    3. Hughes G,
    4. Crook PD
    . 2014. The epidemiology of gonorrhoea in London: a Bayesian spatial modelling approach. Epidemiol Infect 142:211–220. doi:10.1017/S0950268813000745.
    OpenUrlCrossRefWeb of Science
  7. 7.↵
    1. O’Rourke M,
    2. Ison CA,
    3. Renton AM,
    4. Spratt BG
    . 1995. Opa-typing: a high-resolution tool for studying the epidemiology of gonorrhoea. Mol Microbiol 17:865–875. doi:10.1111/j.1365-2958.1995.mmi_17050865.x.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Martin IM,
    2. Ison CA,
    3. Aanensen DM,
    4. Fenton KA,
    5. Spratt BG
    . 2004. Rapid sequence-based identification of gonococcal transmission clusters in a large metropolitan area. J Infect Dis 189:1497–1505. doi:10.1086/383047.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    1. Maiden MC,
    2. Bygraves JA,
    3. Feil E,
    4. Morelli G,
    5. Russell JE,
    6. Urwin R,
    7. Zhang Q,
    8. Zhou J,
    9. Zurth K,
    10. Caugant DA,
    11. Feavers IM,
    12. Achtman M,
    13. Spratt BG
    . 1998. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95:3140–3145. doi:10.1073/pnas.95.6.3140.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Didelot X,
    2. Gardy J,
    3. Colijn C
    . 2014. Bayesian inference of infectious disease transmission from whole-genome sequence data. Mol Biol Evol 31:1869–1879. doi:10.1093/molbev/msu121.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Croucher NJ,
    2. Didelot X
    . 2015. The application of genomics to tracing bacterial pathogen transmission. Curr Opin Microbiol 23:62–67. doi:10.1016/j.mib.2014.11.004.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Didelot X,
    2. Walker AS,
    3. Peto TE,
    4. Crook DW,
    5. Wilson DJ
    . 2016. Within-host evolution of bacterial pathogens. Nat Rev Microbiol 14:150–162. doi:10.1038/nrmicro.2015.13.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Harris SR,
    2. Cartwright EJ,
    3. Török ME,
    4. Holden MT,
    5. Brown NM,
    6. Ogilvy-Stuart AL,
    7. Ellington MJ,
    8. Quail MA,
    9. Bentley SD,
    10. Parkhill J,
    11. Peacock SJ,
    12. Török ME,
    13. Holden MTG,
    14. Brown NM,
    15. Ogilvy-Stuart AL,
    16. Ellington MJ,
    17. Quail MA,
    18. Bentley SD,
    19. Parkhill J,
    20. Peacock SJ
    . 2013. Whole-genome sequencing for analysis of an outbreak of meticillin-resistant Staphylococcus aureus: a descriptive study. Lancet Infect Dis 13:130–136. doi:10.1016/S1473-3099(12)70268-2.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    1. Uhlemann A-C,
    2. Dordel J,
    3. Knox JR,
    4. Raven KE,
    5. Parkhill J,
    6. Holden MT,
    7. Peacock SJ,
    8. Lowy FD
    . 2014. Molecular tracing of the emergence, diversification, and transmission of S. aureus sequence type 8 in a New York community. Proc Natl Acad Sci U S A 111:6738–6743. doi:10.1073/pnas.1401006111.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Didelot X,
    2. Eyre DW,
    3. Cule M,
    4. Ip CL,
    5. Ansari MA,
    6. Griffiths D,
    7. Vaughan A,
    8. O’Connor L,
    9. Golubchik T,
    10. Batty EM,
    11. Piazza P,
    12. Wilson DJ,
    13. Bowden R,
    14. Donnelly PJ,
    15. Dingle KE,
    16. Wilcox M,
    17. Walker AS,
    18. Crook DW,
    19. Peto TE,
    20. Harding RM
    . 2012. Microevolutionary analysis of Clostridium difficile genomes to investigate transmission. Genome Biol 13:R118. doi:10.1186/gb-2012-13-12-r118.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Eyre DW,
    2. Cule ML,
    3. Wilson DJ,
    4. Griffiths D,
    5. Vaughan A,
    6. O’Connor L,
    7. Ip CL,
    8. Golubchik T,
    9. Batty EM,
    10. Finney JM,
    11. Wyllie DH,
    12. Didelot X,
    13. Piazza P,
    14. Bowden R,
    15. Dingle KE,
    16. Harding RM,
    17. Crook DW,
    18. Wilcox MH,
    19. Peto TE,
    20. Walker AS
    . 2013. Diverse sources of C. difficile infection identified on whole-genome sequencing. N Engl J Med 369:1195–1205. doi:10.1056/NEJMoa1216064.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Walker TM,
    2. Ip CL,
    3. Harrell RH,
    4. Evans JT,
    5. Kapatai G,
    6. Dedicoat MJ,
    7. Eyre DW,
    8. Wilson DJ,
    9. Hawkey PM,
    10. Crook DW,
    11. Parkhill J,
    12. Harris D,
    13. Walker AS,
    14. Bowden R,
    15. Monk P,
    16. Smith EG,
    17. Peto TE
    . 2013. Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study. Lancet Infect Dis 13:137–146. doi:10.1016/S1473-3099(12)70277-3.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Walker TM,
    2. Lalor MK,
    3. Broda A,
    4. Saldana Ortega L,
    5. Morgan M,
    6. Parker L,
    7. Churchill S,
    8. Bennett K,
    9. Golubchik T,
    10. Giess AP,
    11. Del Ojo Elias C,
    12. Jeffery KJ,
    13. Bowler IC,
    14. Laurenson IF,
    15. Barrett A,
    16. Drobniewski F,
    17. McCarthy ND,
    18. Anderson LF,
    19. Abubakar I,
    20. Thomas HL,
    21. Monk P,
    22. Smith EG,
    23. Walker AS,
    24. Crook DW,
    25. Peto TEA,
    26. Conlon CP
    . 2014. Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007–12, with whole pathogen genome sequences: an observational study. Lancet Respir Med 2:285–292. doi:10.1016/S2213-2600(14)70027-X.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Grad YH,
    2. Kirkcaldy RD,
    3. Trees D,
    4. Dordel J,
    5. Harris SR,
    6. Goldstein E,
    7. Weinstock H,
    8. Parkhill J,
    9. Hanage WP,
    10. Bentley S,
    11. Lipsitch M
    . 2014. Genomic epidemiology of Neisseria gonorrhoeae with reduced susceptibility to cefixime in the USA: a retrospective observational study. Lancet Infect Dis 14:220–226. doi:10.1016/S1473-3099(13)70693-5.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    1. Ward H,
    2. Ison CA,
    3. Day SE,
    4. Martin I,
    5. Ghani AC,
    6. Garnett GP,
    7. Bell G,
    8. Kinghorn G,
    9. Weber JN
    . 2000. A prospective social and molecular investigation of gonococcal transmission. Lancet 356:1812–1817. doi:10.1016/S0140-6736(00)03234-7.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Martin IM,
    2. Ghani A,
    3. Bell G,
    4. Kinghorn G,
    5. Ison CA
    . 2003. Persistence of two genotypes of Neisseria gonorrhoeae during transmission. J Clin Microbiol 41:5609–5614. doi:10.1128/JCM.41.12.5609-5614.2003.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Bilek N,
    2. Martin IM,
    3. Bell G,
    4. Kinghorn GR,
    5. Ison CA,
    6. Spratt BG
    . 2007. Concordance between Neisseria gonorrhoeae genotypes recovered from known sexual contacts. J Clin Microbiol 45:3564–3567. doi:10.1128/JCM.01453-07.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Choudhury B,
    2. Risley CL,
    3. Ghani AC,
    4. Bishop CJ,
    5. Ward H,
    6. Fenton KA,
    7. Ison CA,
    8. Spratt BG
    . 2006. Identification of individuals with gonorrhoea within sexual networks: a population-based study. Lancet 368:139–146. doi:10.1016/S0140-6736(06)69003-X.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Risley CL,
    2. Ward H,
    3. Choudhury B,
    4. Bishop CJ,
    5. Fenton KA,
    6. Spratt BG,
    7. Ison CA,
    8. Ghani AC
    . 2007. Geographical and demographic clustering of gonorrhoea in London. Sex Transm Infect 83:481–487. doi:10.1136/sti.2007.026021.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Croucher NJ,
    2. Page AJ,
    3. Connor TR,
    4. Delaney AJ,
    5. Keane JA,
    6. Bentley SD,
    7. Parkhill J,
    8. Harris SR
    . 2015. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res 43:e15. doi:10.1093/nar/gku1196.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Ezewudo MN,
    2. Joseph SJ,
    3. Castillo-Ramirez S,
    4. Dean D,
    5. Del Rio C,
    6. Didelot X,
    7. Selden RF,
    8. Shafer WM,
    9. Turingan RS,
    10. Unemo M,
    11. Read TD
    . 2015. Population structure of Neisseria gonorrhoeae based on whole genome data and its relationship with antibiotic resistance. PeerJ 3:e806. doi:10.7717/peerj.806.
    OpenUrlCrossRef
  27. 27.↵
    1. Everitt RG,
    2. Didelot X,
    3. Batty EM,
    4. Miller RR,
    5. Knox K,
    6. Young BC,
    7. Bowden R,
    8. Auton A,
    9. Votintseva A,
    10. Larner-Svensson H,
    11. Charlesworth J,
    12. Golubchik T,
    13. Ip CLC,
    14. Godwin H,
    15. Fung R,
    16. Peto TEA,
    17. Walker AS,
    18. Crook DW,
    19. Wilson DJ
    . 2014. Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus. Nat Commun 5:3956.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Chewapreecha C,
    2. Harris SR,
    3. Croucher NJ,
    4. Turner C,
    5. Marttinen P,
    6. Cheng L,
    7. Pessia A,
    8. Aanensen DM,
    9. Mather AE,
    10. Page AJ,
    11. Salter SJ,
    12. Harris D,
    13. Nosten F,
    14. Goldblatt D,
    15. Corander J,
    16. Parkhill J,
    17. Turner P,
    18. Bentley SD
    . 2014. Dense genomic sampling identifies highways of pneumococcal recombination. Nat Genet 46:305–309. doi:10.1038/ng.2895.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Bilek N,
    2. Ison CA,
    3. Spratt BG
    . 2009. Relative contributions of recombination and mutation to the diversification of the Opa gene repertoire of Neisseria gonorrhoeae. J Bacteriol 191:1878–1890. doi:10.1128/JB.01518-08.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Drummond AJ,
    2. Rambaut A
    . 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214. doi:10.1186/1471-2148-7-214.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Lamelas A,
    2. Harris SR,
    3. Röltgen K,
    4. Dangy J-P,
    5. Hauser J,
    6. Kingsley RA,
    7. Connor TR,
    8. Sie A,
    9. Hodgson A,
    10. Dougan G,
    11. Parkhill J,
    12. Bentley SD,
    13. Pluschke G
    . 2014. Emergence of a new epidemic Neisseria meningitidis serogroup A clone in the African meningitis belt: high-resolution picture of genomic changes that mediate immune evasion. mBio 5:e01974–14. doi:10.1128/mBio.01974-14.
    OpenUrlCrossRef
  32. 32.↵
    1. Didelot X,
    2. Bowden R,
    3. Wilson DJ,
    4. Peto TE,
    5. Crook DW
    . 2012. Transforming clinical microbiology with bacterial genome sequencing. Nat Rev Genet 13:601–612. doi:10.1038/nrg3226.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Biek R,
    2. Pybus OG,
    3. Lloyd-Smith JO,
    4. Didelot X
    . 2015. Measurably evolving pathogens in the genomic era. Trends Ecol Evol 30:306–313. doi:10.1016/j.tree.2015.03.009.
    OpenUrlCrossRef
  34. 34.↵
    1. Grenfell BT,
    2. Pybus OG,
    3. Gog JR,
    4. Wood JL,
    5. Daly JM,
    6. Mumford JA,
    7. Holmes EC
    . 2004. Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303:327–332. doi:10.1126/science.1090727.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Pybus OG,
    2. Rambaut A
    . 2009. Evolutionary analysis of the dynamics of viral infectious disease. Nat Rev Genet 10:540–550. doi:10.1038/nrg2583.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Didelot X,
    2. Nell S,
    3. Yang I,
    4. Woltemate S,
    5. van der Merwe S,
    6. Suerbaum S
    . 2013. Genomic evolution and transmission of Helicobacter pylori in two South African families. Proc Natl Acad Sci U S A 110:13880–13885. doi:10.1073/pnas.1304681110.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Garnett GP,
    2. Mertz KJ,
    3. Finelli L,
    4. Levine WC,
    5. St Louis ME
    . 1999. The transmission dynamics of gonorrhoea: modelling the reported behaviour of infected patients from Newark, New Jersey. Philos Trans R Soc Lond B Biol Sci 354:787–797. doi:10.1098/rstb.1999.0431.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Garnett GP,
    2. Bowden FJ
    . 2000. Epidemiology and control of curable sexually transmitted diseases: opportunities and problems. Sex Transm Dis 27:588–599. doi:10.1097/00007435-200011000-00007.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. White PJ,
    2. Ward H,
    3. Cassell JA,
    4. Mercer CH,
    5. Garnett GP
    . 2005. Vicious and virtuous circles in the dynamics of infectious disease and the provision of health care: gonorrhea in Britain as an example. J Infect Dis 192:824–836. doi:10.1086/432004.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Cohen MS,
    2. Cannon JG
    . 1999. Human experimentation with Neisseria gonorrhoeae: progress and goals. J Infect Dis 179(Suppl 2):S375–S379. doi:10.1086/513847.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Mercer CH,
    2. Sutcliffe L,
    3. Johnson AM,
    4. White PJ,
    5. Brook G,
    6. Ross JD,
    7. Dhar J,
    8. Horner P,
    9. Keane F,
    10. Jungmann E,
    11. Sweeney J,
    12. Kinghorn G,
    13. Garnett GG,
    14. Stephenson JM,
    15. Cassell JA
    . 2007. How much do delayed healthcare seeking, delayed care provision, and diversion from primary care contribute to the transmission of STIs? Sex Transm Infect 83:400–405. doi:10.1136/sti.2006.024554.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Satterwhite CL,
    2. Torrone E,
    3. Meites E,
    4. Dunne EF,
    5. Mahajan R,
    6. Ocfemia MC,
    7. Su J,
    8. Xu F,
    9. Weinstock H
    . 2013. Sexually transmitted infections among US women and men: prevalence and incidence estimates, 2008. Sex Transm Dis 40:187–193. doi:10.1097/OLQ.0b013e318286bb53.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Jombart T,
    2. Cori A,
    3. Didelot X,
    4. Cauchemez S,
    5. Fraser C,
    6. Ferguson N
    . 2014. Bayesian reconstruction of disease outbreaks by combining epidemiologic and genomic data. PLoS Comput Biol 10:e1003457. doi:10.1371/journal.pcbi.1003457.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Peters RP,
    2. Nijsten N,
    3. Mutsaers J,
    4. Jansen CL,
    5. Morré SA,
    6. van Leeuwen AP
    . 2011. Evaluation of sexual history-based screening of anatomic sites for Chlamydia trachomatis and Neisseria gonorrhoeae infection in men having sex with men in routine practice. BMC Infect Dis 11:203. doi:10.1186/1471-2334-11-203.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Mercer CH,
    2. Tanton C,
    3. Prah P,
    4. Erens B,
    5. Sonnenberg P,
    6. Clifton S,
    7. Macdowall W,
    8. Lewis R,
    9. Field N,
    10. Datta J,
    11. Copas AJ,
    12. Phelps A,
    13. Wellings K,
    14. Johnson AM
    . 2013. Changes in sexual attitudes and lifestyles in Britain through the life course and over time: findings from the National Surveys of Sexual Attitudes and Lifestyles (Natsal). Lancet 382:1781–1794. doi:10.1016/S0140-6736(13)62035-8.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    1. Mayer KH,
    2. Krakower DS
    . 2015. Uncoupling epidemiological synergy: new opportunities for HIV prevention for men who have sex with men. Clin Infect Dis 61:288–290. doi:10.1093/cid/civ291.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Wasserheit JN
    . 1992. Epidemiological synergy: interrelationships between human immunodeficiency virus infection and other sexually transmitted diseases. Sex Transm Dis 19:61–77. doi:10.1097/00007435-199219020-00001.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    1. Golden MR,
    2. Stekler J,
    3. Hughes JP,
    4. Wood RW
    . 2008. HIV serosorting in men who have sex with men: is it safe? J Acquir Immune Defic Syndr 49:212–218. doi:10.1097/QAI.0b013e31818455e8.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Malek R,
    2. Mitchell H,
    3. Furegato M,
    4. Simms I,
    5. Mohammed H,
    6. Nardone A,
    7. Hughes G
    . 2015. Contribution of transmission in HIV-positive men who have sex with men to evolving epidemics of sexually transmitted infections in England: an analysis using multiple data sources, 2009–2013. Euro Surveill 20:7–15.
    OpenUrlPubMed
  50. 50.↵
    1. Donnelly P,
    2. Tavaré S
    . 1995. Coalescents and genealogical structure under neutrality. Annu Rev Genet 29:401–421. doi:10.1146/annurev.ge.29.120195.002153.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Wilson DJ,
    2. Falush D,
    3. McVean G
    . 2005. Germs, genomes and genealogies. Trends Ecol Evol 20:39–45. doi:10.1016/j.tree.2004.10.009.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Koelle K,
    2. Rasmussen DA
    . 2012. Rates of coalescence for common epidemiological models at equilibrium. J R Soc Interface 9:997–1007. doi:10.1098/rsif.2011.0495.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Dearlove B,
    2. Wilson DJ
    . 2013. Coalescent inference for infectious disease: meta-analysis of hepatitis C. Philos Trans R Soc Lond B Biol Sci 368:20120314. doi:10.1098/rstb.2012.0314.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Fingerhuth SM,
    2. Bonhoeffer S,
    3. Low N,
    4. Althaus CL
    . 2015. Antibiotic-resistant Neisseria gonorrhoeae spread faster with more treatment, not more sexual partners. PeerJ PrePrints 3:e1805. doi:10.1371/journal.ppat.1005611.
    OpenUrlCrossRef
  55. 55.↵
    1. Snyder LA,
    2. Cole JA,
    3. Pallen MJ
    . 2009. Comparative analysis of two Neisseria gonorrhoeae genome sequences reveals evidence of mobilization of Correia repeat enclosed elements and their role in regulation. BMC Genomics 10:70. doi:10.1186/1471-2164-10-70.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Aziz RK,
    2. Bartels D,
    3. Best AA,
    4. DeJongh M,
    5. Disz T,
    6. Edwards RA,
    7. Formsma K,
    8. Gerdes S,
    9. Glass EM,
    10. Kubal M,
    11. Meyer F,
    12. Olsen GJ,
    13. Olson R,
    14. Osterman AL,
    15. Overbeek RA,
    16. McNeil LK,
    17. Paarmann D,
    18. Paczian T,
    19. Parrello B,
    20. Pusch GD,
    21. Reich C,
    22. Stevens R,
    23. Vassieva O,
    24. Vonstein V,
    25. Wilke A,
    26. Zagnitko O
    . 2008. The RAST server: rapid annotations using subsystems technology. BMC Genomics 9:75. doi:10.1186/1471-2164-9-75.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Zhou Y,
    2. Liang Y,
    3. Lynch KH,
    4. Dennis JJ,
    5. Wishart DS
    . 2011. PHAST: a fast phage search tool. Nucleic Acids Res 39:W347–W352. doi:10.1093/nar/gkr485.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Rice P,
    2. Longden I,
    3. Bleasby A
    . 2000. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet 16:276–277. doi:10.1016/S0168-9525(00)02024-2.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. McKenna A,
    2. Hanna M,
    3. Banks E,
    4. Sivachenko A,
    5. Cibulskis K,
    6. Kernytsky A,
    7. Garimella K,
    8. Altshuler D,
    9. Gabriel S,
    10. Daly M,
    11. DePristo MA
    . 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303. doi:10.1101/gr.107524.110.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Harris SR,
    2. Feil EJ,
    3. Holden MT,
    4. Quail MA,
    5. Nickerson EK,
    6. Chantratita N,
    7. Gardete S,
    8. Tavares A,
    9. Day N,
    10. Lindsay JA,
    11. Edgeworth JD,
    12. de Lencastre H,
    13. Parkhill J,
    14. Peacock SJ,
    15. Bentley SD
    . 2010. Evolution of MRSA during hospital transmission and intercontinental spread. Science 327:469–474. doi:10.1126/science.1182395.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Guindon S,
    2. Dufayard J-F,
    3. Lefort V,
    4. Anisimova M,
    5. Hordijk W,
    6. Gascuel O
    . 2010. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59:307–321. doi:10.1093/sysbio/syq010.
    OpenUrlCrossRefPubMedWeb of Science
  62. 62.↵
    1. Cori A,
    2. Ferguson NM,
    3. Fraser C,
    4. Cauchemez S
    . 2013. A new framework and software to estimate time-varying reproduction numbers during epidemics. Am J Epidemiol 178:1505–1512. doi:10.1093/aje/kwt133.
    OpenUrlCrossRefPubMed
Back to top
Download PDF
Citation Tools
Print

Alerts
Email

Related Articles

Cited By...