Emergence of SARS-CoV-2 Delta Variant and Effect of Nonpharmaceutical Interventions, British Columbia, Canada.

In British Columbia, Canada, initial growth of the SARS-CoV-2 Delta variant was slower than that reported in other jurisdictions. Delta became the dominant variant (>50% prevalence) within ≈7-13 weeks of first detection in regions within the United Kingdom and United States. In British Columbia, it remained at <10% of weekly incident COVID-19 cases for 13 weeks after first detection on March 21, 2021, eventually reaching dominance after 17 weeks. We describe the growth of Delta variant cases in British Columbia during March 1-June 30, 2021, and apply retrospective counterfactual modeling to examine factors for the initially low COVID-19 case rate after Delta introduction, such as vaccination coverage and nonpharmaceutical interventions. Growth of COVID-19 cases in the first 3 months after Delta emergence was likely limited in British Columbia because additional nonpharmaceutical interventions were implemented to reduce levels of contact at the end of March 2021, soon after variant emergence.

Comparative Single-Dose mRNA and ChAdOx1 Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2, Including Variants of Concern: Test-Negative Design, British Columbia, Canada.

BACKGROUND: In British Columbia, Canada, most adults 50-69 years old became eligible for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine in April 2021, with chimpanzee adenoviral vectored vaccine (ChAdOx1) restricted to ≥55-year-olds and second doses deferred ≥6 weeks to optimize single-dose coverage. METHODS: Among adults 50-69 years old, single-dose messenger RNA (mRNA) and ChAdOx1 vaccine effectiveness (VE) against SARS-CoV-2 infection and hospitalization, including variant-specific, was assessed by test-negative design between 4 April and 2 October 2021. RESULTS: Single-dose VE included 11 861 cases and 99 544 controls. Median of postvaccination follow-up was 32 days (interquartile range, 15-52 days). Alpha, Gamma, and Delta variants comprised 23%, 18%, and 56%, respectively, of genetically characterized viruses. At 21-55 days postvaccination, single-dose mRNA and ChAdOx1 VE (95% confidence interval [CI]) was 74% (71%-76%) and 59% (53%-65%) against any infection and 86% (80%-90%) and 94% (85%-97%) against hospitalization, respectively. VE (95% CI) was similar against Alpha and Gamma infections for mRNA (80% [76%-84%] and 80% [75%-84%], respectively) and ChAdOx1 (69% [60%-76%] and 66% [56%-73%], respectively). mRNA VE was lower at 63% (95% CI, 56%-69%) against Delta but 85% (95% CI, 71%-92%) against Delta-associated hospitalization (nonestimable for ChAdOx1). CONCLUSIONS: A single mRNA or ChAdOx1 vaccine dose gave important protection against SARS-CoV-2, including early variants of concern. ChAdOx1 VE was lower against infection, but 1 dose of either vaccine reduced the hospitalization risk by >85% to at least 8 weeks postvaccination. Findings inform program options, including longer dosing intervals.

Single-dose mRNA Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), Including Alpha and Gamma Variants: A Test-negative Design in Adults 70 Years and Older in British Columbia, Canada.

BACKGROUND: Randomized-controlled trials of messenger RNA (mRNA) vaccine protection against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) included relatively few elderly participants. We assess single-dose mRNA vaccine effectiveness (VE) in adults ≥ 70 years old in British Columbia, Canada, where second doses were deferred by up to 16 weeks and where a spring 2021 wave uniquely included codominant circulation of Alpha (B.1.1.7) and Gamma (P.1) variants of concern (VOC). METHODS: Analyses included community-dwelling adults ≥ 70 years old with specimen collection between 4 April (epidemiological week 14) and 1 May (week 17) 2021. Adjusted VE was estimated by test-negative design. Cases were reverse-transcription polymerase chain reaction (RT-PCR) test-positive for SARS-CoV-2, and controls were test-negative. Vaccine status was defined by receipt of a single-dose ≥ 21 days before specimen collection, but a range of intervals was assessed. Variant-specific VE was estimated against viruses genetically characterized as Alpha, Gamma or non-VOC lineages. RESULTS: VE analyses included 16 993 specimens: 1226 (7%) test-positive cases and 15 767 test-negative controls. Of 1131 (92%) genetically characterized viruses, 509 (45%), 314 (28%), and 276 (24%) were Alpha, Gamma, and non-VOC lineages, respectively. At 0-13 days postvaccination, VE was negligible at 14% (95% confidence interval [CI], 0-26) but increased from 43% (95% CI, 30-53) at 14-20 days to 75% (95% CI, 63-83) at 35-41 days postvaccination. VE at ≥ 21 days postvaccination was 65% (95% CI, 58-71) overall: 72% (95% CI, 58-81), 67% (95% CI, 57-75), and 61% (95% CI, 45-72) for non-VOC, Alpha, and Gamma variants, respectively. CONCLUSIONS: A single dose of mRNA vaccine reduced the risk of SARS-CoV-2 by about two-thirds in adults ≥ 70 years old, with protection only minimally reduced against Alpha and Gamma variants.

Two-Dose Severe Acute Respiratory Syndrome Coronavirus 2 Vaccine Effectiveness With Mixed Schedules and Extended Dosing Intervals: Test-Negative Design Studies From British Columbia and Quebec, Canada.

BACKGROUND: The Canadian coronavirus disease 2019 (COVID-19) immunization strategy deferred second doses and allowed mixed schedules. We compared 2-dose vaccine effectiveness (VE) by vaccine type (mRNA and/or ChAdOx1), interval between doses, and time since second dose in 2 of Canada's larger provinces. METHODS: Two-dose VE against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or hospitalization among adults ≥18 years, including due to Alpha, Gamma, and Delta variants of concern (VOCs), was assessed ≥14 days postvaccination by test-negative design studies separately conducted in British Columbia and Quebec, Canada, between 30 May and 27 November (epi-weeks 22-47) 2021. RESULTS: In both provinces, all homologous or heterologous mRNA and/or ChAdOx1 2-dose schedules were associated with ≥90% reduction in SARS-CoV-2 hospitalization risk for ≥7 months. With slight decline from a peak of >90%, VE against infection was ≥80% for ≥6 months following homologous mRNA vaccination, lower by ∼10% when both doses were ChAdOx1 but comparably high following heterologous ChAdOx1 + mRNA receipt. Findings were similar by age group, sex, and VOC. VE was significantly higher with longer 7-8-week versus manufacturer-specified 3-4-week intervals between mRNA doses. CONCLUSIONS: Two doses of any mRNA and/or ChAdOx1 combination gave substantial and sustained protection against SARS-CoV-2 hospitalization, spanning Delta-dominant circulation. ChAdOx1 VE against infection was improved by heterologous mRNA series completion. A 7-8-week interval between first and second doses improved mRNA VE and may be the optimal schedule outside periods of intense epidemic surge. Findings support interchangeability and extended intervals between SARS-CoV-2 vaccine doses, with potential global implications for low-coverage areas and, going forward, for children.

Cov2clusters: genomic clustering of SARS-CoV-2 sequences.

BACKGROUND: The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada. RESULTS: We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data. CONCLUSIONS: Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

AMPlify: attentive deep learning model for discovery of novel antimicrobial peptides effective against WHO priority pathogens.

BACKGROUND: Antibiotic resistance is a growing global health concern prompting researchers to seek alternatives to conventional antibiotics. Antimicrobial peptides (AMPs) are attracting attention again as therapeutic agents with promising utility in this domain, and using in silico methods to discover novel AMPs is a strategy that is gaining interest. Such methods can sift through large volumes of candidate sequences and reduce lab screening costs. RESULTS: Here we introduce AMPlify, an attentive deep learning model for AMP prediction, and demonstrate its utility in prioritizing peptide sequences derived from the Rana [Lithobates] catesbeiana (bullfrog) genome. We tested the bioactivity of our predicted peptides against a panel of bacterial species, including representatives from the World Health Organization's priority pathogens list. Four of our novel AMPs were active against multiple species of bacteria, including a multi-drug resistant isolate of carbapenemase-producing Escherichia coli. CONCLUSIONS: We demonstrate the utility of deep learning based tools like AMPlify in our fight against antibiotic resistance. We expect such tools to play a significant role in discovering novel candidates of peptide-based alternatives to classical antibiotics.

Rapid Increase in SARS-CoV-2 P.1 Lineage Leading to Codominance with B.1.1.7 Lineage, British Columbia, Canada, January-April 2021.

Several severe acute respiratory syndrome coronavirus 2 variants of concern (VOCs) emerged in late 2020; lineage B.1.1.7 initially dominated globally. However, lineages B.1.351 and P.1 represent potentially greater risk for transmission and immune escape. In British Columbia, Canada, B.1.1.7 and B.1.351 were first identified in December 2020 and P.1 in February 2021. We combined quantitative PCR and whole-genome sequencing to assess relative contribution of VOCs in nearly 67,000 infections during the first 16 weeks of 2021 in British Columbia. B.1.1.7 accounted for <10% of screened or sequenced specimens early on, increasing to >50% by week 8. P.1 accounted for <10% until week 10, increased rapidly to peak at week 12, and by week 13 codominated within 10% of rates of B.1.1.7. B.1.351 was a minority throughout. This rapid expansion of P.1 but suppression of B.1.351 expands our understanding of population-level VOC patterns and might provide clues to fitness determinants for emerging VOCs.

Self-Collected Saline Gargle Samples as an Alternative to Health Care Worker-Collected Nasopharyngeal Swabs for COVID-19 Diagnosis in Outpatients.

We assessed the performance, stability, and user acceptability of swab-independent self-collected saliva and saline mouth rinse/gargle sample types for the molecular detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in adults and school-aged children. Outpatients who had recently been diagnosed with COVID-19 or were presenting with suspected COVID-19 were asked to have a nasopharyngeal (NP) swab collected and provide at least one self-collected sample type. Participants were also asked about sample acceptability using a five-point Likert scale. For those previously diagnosed with COVID-19, all samples underwent real-time PCR testing using a lab-developed assay, and the majority were also tested using an FDA-authorized assay. For those presenting with suspected COVID-19, only those with a positive nasopharyngeal swab sample went on to have other samples tested. Saline mouth rinse/gargle and saliva samples were tested daily at time zero, day 1, and day 2 to assess nucleic acid stability at room temperature. Fifty participants (aged 4 to 71 years) were included; of these, 40 had at least one positive sample and were included in the primary sample yield analysis. Saline mouth rinse/gargle samples had a sensitivity of 98% (39/40), while saliva samples had a sensitivity of 79% (26/33). Both saline mouth rinse/gargle and saliva samples showed stable viral RNA detection after 2 days of room temperature storage. Mouth rinse/gargle samples had the highest (mean, 4.9) and health care worker (HCW)-collected NP swabs had the lowest acceptability scores (mean, 3.1). In conclusion, saline mouth rinse/gargle samples demonstrated higher combined user acceptability ratings and analytical performance than saliva and HCW-collected NP swabs. This sample type is a promising swab-independent option, particularly for outpatient self-collection in adults and school-aged children.

Outcomes of Critically Ill Patients Who Have Serotype 5 Invasive Pneumococcal Disease.

PURPOSE: To determine whether invasive pneumococcal disease (IPD) due to serotype 5, which occurred as a local outbreak in 2006 to 2007, is associated with intensive care unit (ICU) admission, hospital mortality, or organ supports in those who are critically ill. MATERIALS AND METHODS: Retrospective review of patients who presented with IPD to 2 tertiary hospitals in Vancouver, Canada, from July 2004 to June 2007. We compared patient characteristics, interventions, and outcomes between patients who had serotype 5 and other serotypes using bivariate and multivariate analyses. RESULTS: A total of 149 patients had serotype 5 and 106 had nonserotype 5. Patients with serotype 5 were younger, had lower prevalence of comorbid diseases, and had higher rates of substance use than patients with nonserotype 5. There were no differences in chest tube placement for complications of pneumonia or in ICU admission. Frequency of necrotizing pneumonia and hospital mortality were lower in the serotype 5 group. For the 71 patients with IPD who were admitted to ICU, there was no difference in severity of illness, ICU length of stay, or ICU mortality between the groups. There was also no difference in organ supports except that the serotype 5 group was more likely to receive vasopressors. CONCLUSION: Serotype 5 in patients who have IPD is associated with no difference in ICU admission but with increased use of vasopressors and lower hospital mortality.

Mycoplasma genitalium Antibiotic Resistance-Mediating Mutations in Canadian Women With or Without Chlamydia Trachomatis Infection.

Testing remnant Aptima specimens from women infected with Chlamydia trachomatis detected 13.4% (53/396) with Mycoplasma genitalium compared with 5.4% (22/406) in matched C. trachomatis-negative women. Overall, 9.4% (provincial ranges of 3-20%) were infected with M. genitalium and resistance mediating mutations were found in 47.3% (26/55) to macrolides and 1.9% (1/53) to fluoroquinolones by sequencing.

Characterization of isolates of Eisenbergiella tayi, a strictly anaerobic Gram-stain variable bacillus recovered from human clinical materials in Canada.

Eisenbergiella gen. nov. was proposed in 2014 to describe an obligate anaerobic, structurally Gram-positive but Gram-stain-negative-appearing bacillus recovered from the blood culture of an elderly Israeli man. Here, we describe features for eight blood culture isolates as well one appendix-derived isolate, recovered from seven patients located in two Canadian provinces, which by 16S rRNA gene sequencing, were identifiable as Eisenbergiella tayi, the sole validly- named species in this genus. After whole genome sequencing, isolates were found to be essentially identical (96.8-98.7% identity) to each other and to E. tayi DSM 26961T, after comparison using the ANIb tool and in silico DNA-DNA hydridization. All isolates were observed to have remarkably large genomes (7.1-8.3 Mb) with a G + C content of 46.5%-46.9%.

Epidemiologic and Genotypic Review of Carbapenemase-Producing Organisms in British Columbia, Canada, between 2008 and 2014.

Carbapenemase-producing organisms (CPOs) are a serious emerging problem for health care facilities worldwide. Owing to their resistance to most antimicrobial therapies, CPOs are difficult to treat and pose a challenge for infection prevention and control. Since 2010, lab-based surveillance for CPOs and PCR-based testing were implemented in British Columbia (BC), Canada. A review of CPOs in BC from 2008 to March 2014 was done to characterize the resistance mechanisms and possible clonal strain transmission and to compare pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and plasmid restriction fragment length polymorphism (RFLP) as molecular typing tools. During this study period, a total of 177 CPO cases were identified. Patient demographics and travel history were reviewed, and a descriptive analysis was carried out. PFGE profiles, MLST, and plasmid RFLP analysis for a subset of Escherichia coli, Klebsiella pneumoniae, and Enterobacter species isolates were obtained and analyzed. Our findings demonstrate that CPOs have been increasing in number in BC over time, from 1 isolate/year retrospectively identified in 2008 and 2009 to 82 isolates in 2013 and 30 isolates in the first quarter of 2014. Overall, K. pneumoniae isolates lack clonality, although some seemingly related clusters have been found. Plasmid analysis showed evidence of the spread of plasmids carrying carbapenemase-encoding genes between the examined isolates. Analysis of Enterobacter cloacae isolates revealed a more clonal nature of these CPOs in BC. The presence of related clusters provides evidence of interpatient organism transmission both within and between institutions. Although in our study, NDM-harboring E. cloacae isolates appeared to spread clonally, the spread of carbapenem resistance in K. pneumoniae seems to be plasmid mediated.

Review of 16S and ITS Direct Sequencing Results for Clinical Specimens Submitted to a Reference Laboratory.

We evaluated the performance of 16S and internal transcribed spacer (ITS) region amplification and sequencing of rDNA from clinical specimens, for the respective detection and identification of bacterial and fungal pathogens. Direct rDNA amplification of 16S and ITS targets from clinical samples was performed over a 4-year period and reviewed. All specimens were from sterile sites and submitted to a reference laboratory for evaluation. Results of 16S and ITS were compared to histopathology, Gram and/or calcofluor stain microscopy results. A total of 277 16S tests were performed, with 64 (23%) positive for the presence of bacterial DNA. Identification of an organism was more likely in microscopy positive 16S samples 14/21 (67%), compared to 35/175 (20%) of microscopy negative samples. A total of 110 ITS tests were performed, with 14 (13%) positive. The yield of microscopy positive ITS samples, 9/44 (21%), was higher than microscopy negative samples 3/50 (6%). Given these findings, 16S and ITS are valuable options for culture negative specimens from sterile sites, particularly in the setting of positive microscopy findings. Where microscopy results are negative, the limited sensitivity of 16S and ITS in detecting and identifying an infectious agent needs to be considered.

Characterization of Clostridium difficile Strains in British Columbia, Canada: A Shift from NAP1 Majority (2008) to Novel Strain Types (2013) in One Region.

Background. Clostridium difficile is a major cause of gastrointestinal illness. Epidemic NAP1 strains contain toxins A and B, a deletion in repressor tcdC, and a binary toxin. Objectives. To determine the molecular epidemiology of C. difficile in British Columbia and compare between two time points in one region. Methods. C. difficile isolates from hospital and community laboratories (2008) and one Island Health hospital laboratory (2013) were characterized by pulsed-field gel electrophoresis, PCR-ribotyping, toxin possession, tcdC genotype, and antimicrobial susceptibility. Results. In 2008, 42.7% of isolates had NAP1 designation. Hospital-collected isolates were associated with older patients and more NAP1 types. Unlike other isolates, most NAP1 isolates possessed binary toxin and a 19 bp loss in tcdC. All isolates were susceptible to metronidazole and vancomycin. A 2013 follow-up revealed a 28.9% decrease in NAP1 isolates and 20.0% increase in isolates without NAP designation in one region. Then, community-associated cases were seen in younger patients, while NAP types were evenly distributed. Isolates without NAP designation did not cluster with a PFGE pattern or ribotype. Conclusions. Evaluation of C. difficile infections within British Columbia revealed demographic associations, epidemiological shifts, and characteristics of strain types. Continuous surveillance of C. difficile will enable detection of emerging strains.

Rapid Detection of KPC, NDM, and OXA-48-Like Carbapenemases by Real-Time PCR from Rectal Swab Surveillance Samples.

We describe a multiplex real-time PCR assay for use on the ABI 7500 Fast TaqMan platform to detect all currently described Klebsiella pneumoniae carbapenemases (KPC), New Delhi metallo-ß-lactamases (NDM), and the OXA-48-like family of carbapenemases from bacterial culture lysates or sample enrichment broth lysates.

De novo synthetic antimicrobial peptide design with a recurrent neural network.

Antibiotic resistance is recognized as an imminent and growing global health threat. New antimicrobial drugs are urgently needed due to the decreasing effectiveness of conventional small-molecule antibiotics. Antimicrobial peptides (AMPs), a class of host defense peptides, are emerging as promising candidates to address this need. The potential sequence space of amino acids is combinatorially vast, making it possible to extend the current arsenal of antimicrobial agents with a practically infinite number of new peptide-based candidates. However, mining naturally occurring AMPs, whether directly by wet lab screening methods or aided by bioinformatics prediction tools, has its theoretical limit regarding the number of samples or genomic/transcriptomic resources researchers have access to. Further, manually designing novel synthetic AMPs requires prior field knowledge, restricting its throughput. In silico sequence generation methods are gaining interest as a high-throughput solution to the problem. Here, we introduce AMPd-Up, a recurrent neural network based tool for de novo AMP design, and demonstrate its utility over existing methods. Validation of candidates designed by AMPd-Up through antimicrobial susceptibility testing revealed that 40 of the 58 generated sequences possessed antimicrobial activity against Escherichia coli and/or Staphylococcus aureus. These results illustrate that AMPd-Up can be used to design novel synthetic AMPs with potent activities.

Cephalosporin and azithromycin susceptibility in Neisseria gonorrhoeae isolates by site of infection, British Columbia, 2006 to 2011.

BACKGROUND: Widespread resistance of Neisseria gonorrhoeae to penicillin, tetracycline, and fluoroquinolones has challenged effective treatment and control; recent international case reports of cefixime, ceftriaxone, and azithromycin resistance suggest that the remaining treatment options are now additionally threatened. To explore trends in antimicrobial susceptibility of N. gonorrhoeae, we reviewed provincial laboratory data from British Columbia, 2006 to 2011. METHODS: Susceptibility testing was performed for all N. gonorrhoeae isolates detected in-house or forwarded to the reference laboratory. Resistance or intermediate resistance (nonsusceptibility) was defined by standard breakpoints for penicillin, tetracycline, ciprofloxacin, and spectinomycin. Elevated minimum inhibitory concentrations (MICs) at serial dilutions of 0.064 µg/mL or greater were explored for cefixime/ceftriaxone and 0.5 µg/mL or greater for azithromycin. Nonsusceptibility/elevated MIC was compared by year, site of infection, sex, and age. RESULTS: A total of 1837 isolates representing 22% of all reported gonorrhea cases were analyzed. Nonsusceptibility to penicillin was established at baseline. Nonsusceptibility to tetracycline and ciprofloxacin increased over the study period, reaching 96% and 36%, respectively, in 2011. Sixteen isolates (1%) had a cefixime MIC of 0.25 µg/mL (none ≥0.5), none had a ceftriaxone MIC of 0.25 µg/mL or greater, and 15 (1%) had an azithromycin MIC of 2.0 µg/mL or greater. Elevated MIC of these agents showed an increasing trend over time. Nonsusceptibility and elevated MIC were consistently highest at the rectal and pharyngeal sites and higher in isolates from males, including when stratified to the pharyngeal site. INTERPRETATION: Increases in elevated MIC of cefixime/ceftriaxone/azithromycin were superimposed on a background of established resistance to penicillin, tetracycline, and ciprofloxacin and may signal impending gonococcal resistance to first-line treatments. Ongoing surveillance will inform timely shifts in treatment recommendations.

The utility of SARS-CoV-2 genomic data for informative clustering under different epidemiological scenarios and sampling.

OBJECTIVES: Clustering pathogen sequence data is a common practice in epidemiology to gain insights into the genetic diversity and evolutionary relationships among pathogens. We can find groups of cases with a shared transmission history and common origin, as well as identifying transmission hotspots. Motivated by the experience of clustering SARS-CoV-2 cases using whole genome sequence data during the COVID-19 pandemic to aid with public health investigation, we investigated how differences in epidemiology and sampling can influence the composition of clusters that are identified. METHODS: We performed genomic clustering on simulated SARS-CoV-2 outbreaks produced with different transmission rates and levels of genomic diversity, along with varying the proportion of cases sampled. RESULTS: In single outbreaks with a low transmission rate, decreasing the sampling fraction resulted in multiple, separate clusters being identified where intermediate cases in transmission chains are missed. Outbreaks simulated with a high transmission rate were more robust to changes in the sampling fraction and largely resulted in a single cluster that included all sampled outbreak cases. When considering multiple outbreaks in a sampled jurisdiction seeded by different introductions, low genomic diversity between introduced cases caused outbreaks to be merged into large clusters. If the transmission and sampling fraction, and diversity between introductions was low, a combination of the spurious break-up of outbreaks and the linking of closely related cases in different outbreaks resulted in clusters that may appear informative, but these did not reflect the true underlying population structure. Conversely, genomic clusters matched the true population structure when there was relatively high diversity between introductions and a high transmission rate. CONCLUSION: Differences in epidemiology and sampling can impact our ability to identify genomic clusters that describe the underlying population structure. These findings can help to guide recommendations for the use of pathogen clustering in public health investigations.