ABSTRACT
Objectives Mathematical modelling and genomic analyses are powerful methods for investigating the transmission dynamics of Neisseria gonorrhoeae, however, often make the implicit assumption that N. gonorrhoeae isolates at different anatomical sites within the same individual are the same strain.
Methods In this study, two approaches were used to explore genetic diversity. First, we examined a collection of stored, clinical N. gonorrhoeae isolates sourced from multiple anatomical sites of single individuals attending a sexual health clinic in Melbourne from 2011-2019. Second, we obtained multiple colony picks from primary clinical samples from individuals attending a sexual health clinic in Melbourne from 2019-2020. Whole genome sequencing and a variety of bioinformatics approaches were used to determine both within-host and within-sample genetic diversity.
Results Thirty-seven individuals were identified that had cultured N. gonorrhoeae from two or more anatomical sites (urogenital, anorectal, or oropharyngeal), with a final dataset of 105 isolates. In 35/37 (94.6%) individuals, infections were highly similar at the genetic level, with identical MLST and NG-MAST profiles. Pairwise comparisons of isolates within each individual indicated that the maximum within-host pairwise SNP distance was 13 SNPs (median = 1, IQR: 0-3). Notably, four distinct multi-individual phylogenetic clusters were identified, where the maximum pairwise SNP distance was 19 SNPs (median = 6, IQR = 2-11). Similarly, comparisons of isolates within each primary sample indicated that the maximum pairwise SNP distance was 8 SNPs (median = 2, IQR:1-3).
Conclusions This study suggests that in most cases of multi-site infection, the same strain of N. gonorrhoeae causes the infection at each anatomical site. However, WGS data alone cannot differentiate between the same infecting strain or (re)infections from the same transmission network. These data guide recommendations regarding optimal bioinformatic approaches to infer genetic relatedness of N. gonorrhoeae and will help inform future studies of gonorrhoea transmission and epidemiology.
KEY MESSAGES
In most cases of multi-site infection, the same strain of N. gonorrhoeae causes infection at each anatomical site (urogenital, anorectal, and oropharyngeal) within an individual.
It is not possible, using genomics alone, to determine if an infection is from a single individual or multiple individuals within the same sexual network.
Current literature reports an array of different bioinformatic methods for describing genetic diversity, highlighting that SNP thresholds are not directly transferable between studies.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
D.A.W, E.P.F.C and C.K.F are supported by an Investigator Grant from the National Health and Medical Research Council (NHMRC) of Australia (GNT1174555, GNT1172873 and GNT1172900, respectively).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The study was approved by the Alfred Hospital Ethics Committee (399/19).
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
Sequence data (reads) have been submitted to the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov) under the BioProject Accession: PRJNA520805. All supporting data have been provided within the article or supplementary data files.
