Use of whole-genome sequencing to identify clusters of Shigella flexneri associated with sexual transmission in men who have sex with men in England: a validation study using linked behavioural data

Since the 1970s, shigellosis has been reported as a sexually transmissible infection, and in recent years, genomic data have revealed the breadth of Shigella spp. transmission among global networks of men who have sex with men (MSM). In 2015, Public Health England (PHE) introduced routine whole-genome sequencing (WGS) of Shigella spp. to identify transmission clusters. However, limited behavioural information for the cases hampers interpretation. We investigated whether WGS can distinguish between clusters representing sexual transmission in MSM and clusters representing community (non-sexual) transmission to inform infection control. WGS data for Shigella flexneri from August 2015 to July 2017 were aggregated into single linkage clusters based on SNP typing using a range of SNP distances (the standard for Shigella surveillance at PHE is 10 SNPs). Clusters were classified as 'adult male', 'household', 'travel-associated' or 'community' using routine demographic data submitted alongside laboratory cultures. From August 2015 to March 2017, PHE contacted those with shigellosis as part of routine public-health follow-up and collected exposure data on a structured questionnaire, which for the first time included questions about sexual identity and behaviour. The questionnaire data were used to determine whether clusters classified as 'adult male' represented likely sexual transmission between men, thereby validating the use of the SNP clustering tool for informing appropriate public-health responses. Overall, 1006 S. flexneri cases were reported, of which 563 clustered with at least one other case (10-SNP threshold). Linked questionnaire data were available for 106 clustered cases, of which 84.0 % belonged to an 'adult male' cluster. At the 10-SNP threshold, 95.1 % [95 % confidence interval (CI) 88.0-98.1%] of MSM belonged to an 'adult male' cluster, while 73.2 % (95 % CI 49.1-87.5%) of non-MSM belonged to a 'community' or 'travel-associated' cluster. At the 25-SNP threshold, all MSM (95 % CI 96.0-100%) belonged to an 'adult male' cluster and 77.8 % (95 % CI 59.2-89.4%) of non-MSM belonged to a 'community' or 'travel-associated' cluster. Within one phylogenetic clade of S. flexneri, 9 clusters were identified (7 'adult male'; 2 'community') using a 10-SNP threshold, while a single 'adult male' cluster was identified using a 25-SNP threshold. Genotypic markers of azithromycin resistance were detected in 84.5 % (294/348) of 'adult male' cases and 20.9 % (9/43) of cases in other clusters (10-SNP threshold), the latter of which contained gay-identifying men who reported recent same-sex sexual contact. Our study suggests that SNP clustering can be used to identify Shigella clusters representing likely sexual transmission in MSM to inform infection control. Defining clusters requires a flexible approach in terms of genetic relatedness to ensure a clear understanding of underlying transmission networks