Introduction
Invasive meningococcal disease (IMD) is caused by Neisseria meningitidis, the meningococcus. Despite improved treatment and the development of effective vaccines, mortality remains high even in high-income countries. Many surviving patients suffer sequelae such as skin scarring, limb amputation, hearing loss or learning disabilities [Reference Brandtzaeg, Van Deuren and Walker1]. Meningococci are classified into 12 serogroups based on the structure of the polysaccharide capsule. IMD is usually caused by serogroups A, B, C, W, X or Y, which all have epidemic potential. Uncapsulated meningococci seldom cause invasive disease, but are commonly detected in the oropharynx in healthy individuals [Reference Caugant and Brynildsrud2].
Meningococcal carriage is a prerequisite for developing IMD and essential for transmission. N. meningitidis is only found in humans and is transmitted through respiratory secretions resulting in transient nasopharyngeal colonisation. Carriage is usually asymptomatic and peaks in adolescents and young adults [Reference Christensen3]. Smoking, kissing, overcrowding and frequenting bars are known risk factors for meningococcal carriage [Reference Maclennan4]. Only a small fraction of carriers will develop invasive disease, usually shortly after acquisition of the bacterium.
The incidence of IMD in Norway has been below 0.5 per 100 000 in the general population in the last 3 years [5]. After a peak in 2010 with 5.3 cases per 100 000, the mean incidence in 15–19-year-olds was 1.2 per 100 000 in the past 3 years. Most cases in Norwegian teenagers have been associated with the month-long ‘russ celebration’, a tradition with heavy drinking and partying among graduates from upper secondary school [Reference Taylor6]. Serogroup Y dominated (75%) in 15–19-year-olds, even though the two fatal cases in 2017–2019 were caused by serogroups C and W, clonal complex (cc) 11 [5]. In the same period, the mean incidence in other age groups was 0.8 per 100 000 in children <5 years and 0.3 per 100 000 in both 5–14-year-olds and adults >19 years. Serogroup B dominated in children <5 years, while cases in the 5–14-year-olds and adults >19 years were caused by serogroups B, C, W and Y.
Meningococcal vaccination is not part of the national immunisation programme in Norway. A substantial increase in the incidence of IMD was observed in teenagers involved in the russ celebration in 2009–2010. Therefore, the Norwegian Institute of Public Health (NIPH) started recommending meningococcal ACWY conjugate vaccine (MCV4) to 17–19-year-olds engaged in activities that increase the risk of IMD, such as smoking, sharing drinking bottles and participating in youth gatherings, in 2011. The recommendation was extended to 16–19-year-olds from 2012. The school health services usually administer vaccination at the student's own cost. Vaccination is usually offered during the last year of upper secondary school. Uptake of vaccination among graduating students aged 18–19 years has increased from 27% in 2015 to 51% in 2019 [7].
The aims of this study were to investigate the prevalence of meningococcal carriage in Norwegian teenagers, to characterise circulating carriage strains of N. meningitidis and to identify risk factors for meningococcal carriage in this age group. These data are needed to improve public health recommendations regarding IMD in teenagers and to evaluate if meningococcal vaccines should be implemented in the national immunisation programme.
Materials and methods
Study design and inclusion of participants
Students in lower (grades 8–10) and upper (grades 11–14) secondary schools in three counties in southeastern Norway were recruited in a cross-sectional study. Only schools with >300 students were invited to participate. The study was conducted in two sampling periods, October–November 2018 (county 1) and February–April 2019 (counties 2 and 3), due to laboratory capacity. Students and parents received information about the study through SMS or e-mail from the school administrations. Students aged ≥16 years consented on their own behalf, while parental consent was requested for younger students. The study was approved by the Regional Committee for Medical and Health Research Ethics, Southeast Norway (reference number 2018/465).
Sampling and bacterial identification
Swabbing of the posterior pharyngeal wall and one tonsil was performed on school grounds using a sterile cotton swab (Copan Diagnostics, CA, USA). Samples were plated on site on chocolate agar (with lincomycin 2.0 mg, colistin sulphate 12.0 mg, amphotericin B 2.0 mg and trimethoprim lactate 13.0 mg). Plates were transported to the laboratory in styrofoam containers at room temperature and then incubated for 24–48 h at 35 °C with 5% CO2 within 6 h of collection. Single colonies suspected to represent N. meningitidis were harvested and species identification was confirmed with MALDI-TOF mass spectrometry (Bruker Daltonik GmbH, Bremen, Germany). Meningococcal isolates were stored at −80 °C in Greaves solution.
Molecular characterisation of N. meningitidis isolates
DNA from the confirmed isolates was extracted using MagNA Pure 96 (Roche Life Science, Basel, Switzerland). Whole genome sequencing (WGS) was performed using the MiSeq platform (Illumina Inc., San Diego, CA, USA) as described previously [Reference Bårnes8]. Genogroups, multilocus sequence types, PorA and FetA types were identified using the PubMLST database [Reference Jolley, Bray and Maiden9]. Isolates with deletions or stop codons in the capsule locus were identified as non-groupable (NG). Those lacking the capsule operon were classified as capsule null (cnl).
Phylogenetic analysis
The genomes of the meningococcal carriage isolates were compared to the genomes of all invasive isolates in Norway submitted to the National Reference laboratory between January 2018 and July 2019 (n = 35). Phylogenetic trees were created by using the Neighbour-Joining algorithm on the pairwise distances across the 1605 loci defined in the core genome multilocus sequence typing (cgMLST) scheme v 1.0 [Reference Bratcher10]. In this calculation, incomplete loci were disregarded for the purpose of pairwise distance calculation. For the dominating invasive ccs in the study period, cc11 and cc23, minimum spanning trees were created based on the aforementioned allelic distance profiles and visualised using GrapeTree [Reference Zhou11]. For the evaluation of differences between closely related invasive-carrier clusters, we evaluated all allelic differences across the 3050 loci defined in PubMLST (wgMLST) at the time of analysis (November 2019). Closely related invasive-carrier clusters were defined as a difference of <20 genes.
Assessment of risk factors for meningococcal carriage and definition of variables
The participants were asked to complete an electronic questionnaire that assessed smoking habits, exposure to passive smoking at home, use of Swedish snus (smokeless tobacco), intimate kissing, sharing of drinking bottles, attendance of youth gatherings and parties, participation in the russ celebration, recent throat infection, recent use of antibiotics, parental background, parental education and the number of persons in the household. Using the unique personal identification number assigned to all residents of Norway, data from each participant were linked to their records of meningococcal vaccination in The Norwegian Immunisation Registry, SYSVAK [7]. The participants were considered vaccinated if they had received a meningococcal vaccine more than 2 weeks and less than 5 years before sampling (i.e. the expected duration of protection for MCV4 [Reference Baxter12]). Age was defined as age in years at the time of sampling. In Norway, the attended school grade is based on birth cohort and therefore corresponds well with age.
Study population
Among the 106 schools invited to recruit students, 21 lower and 24 upper secondary schools participated in the study. In total, 2511 students consented and throat swabs were obtained from 2296 participants; 1354 students in October–November 2018 (59.0%) and 942 students in February–April 2019 (41.0%) (Fig. 1). The median number of students with throat swabs per school was 51 (range 10–80). All 2296 participants were included when assessing carriage prevalence. Among the 2296 participants from whom throat swabs were collected, 137 were excluded because of incomplete questionnaire data and/or unknown meningococcal vaccination status and 2159 (93.8%) individuals were included in the analyses of risk factors (Fig. 1).
Statistical methods
We calculated carriage rates with 95% confidence intervals (CIs) based on cluster-robust standard error. When testing whether the carriage rate differed between sampling periods or genders, the Pearson χ 2 test statistic was corrected with the second-order correction of Rao and Scott and converted into an F statistic [Reference Rao and Scott13]. The associations between potential risk factors and the risk of carriage were studied using a logistic model estimating odds ratios (ORs) and 95% CIs. The generalised estimating equations approach was used to fit the model in order to account for dependencies within schools. An exchangeable correlation structure was assumed. We used cluster-robust estimates of the standard errors. The multivariable model included smoking habits, exposure to passive smoking at home, use of Swedish snus, intimate kissing, sharing of drinking bottles, attendance of youth gatherings and parties, participation in the russ celebration, recent throat infection and parental background. For the variables assessing smoking habits and the use of snus, the response categories ‘daily’ and ‘occasionally’ were merged in the multivariable analysis. In addition, we adjusted for sampling period, gender, school grade (with grades 13 and 14 combined) and record of previous vaccination with MCV4. Use of antibiotics, parental education and number of persons in the household were not included in the multivariable analysis due to high degree of missing and/or inadequate/low quality of data. P-values <0.05 were considered statistically significant. All analyses were performed with Stata/SE 15.0 (Stata-Corp, College Station, Texas, USA).
Results
Characteristics of the participants
The median age of the 2296 participants from whom throat swabs were obtained was 16 years (range 12–24 years) and 68.5% attended upper secondary school. In lower secondary schools (grades 8–10), the majority of students (96.3%) were 13–15 years (range 12–16 years) whereas the majority of students (93.1%) in upper secondary schools (grades 11–14) were 16–18 years (range 15–24 years). Overall, 61.5% of participants were female.
Carriage rate of N. meningitidis
In total, 167 of the 2296 participants were identified as carriers of N. meningitidis. The overall carriage rate was 7.3% (95% CI 5.5–9.6%). There was no significant difference in carriage rate between the first and second sampling periods (7.8% (95% CI 5.3–11.2) vs. 6.6% (95% CI 4.3–9.8), P = 0.556), or between boys and girls (7.8% (95% CI 5.1–11.6) vs. 6.9% (95% CI 5.2–9.2), P = 0.540). The majority of the meningococcal isolates (91.6%) were found in students attending upper secondary schools. The carriage rate in lower secondary school was 1.9% (95%CI 0.9–3.9) and in upper secondary school 9.7% (95% CI 7.5–12.6). The highest carriage rate was found among the 18-year-olds (16.4% (95% CI 12.7–21.0%)) (Fig. 2).
Molecular characterisation of carriage isolates
Among the 167 carriage isolates, cnl dominated (40.1%), followed by genogroups Y (33.5%), NG (11.4%), B (9.0%), X (2.4%), C (1.8%) and W (1.8%) (Fig. 3a). In October–November 2018, cnl isolates dominated (45.7%) followed by genogroup Y (21.9%), whereas in February-April 2019 genogroup Y dominated (53.2%) followed by cnl (30.7%) (Fig. 3b and 3c). The distributions of genogroup Y and cnl between sampling periods were not significantly different (P = 0.06 and 0.12 for genogroup Y and cnl, respectively).
Among the 167 carriage isolates, 157 were assigned to 17 previously defined ccs (Table 1). Ten isolates representing six sequence types (STs) were unassigned to a cc. The most common ccs were cc23 (35.9%) and cc198 (32.3%). There were 28 previously defined STs and 11 new STs represented among the 167 isolates; ST-23 (33.5%) and ST-823 (26.4%) dominated. The majority of cc23 were ST-23 genogroup Y (86.7%) and the majority of cc198 were ST-823 cnl (81.5%). These two clones were found in 22 of 45 and 21 of 45 schools, respectively. Genogroups B, C, W and X isolates belonged to multiple STs. Among the 167 isolates, there were 48 PorA types (nine new variants), with P1.5-2, 10-1 (29.3%) and P1.18, 25–44 (22.2%) dominating and 20 FetA types (three new variants) with predominance of F4-1 (34.1%) and F1-43 (18.6%). The fetA gene was missing in 14 isolates (8.4%).
cc, clonal complex; ST, sequence type; PorA, porin A; FetA, ferric enterobactin receptor; NG, non-groupable; cnl, capsule null locus; u.a., unassigned; –, gene not present.
a New variants.
A phylogenetic tree was generated to illustrate the relationships among the 167 carrier isolates and 35 recent invasive Norwegian isolates (Fig. 4). Of the carriage isolates, 38.9% were both capsulated and belonged to ccs presently causing invasive disease in Norway. Isolates of cc23 had a median allelic distance of 120.0 (range 4.0–732.0, s.d. 142.2) and cc198 isolates a median distance of 149.0 (range 3.0–827.0, s.d. 169.6), showing that the latter group was slightly more diverse. There was no apparent association between cc and school number (P = 0.49), county (P = 0.31) or vaccination status (P = 0.49) for the carriage isolates, as evaluated using the χ 2-test. The three C:P1.5,2:F3-3:ST-11 carriage isolates were very closely related to the invasive isolate with the same composition, differing in less than 20 genes in wgMLST (Fig. 5, panel a, upper right corner). The cc23 carriage and invasive isolates appeared more genetically distant (Fig. 5, panel b), although six carriage isolates also differed by less than 20 genes from an invasive isolate.
Risk factors associated with carriage of N. meningitidis
Of the 2159 participants assessed for risk factors, 157 were carriers of N. meningitidis. In general, upper secondary school students reported higher exposure to risk factors than lower secondary school students (Table 2). Occasional cigarette smoking was reported by 18.4% of the upper secondary school students, but only 2.3% of lower secondary school students. Daily smoking was rare in both upper and lower secondary schools students (⩽0.5%). In upper secondary school student, 11.8% used Swedish snus daily and 9.3% occasionally, while in lower secondary school students, 0.6% used Swedish snus daily and 1.3% occasionally. Less than 5% of participants had used e-cigarettes (vaping) or waterpipe, of which few were daily users (⩽0.6% and ⩽0.1%, respectively). While 47.1% of upper secondary school students participated in the russ celebration, only 1.2% of lower secondary school students did. Vaccination with MCV4 was registered for 0.6% and 25.1% of lower and upper secondary school students, respectively.
a Grade 14 comprises students with vocational education taking supplementary studies to qualify for higher education
b Variables not included in multivariable analysis
Both use of Swedish snus (OR 1.56, 95% CI 1.07–2.27) and participation in the russ celebration (OR 2.85, 95% CI 1.62–5.02) were associated with a higher risk of carriage of N. meningitidis (Table 3). We also observed a positive association between carriage and the numbers of persons kissed as well as the number of times attending parties or youth gatherings. Compared to not having kissed, ORs for number of persons kissed the last month were 2.73 (95% CI 1.58–4.70) for one person, 3.03 (95% CI 1.44–6.36) for two persons and 2.76 (95% CI 1.49–5.10) for more than two persons. Compared to not having attended parties or youth gatherings, ORs for attendance the last three months were 2.14 (95% CI 1.01–4.54) for 1–3 times, 2.20 (95% CI 1.03–4.71) for 4–6 times, 2.84 (95% CI 1.27–6.36) for 7–10 times and 3.50 (95% CI 1.45–8.48) for >10 times. There was no significant association between carriage and sampling period, gender, throat pain or upper respiratory infection, active or passive cigarette smoking, use of waterpipe or e-cigarettes, sharing drinking bottles or parental background. Carriage rate increased with school grade (hence age) in the univariate analysis, but not when adjusting for other risk factors. Previous vaccination with MCV4 was not associated with carriage, neither carriage of all genogroups (Table 3) nor of vaccine-specific genogroups (OR 1.00; 95% CI 0.53–1.88).
OR, odds ratio; CI, confidence interval.
a Grade 14 comprises students with vocational education taking supplementary studies to qualify for higher education.
* P values <0.05 are considered statistically significant.
Discussion
The overall carriage rate of N. meningitidis was 7.3%. Meningococcal carriage was associated with the use of Swedish snus, intimate kissing, attending parties or youth gatherings and participation in the russ celebration. Age, cigarette smoking and MCV4 vaccination was not associated with carriage. Isolates belonging to cnl cc198 and genogroup Y cc23 dominated. Almost 40% of the carriage isolates were similar to invasive isolates currently causing IMD in Norway.
The main strengths of this study were that over 90% of the throat swabs were collected by one study team member, plated on site and analysed at the meningococcal reference laboratory, which probably increased the quality of the samples. The study also had potential limitations. A higher proportion of upper secondary school students were included since participants over 16 years age, not needing parental consent, could be recruited on site. This may have overestimated carriage rate in the overall study population. Assessment of risk factors might have been affected by information bias as students who completed questionnaires in the company of parents or peers could have been influenced to either under- or over-report risk-behaviour.
Meningococcal carriage has decreased in Norway in the last 30 years, from 28.3% in 15–24-year-olds in 1991 [Reference Caugant14] to 7.3% in the present study. Studies in Europe have shown variable carriage rates in the past decade. Rates below 10% was found in an ongoing study in Swedish university students [Reference Stenmark15], in Italy in 14–19-year-olds 2012–2013 [Reference Neri16] and in 14–21-year-olds in 2016 [Reference Terranova17] and in Turkey in 10–24-year-olds in 2015 [Reference Tekin18]. Higher carriage rates than in our study have been observed in the UK in 18–19-year-olds in 2015–2016 (14–46%) [Reference Oldfield19] and in the Netherlands in 13–23-year-olds in 2013–2014 (16%) [Reference Van Ravenhorst20]. The latter two studies included university students, which could explain the higher carriage rates. Despite increasing carriage prevalence with age in our study, age was not significantly associated with carriage when adjusting for other risk factors. Risk-behaviour rather than age influenced the risk of carriage, in concordance with the Dutch study [Reference Van Ravenhorst20].
To our knowledge, our study is the first to report the use of Swedish snus as a possible risk factor for meningococcal carriage. Swedish snus is a form of moist and ground smokeless tobacco put under the upper lip either in loose powder or portioned in small cellulose pouches. Swedish snus is prohibited in all EU/EEA countries except Sweden and Norway. However, Swedish snus is also marketed in the USA and 80% of American 12–17-year-olds reported to have tried snus pouches in 2013–2015 [Reference Ambrose21]. Use of Swedish snus has increased over the last 20 years in Norway and Sweden, and 25% of 16–24-year-olds in Norway used Swedish snus daily or occasionally in 2018 [22]. This corresponds to data from 16 to 24-year-olds in our study. Nicotine plasma levels have been reported to be higher after using Swedish snus compared to cigarettes and nicotine gum [Reference Digard23]. Nicotine is believed to have immunomodulatory effects and may suppress antibody production [Reference Geng24]. Moreover, nicotine increases biofilm formation of common oral bacteria [Reference El-Ezmerli and Gregory25]. Both capsulated and uncapsulated strains of N. meningitidis have shown to produce biofilm in human bronchial epithelial cells [Reference Neil, Shao and Apicella26]. A combination of more optimal conditions for bacterial colonisation through induction of biofilm and a weaker mucosal immune defence, might explain the increased risk of meningococcal carriage in people using Swedish snus.
Participation in the russ celebration increased the risk of carriage almost threefold in our study. Bacterial transmission is probably high during this event through overcrowding, sharing of bottles, intimate kissing and excessive alcohol consumption [Reference Taylor6]. IMD cases are annually linked to the russ celebration in Norway. Similarly, first year university students in the UK and USA, with similar risk behaviour as in the russ celebration, also have a higher risk of meningococcal carriage and IMD [Reference Maclennan4, Reference Breakwell27]. High alcohol consumption has been associated with carriage [Reference Imrey28]. This relationship was not assessed in our study, but needs further evaluation.
A relationship between cigarette smoking and carriage has previously been reported in Norway [Reference Caugant14], in a carriage study in the UK from 1999 [Reference Maclennan4], as well as in a recent study in South Australia from 2019 [Reference Marshall29]. In contrast, such a relationship was not found in the present study, or in university students in USA in 1992–1993 [Reference Imrey28]. The studies from Norway and the UK were conducted in the 1990s, before the European Council recommended comprehensive smoke-free laws in EU member states and when smoking rates were high [30]. The low level of smoking in our study could explain the lack of an association between carriage and smoking.
Along with an increasing uptake of MCV4 among the graduating students involved in the russ celebration, the incidence of IMD in this group has decreased in the last decade. However, vaccination with MCV4 did not have an impact on carriage in our study. The effects of MCV4 on carriage and herd immunity are uncertain [Reference Balmer31] and needs further assessment.
The dominant clones among the 167 isolates of N. meningitidis were observed in many of the schools and in all three counties. The majority of the isolates were cnl or NG, in accordance with recent studies in adolescents and young adults in Europe [Reference Stenmark15, Reference Neri16, Reference Oldfield19, Reference Van Ravenhorst20, Reference Jeppesen32]. In our study, cnl cc198 dominated, as recently reported in a study among Italian teenagers [Reference Neri16] and in the ongoing study in Swedish university students [Reference Stenmark15]. The lack of capsular expression is thought to increase the meningococcus' ability to colonise the nasopharyngeal mucosa [Reference Tzeng, Thomas and Stephens33]. Cnl strains are usually apathogenic since the meningococcal capsule is a major virulence factor [Reference Vogel and Frosch34]. In rare occasions, however, cnl cc198 strains have caused IMD in immunocompetent individuals [Reference Hoang35–Reference Johswich37].
Genogroup Y cc23 dominated among the isolates expressing capsular genes. Similar findings were reported in a previous carriage study in Norway among 13–14-year-olds in 1989 [Reference Rosenqvist, Conde-Glez, Morse, Rice, Sparling and Calderón38] and from the ongoing carriage study in Swedish university students [Reference Stenmark15]. Carriage of genogroup Y has been rare in adolescents and young adults in other European countries the last years [Reference Neri16, Reference Tekin18–Reference Van Ravenhorst20]. IMD caused by genogroup Y has increased in several European countries in the last decades [Reference Ladhani39–Reference Bijlsma42]. This genogroup previously caused IMD mainly in adults [Reference Hellenbrand40], but it has also affected adolescents and young adults the last decade [Reference Säll41, Reference Bröker43]. In 2018, 50% of the overall IMD cases in Norway and five of seven cases in teenagers were caused by genogroup Y cc23.
Whereas genogroup B has dominated among isolates expressing capsular genes in carriage studies among adolescents and young adults in Europe the last 10 years [Reference Oldfield19, Reference Van Ravenhorst20, Reference Tryfinopoulou44], only 10% of the carriage isolates in our study were genogroup B. Correspondingly, there has not been any genogroup B IMD cases in teenagers in Norway in the previous 5 years. The relatively high circulation of serogroup B in Norway in the late 1990s [Reference Yazdankhah45] might have led to a higher level of natural immunity in the population, which may explain our findings.
We found four ST-11 carriage isolates in our study. Analyses of the three C:P1.5,2:F3-3:ST-11 isolates using cgMLST revealed a close similarity to a Norwegian invasive isolate from the same time period in a patient residing in the same county. After an increase in IMD caused by genogroup C ST-11 in Europe in the late 1990s, implementation of meningococcal C conjugate vaccines in many immunisation programmes across Europe have been effective in reducing the incidence [Reference Whittaker46]. Genogroup C IMD has been rare in Norway in the last 5 years, even though meningococcal vaccines are not included in the national immunisation programme.
Only one of the three genogroup W isolates represented the hypervirulent clone W:P1.5,2:F1-1:ST11. This clone has been seen among Norwegian invasive isolates the last years. A long-term increase such as seen in West Africa [Reference Koumare47] and South America [Reference Abad48] in the early 2000s and in some European countries in recent years [Reference Krone49] has not been observed in Norway.
Genogroup X isolates, which has also been reported among teenagers in other European countries [Reference Neri16, Reference Van Ravenhorst20, Reference Jeppesen32, Reference Tryfinopoulou44], were found in three individuals in our study. IMD caused by genogroup X is uncommon in Western Europe [Reference Whittaker46], but an increase has been observed in Eastern Europe recently [Reference Bai50]. Differences between countries and fluctuations in genogroup distribution may be explained by emergence of hypervirulent clones, natural immunity, exposure to different environments and behavioural risk factors, as well as variations in vaccination policy.
Conclusion
Carriage of N. meningitidis was 7.3% with a peak of 16.4% in 18-year-olds and is probably due to age-related risk behaviour, supporting current national recommendations regarding meningococcal vaccination of teenagers. The finding that Swedish snus might be a new risk factor for meningococcal carriage needs further investigation. Even though most of the circulating isolates of N. meningitidis lacked capsule-expressing genes suggesting a low risk for development of IMD, more than one-third of the circulating isolates have invasive potential. Analyses of the degree of natural protection as well as health economic evaluations are necessary to evaluate the demand for a meningococcal vaccination programme for teenagers in Norway.
Acknowledgements
We thank the students who contributed to this study and their parents, the school administrations and the school nurses. We also thank the clinical study team Anne-Cathrine Kristoffersen, Jeanette Stålcrantz, Inger Lise Haugen, Jorunn Hill, Mark Woods, Jacob Dag Berild, Are Berg and the meningococcal reference laboratory with Gunnhild Rødal, Inger Marie Saga, Jan Oksnes, Anusha Balasingham, Nadia Debech, Ingerid Kirkeleite, Jill Andersen and Torbjørn Bruvik. This work was supported by the Norwegian Institute of Public Health, Oslo, Norway. This work was performed on the TSD (Tjeneste for Sensitive Data) facilities, owned by the University of Oslo, operated and developed by the TSD service group at the University of Oslo, IT-Department (tsd- [email protected]). This publication made use of the Neisseria Multi Locus Sequence Typing website (https://pubmlst.org/Neisseria/) sited at the University of Oxford [Reference Jolley, Bray and Maiden9]. The development of this site has been funded by the Wellcome Trust and European Union. This research received no specific grant from any funding agency, commercial or not-for-profit sectors
Author contributions
SVW, LMN, DAC and GT initiated and designed the study. GT, LMN and SVW participated in the sampling. IL, SVW and GT performed the statistical analyses. TB was head of data management. OBB performed phylogenetic analyses. SVW drafted the manuscript. All authors contributed to the interpretation of the data, to writing and revising the manuscript and to approval of the final manuscript.
Conflict of interest
The authors declare that there are no conflicts of interest.