# Development of a PCR algorithm to detect and characterize Neisseria meningitidis carriage isolates in the African meningitis belt

Kanny Diallo, Mamadou D. Coulibaly, Lisa S. Rebbetts, Odile B. Harrison, Jay Lucidarme, Kadidja Gamougam, Yenenesh K. Tekletsion, Akalifa Bugri, Aliou Toure, Bassira Issaka, et al.

► To cite this version:

Kanny Diallo, Mamadou D. Coulibaly, Lisa S. Rebbetts, Odile B. Harrison, Jay Lucidarme, et al. Development of a PCR algorithm to detect and characterize Neisseria meningitidis carriage isolates in the African meningitis belt. PLoS ONE, 2018, 13 (12), pp.e0206453. 10.1371/journal.pone.0206453. pasteur-03238217

HAL Id: pasteur-03238217

https://riip.hal.science/pasteur-03238217

Submitted on 27 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Development of a PCR algorithm to detect and characterize *Neisseria meningitidis* carriage isolates in the African meningitis belt

Kanny Diallo¹,², Mamadou D. Coulibaly¹, Lisa S. Rebbetts², Odile B. Harrison², Jay Lucidarme³, Kadidja Gamougam⁴, Yenenesh K. Tekletsion⁵, Akalifa Bugri⁶†, Aliou Toure¹, Bassira Issaka⁷, Marietou Dieng⁸, Caroline Trotter⁹, Jean-Marc Collard⁷, Samba O. Sow¹, Xin Wang¹⁰, Leonard W. Mayer¹⁰, Ray Borrow³, Brian M. Greenwood¹¹, Martin C. J. Maiden², Olivier Manigart[11*], for the MenAfriCar Consortium¹¹

¹ Centre pour le Développement des Vaccins (CVD), Bamako, Mali, ² University of Oxford (Department of Zoology), Oxford, United Kingdom, ³ Public Health England, (PHE–Vaccine Evaluation Unit), Manchester, United Kingdom, ⁴ Centre de Support en Santé Internationale (CSSI), Ndjamena, Chad, ⁵ Armauer Hansen Research Institute (AHRI), Addis Ababa, Ethiopia, ⁶ Navrongo Health Research Centre (NHRC), Navrongo, Ghana, ⁷ Centre de Recherche Médicale et Sanitaire (CERMES), Niamey, Niger, ⁸ Institut de Recherche pour le Développement (IRD), Dakar, Senegal, ⁹ University of Cambridge (Disease Dynamics Unit -Department of Veterinary Medicine), Cambridge, United Kingdom, ¹⁰ Centers for Disease Control and Prevention, Division of Bacterial Diseases, Atlanta, United States of America, ¹¹ London School of Hygiene and Tropical Medicine (LSHTM), London, United Kingdom

† Deceased.

¶ Membership of the MenAfriCar Consortium is provided in the Acknowledgments

* olivier.manigart@lshtm.ac.uk

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Diallo K, Coulibaly MD, Rebbetts LS, Harrison OB, Lucidarme J, Gamougam K, et al. (2018) Development of a PCR algorithm to detect and characterize *Neisseria meningitidis* carriage isolates in the African meningitis belt. PLoS ONE 13(12): e0206453. https://doi.org/10.1371/journal.pone.0206453

**Editor:** Ulrich Melcher, Oklahoma State University, UNITED STATES

**Received:** February 13, 2018

**Accepted:** October 12, 2018

**Published:** December 5, 2018

**Copyright:** This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was funded by the Meningitis Research Foundation (www.meningitis.org) under the auspices of the MenAfriCar project. The MenAfriCar Consortium was supported by grants from the Wellcome Trust (086546) and the Bill & Melinda Gates Foundation (OPP51251). The funders had no role in study design, data collection

## Abstract

Improved methods for the detection and characterization of carried Neisseria meningitidis isolates are needed. We evaluated a multiplex PCR algorithm for the detection of a variety of carriage strains in the meningitis belt. To further improve the sensitivity and specificity of the existing PCR assays, primers for gel-based PCR assays (sodC, H, Z) and primers/probe for real-time quantitative PCR (qPCR) assays (porA, cnl, sodC, H, E, Z) were modified or created using Primer Express software. Optimized multiplex PCR assays were tested on 247 well-characterised carriage isolates from six countries of the African meningitis belt. The PCR algorithm developed enabled the detection of N. meningitidis species using gel-based and real-time multiplex PCR targeting porA, sodC, cnl and characterization of capsule genes through sequential multiplex PCR assays for genogroups (A, W, X, then B, C, Y and finally H, E and Z). Targeting both porA and sodC genes together allowed the detection of meningococci with a sensitivity of 96% and 89% and a specificity of 78% and 67%, for qPCR and gel-based PCR respectively. The sensitivity and specificity ranges for capsular genogrouping of N. meningitidis are 67%-100% and 98%-100% respectively for gel-based PCR and 90%-100% and 99%-100% for qPCR. We developed a PCR algorithm that allows simple, rapid and systematic detection and characterisation of most major and minor N. meningitidis capsular groups, including uncommon capsular groups (H, E, Z).
and analysis, decision to publish, or preparation of
the manuscript.

**Competing interests:** RB reports performing contract research on behalf of Public Health England for Pfizer, Sanofi Pasteur, Serum Institute of India, and GSK. CLT and OM report receiving a consulting payment from GSK. All other authors report no potential conflicts. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

**Introduction**

*Neisseria meningitidis* has been responsible for numerous outbreaks of septicaemia and menin-
gitis in countries of the African Sahel and sub-Sahel, a region known as the “African meningi-
tis belt”[1]. The devastating nature of this disease, which predominately affects children and
young adults, makes its control an international health priority. *N. meningitidis* is a common
commensal bacterium capable of harmlessly colonizing the upper respiratory tract in a
dynamic process referred to as “carriage” [2]. Mucosal carriage is a necessary prerequisite for
meningococcal disease. Carriage of *N. meningitidis* causes disease when it breaches the host
epithelial barriers and invades the bloodstream or penetrates the subarachnoidal cavity and
infects the meninges [3].

*N. meningitidis* are classified into 12 serogroups (A, B, C, E, H, I, K, L, W, X, Y and Z), based on the structure and chemical composition of their capsular polysaccharide [4]. The serogroups associated with disease in the African meningitis belt are A, W, X, Y and C. With serogroup A epidemics being eliminated following the introduction of the serogroup A conjugate vaccine MenAfriVac[5-7], non-A meningococcal serogroups become the major cause of meningococcal disease in sub Sahelian Africa. Recently, large serogroup C epidemics occurred in Niger and Nigeria, [8] demonstrating the epidemic potential of non- serogroup A meningococci. Timely control and prevention of the epidemics have been a challenge due to vaccine shortages, which has resulted in tremendous pressure on hospital facilities and public health authorities. Serogroup B, a common cause of meningococcal disease in many high-income countries [9], is rarely detected in the meningitis belt. Serogroups E, I, H, K, L, Z are less commonly identified in invasive disease. Queries on the more than 38,633 records on PubMLST. org/neisseria website have demonstrated that serogroups E and Z are found both in invasive and carriage specimens, serogroups K, L, H are only found in carriage samples but very rarely, and serogroup I was not found at the time of writing. Furthermore, previous carriage studies in Europe have identified a high proportion of *N. meningitidis* isolates, among which, up to 30-45% did not express their capsule genes [10] and 16-20% lacked regions A, B and C of the capsular polysaccharide synthesis (*cps*) genes, known as a capsule null (*cnl*) genotype [11].

Several meningococcal cross-sectional carriage studies have been conducted in the African meningitis belt in recent years, including surveys in Burkina Faso [12, 13], and in seven other countries (Chad, Ethiopia, Ghana, Mali, Niger, Nigeria and Senegal) as part of the MenAfriCar project [14, 15]. Conventional microbiological methods have been used in these studies to detect and characterise meningococci, but these methods are time consuming, expensive and susceptible to errors. Consequently, many African laboratories have been increasingly adopt-
ing molecular biology tools to detect and characterise *N. meningitidis* and its capsular geno-
types as they are more sensitive, specific, and less time consuming than bacteriological
methods. DNA-based methods such as Polymerase Chain Reaction (PCR) have greatly
improved laboratory diagnosis [16]; many studies have developed reliable PCR assays for the
detection of *N. meningitidis* [17–19] and characterization of its capsular genogroup using spe-
cific genes. However, most studies have focused on the common disease serogroups, A, B, C,
W, X and Y [20–23]. Fewer studies have investigated PCR for the detection of *cnl* or ser-
ogroups E, H, and Z which are mostly found in carriage [20, 24]. Some research groups have
developed gel-based PCR assays for the detection and genogrouping of *N. meningitidis*
whereas others have used real-time quantitative PCR (qPCR) protocols. The choice of one
technique over another is often dictated by budgetary considerations, as qPCR is more expen-
sive but less subject to contamination. Evaluation of the impact of meningococcal conjugate
vaccines requires sensitive and specific PCR assays that cover all serogroups to accurately and
thoroughly assess the epidemiology of *N. meningitidis* carriage and disease, including the
potential emergence of less common groups through strain replacement and/or capsule switching [25].

The present study was conducted within the MenAfriCar project, which collected pharyngeal swabs from over 60,000 subjects in seven countries within the African meningitis belt [15]. Our objectives were to improve *N. meningitidis* carriage assessment by developing and optimising a simple, accurate and standardised multiplex PCR algorithm for the detection of *N. meningitidis* species and its capsular genogroups including: (i) the most prevalent genogroups A, B, C, X, Y and W (ii) carriage-associated genogroups such as E, H, and Z; and (iii) unencapsulated capsule null meningococci. Both traditional gel-based PCR and qPCR that can be adapted to African laboratories were validated using carriage isolates collected as part of MenAfriCar studies.

## Material and methods

### Literature search for available primers/probes

A literature search was undertaken to identify PCR methods with primers and probes published on *N. meningitidis* detection and characterization. The NCBI PubMed database was used with search terms “*Neisseria meningitidis*” and “PCR” (date of access: November 2011). The different techniques and detection algorithms were evaluated in order to define the most appropriate ones for our context. A focus was placed on summarising the best gel-based and qPCR assays already existing for detection and characterization of meningococci.

### In silico specificity screen—Evaluation of selected primers and probes

A total of 112 primers and 16 probes reported in the literature were chosen for further analysis. Firstly, nucleotide mismatches were identified through the comparison with consensus sequences derived from target genes obtained from the pubMLST.org/neisseria website (accessed on September 2012). Primer and probe sequences were compared to all available allelic variants of the target genes to determine whether they were conserved in the majority of alleles available at that time. Technical checks on all primers and probes were then undertaken using the *Primer Express* software design guidelines (Applied Biosystem International (ABI), Villon sur Yvette, France): melting point temperatures (Tm) and potential secondary structures were assessed with optimal Tm varying between 58–60°C and 68–70°C for primers and probes respectively. Finally, nucleotide BLAST searches were undertaken using the NCBI BLAST tool to identify possible cross-reactions with other bacteria [26]. Primers and probes with more than 3 base pair (bp) differences compared to the consensus sequences were not used. Secondary structures or self-dimers were minimised, to avoid their potential detrimental effect on the efficiency of the reaction, but their presence was not a systematic criterion for rejection. Sequences were retained if they displayed 100% coverage according to the NCBI BLAST tool. When the result of these evaluations was unsatisfactory, new primers and/or probes were created.

### Positive controls

Positive controls consisted of purified DNA from genogroups: A (isolate name: Z2491); W (A22); X (860060); Y (71/94); B (EG329, BZ198); C (L93/4286); Z (0084/93); E (297–0); H (H/ASH/87); and capsule null (0083/93) [27] as well as two genogroup E isolates (M07746, M22445). Reference strains for genogroups A and X, were purchased from the American Type Culture Collection (ATCC). Another confirmed capsule null isolate from the first Malian
MenAfriCar cross-sectional study (sample 01-010664-XS1) (18) was used as a *cnl* positive control.

## Carriage isolates

*N. meningitidis* isolates were collected during the first MenAfriCar cross-sectional studies as previously described [15]. Briefly, pharyngeal swabs were obtained from participants in the rural and urban sites of six centres in Chad, Ethiopia, Ghana, Mali, Niger and Senegal. Suspected *N. meningitidis* were subjected to oxidase test, Gram stain, biochemical tests and sero-agglutination [28]. Boiled cell suspensions to release high concentration of DNA were prepared in each center as part of the MenAfriCar routine protocol from bacteria grown overnight on blood agar plates [15]. Extracted DNA samples from all oxidase positive bacteria were also sent to the University of Oxford [14] for molecular characterization. A subset of characterized samples was chosen to constitute a representative panel of the most prevalent *N. meningitidis* genogroups as well as non *N. meningitidis* from each country. Next, the selected boiled cell suspension extracted DNA samples were sent to the Center for Vaccine Development in Mali (CVD-Mali) where PCR validation was conducted. A total of 247 extracted DNA samples were received, with 42, 38, 18, 46, 44 and 59 samples from Chad, Ethiopia, Ghana, Niger, Mali and Senegal, respectively.

## Characterisation of *N. meningitidis* at the University of Oxford, UK

*Neisseria* species were analysed using the *rplF* assay [29]. All samples were also tested with a capsule null (*cnl*) sequencing assay [11] or sequencing of a fragment of the *ctrA* gene allowing the characterisation of genogroup A, E, K, L, X and Z. Samples that were positive for *N. meningitidis* and negative for *cnl* were tested by a qPCR multiplex assay to determine genogroups A, B, C, W, X and Y [23]. These results were considered as the reference or gold standard for this study since the assays described above were already validated and published.

## Novel gel-based PCR assay design and optimization in Bamako, Mali

New primers were designed for *sodC*, genogroup Z, and H using the online tools Primer 3 and Primer quest (Table 1)[30, 31]. A ready-to-load Master Mix (rlMMx) from Solis Biodyne, stable at room temperature and therefore very convenient for transport and testing in Africa, was employed. rlMMx with different concentrations of MgCl₂ were tested to find the optimal concentration for the assay. Thermal cycling conditions were the same for all tests: 1 cycle of 3 minutes (min) at 95°C; 35 cycles of 30 seconds (s) at 95°C; 30 s at 55°C; 30 s at 72°C. This was followed by 1 cycle of 10 min at 72°C and a final conservation step at 4°C.

## Novel quantitative real-time PCR (qPCR) assay design and optimization in Bamako, Mali

Although our assay was developed to determine whether a specimen is positive or negative for meningococci and their genogroups, we keep the term “qPCR” to comply with the MIQE guidelines [32]. Quantification of meningococci would not add any information in our context since we were working with DNA samples extracted from bacterial culture. Primers and probes were designed for the *N. meningitidis* species-specific gene *porA* and for the non-species-specific *cnl*, and a modified species-specific *sodC* probe was designed based on a previously published method [18]. Genogroups were determined using the 2^nd^ (A, W, X) and 3^rd^ (B, C, Y) qPCR multiplexes developed by Wang *et al.* [23]. Primers and probes for these genogroups were kindly provided by the Centers for Disease Control and Prevention (CDC).
<table><caption>Table 1. Primer sequences for gel-based PCR.</caption><thead><tr><th>Gene</th><th>Primer name</th><th>Sequence</th><th>Start/end</th><th>Size</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="2">porA</td><td>porA 2F</td><td>GCG GTT TTG CCG GGA ACT AT</td><td rowspan="2">755-1017</td><td rowspan="2">251</td><td rowspan="2">Bennett, DE, 2006</td></tr><tr><td>porA 15R</td><td>AGT GGC GGC AAT TTC GGT CGT ACT</td></tr><tr><td rowspan="2">galE-cnl-tex</td><td>GH26R (cnl aF)</td><td>GGT CGT CTG AAA GCT TGC CTT GCT C</td><td>250-226</td><td rowspan="2">432</td><td rowspan="2">Claus, H, 2002</td></tr><tr><td>HC344 (cnl aR)</td><td>GGA TTG GAC GAG CGA GAC</td><td>72-55</td></tr><tr><td rowspan="2">sodC</td><td>sodC-F2</td><td>GCG GTT AGT GCA GTA TGT TCA G</td><td rowspan="2">113-641</td><td rowspan="2">530</td><td rowspan="2">Present study</td></tr><tr><td>sodC-R2</td><td>TAA TCA CGC CAC ATG CCA TA</td></tr><tr><td rowspan="2">mynB/csaB</td><td>98-28 (AF)</td><td>CGC AAT AGG TGT ATA TAT TCT TCC</td><td rowspan="2">270-664</td><td rowspan="2">395</td><td rowspan="2">Taha, MK, 2000</td></tr><tr><td>98-29 (AR)</td><td>CGT AAT AGT TTC GTA TGC CTT CTT</td></tr><tr><td rowspan="2">siaD (csw)</td><td>98-32 (W135F)</td><td>CAG AAA GTG AGG GAT TTC CAT A</td><td rowspan="2">911-1030</td><td rowspan="2">120</td><td rowspan="2">Taha, MK, 2000</td></tr><tr><td>98-33 (W135R)</td><td>CAC AAC CAT TTT CAT TAT AGT TAC TGT</td></tr><tr><td rowspan="2">Xcba (csxA)</td><td>X-10 (F)</td><td>ACA GCC CAT AAA ACA CCC GTA TCA TC</td><td rowspan="2">142-343</td><td rowspan="2">202</td><td rowspan="2">Yaro, S, 2012</td></tr><tr><td>X-11 (R)</td><td>GTG ATT GGA ATC TTC CAA TAT CGG T</td></tr><tr><td rowspan="2">siaD (csb)</td><td>98-20 (BF)</td><td>GCA TGC TGG AGG AAT AAG CAT TAA</td><td rowspan="2">645-1099</td><td rowspan="2">455</td><td rowspan="2">Taha, MK, 2000</td></tr><tr><td>98-19 (BR)</td><td>GGA TCA TTT CAG TGT TTT CCA CCA</td></tr><tr><td rowspan="2">siaD (csc)</td><td>98-17 (CF)</td><td>TCA AAT GAG TTT GCG AAT AGA AGG T</td><td rowspan="2">789-1054</td><td rowspan="2">266</td><td rowspan="2">Taha, MK, 2000</td></tr><tr><td>98-18 (CR)</td><td>CAA TCA CGA TTT GCC CAA TTG AC</td></tr><tr><td rowspan="2">siaD (csy)</td><td>98-34 (YF)</td><td>CTC AAA GCG AAG GCT TTG GTT A</td><td rowspan="2">911-1030</td><td rowspan="2">120</td><td rowspan="2">Taha, MK, 2000</td></tr><tr><td>98-35 (YR)</td><td>CTG AAG CGT TTT CAT TAT AAT TGC TAA</td></tr><tr><td rowspan="2">cshC</td><td>cshC-F1</td><td>GTG CCG ATA TTG CCT CAG AT</td><td rowspan="2">2378-2533</td><td rowspan="2">156</td><td rowspan="2">Present study</td></tr><tr><td>cshC-R1</td><td>CTT CGG ATG GGA ACT TGA AA</td></tr><tr><td rowspan="2">cap29EH</td><td>39429 (29E aF)</td><td>TTG GCG GTT GAA ACC TTA C</td><td rowspan="2">22-715</td><td rowspan="2">694</td><td rowspan="2">Zhu, H, 2011</td></tr><tr><td>39430 (29E aR)</td><td>GCG TAT CAT GCT CCA TTA CCA</td></tr><tr><td rowspan="2">cszC</td><td>cszC-Fwd (ZF)</td><td>AGG TTC ATC TGC TGG GAT TAC GCT</td><td rowspan="2">3506-3806</td><td rowspan="2">301</td><td rowspan="2">Present study</td></tr><tr><td>cszC-Rev (ZR)</td><td>AAG CGA TTA ATG GCC TGT TGC TGG</td></tr></tbody></table>

Sequences, positions and amplicons size of each primer pair used in the 4 multiplexes. In bold are base pairs that showed a mismatch with at least one of the alleles present in pubMLST, but not with the majority of them.

https://doi.org/10.1371/journal.pone.0206453.t001

Primers and probes were also designed for genogroups E, H, and Z using the Primer Express software from ABI. Minor Groove Binding (MGB) hydrolysis probes were designed, as opposed to Black Hole Quencher (BHQ) hydrolysis probes which are more difficult to manufacture. The *sodC* probe was modified in Primer Express to make it 12 base pairs shorter and accommodate an MGB quencher (Table 2). All primers/probe were tested at the following concentrations in a 25 µl solution: 900nM/250nM; 600nM/200nM; 300nM/100nM and 100nM/50nM respectively. Optimal concentrations in a TaqMan Gene Expression MasterMix (reference: 4369016 from ABI) were identified for each primers/probe set and used subsequently in the different assays. Cycling conditions were similar for all tests: 1 cycle of 2 min at 50°C, 1 cycle of 10 min at 95°C followed by 50 cycles of 15 s at 95°C and 1 min at 60°C. The ABI 7500 fast cycler was used to perform the reaction and the results analysed using the 7500 Fast software (see dx.doi.org/10.17504/protocols.io.rjsd4ne for detailed information).

## Study design and statistical analyses

Initially, 16 positive controls, representing all target genogroups (A, B, C, W, X, Y, H, E, Z, cnl), were used for the optimization of both gel-based and qPCR in monoplex and appropriate multiplexes. Optimal qualitative results from monoplex and multiplex assays were evaluated for gel-based PCR by visualizing the quality and strength of the bands, whereas optimal
<table><caption>Table 2. Primer and probe sequences for qPCR.</caption><thead><tr><th>Gene</th><th>Primer/Probe name</th><th>Sequence</th><th>Start/end</th><th>Concentration (nM)</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="3">porA</td><td>porA_fwd_1</td><td>GCC GGC GTT GAT TAT GAT TT</td><td>1126-1145</td><td>600</td><td rowspan="3">Present study</td></tr><tr><td>porA_rev_1</td><td>AGT TGC CGA TGC CGG TAT T</td><td>1210-1192</td><td>600</td></tr><tr><td>porA_Pb_1 FAM (MGB)</td><td>CTT CCG CCA TCG TGT C</td><td>1157-1172</td><td>200</td></tr><tr><td rowspan="3">sodC</td><td>Nm sodC Fwd 351</td><td>GCA CAC TTA GGT GAT TTA CCT GCA T</td><td>446-470</td><td>600</td><td rowspan="3">Dolan, T, 2011</td></tr><tr><td>Nm sodC Rev 478</td><td>CCA CCC GTG TGG ATC ATA ATA GA</td><td>551-573</td><td>600</td></tr><tr><td>sodC Pb387 NED (MGB)</td><td>CAT GAT GGC ACA GCA A</td><td>482-497</td><td>100</td></tr><tr><td rowspan="4">galE/cnl/tex</td><td>cnl_fwd_2</td><td>GAA TTG CAT AGG TTA TCC AAA ATC AC</td><td>110-85</td><td>900</td><td rowspan="4">Present study</td></tr><tr><td>cnl_fwd_2b</td><td>GAG TTG CAC AGA TTA TCC AGA ATC AC</td><td>110-85</td><td>900</td></tr><tr><td>cnl_rev_4</td><td>TTT GCC CGA TAC AAT CTG AAA G</td><td>85-106</td><td>900</td></tr><tr><td>cnl_Pb_2 VIC (MGB)</td><td>ATA AAA CCG GTG CCG CC</td><td>38-22</td><td>250</td></tr><tr><td rowspan="3">csw/synG</td><td>F857</td><td>TAT TTA TGG AAG GCA TGG TGT ATG</td><td>935-958</td><td>600</td><td rowspan="3">Wang, X, 2011</td></tr><tr><td>R964</td><td>TTG CCA TTC CAG AAA TAT CAC C</td><td>1063-1042</td><td>600</td></tr><tr><td>Pb907i FAM (BHQ)</td><td>AAA TAT GGA GCG AAT GAT TAC AGT AAC TAT AAT GAA</td><td>985-1020</td><td>200</td></tr><tr><td rowspan="3">csaB/sacB</td><td>F2531</td><td>AAA ATT CAA TGG GTA TAT CAC GAA GA</td><td>934-959</td><td>900</td><td rowspan="3">Mothershed, EA, 2004</td></tr><tr><td>R2624</td><td>ATA TGG TGC AAG CTG GTT TCA ATA G</td><td>1025-1001</td><td>900</td></tr><tr><td>Pb2591i HEX (BHQ)</td><td>CTA AAA GTA GGA AGG GCA CTT TGT GGC ATA AT</td><td>992-961</td><td>250</td></tr><tr><td rowspan="3">csxB/xcbB</td><td>F173</td><td>TGT CCC CAA CCG TTT ATT GG</td><td>429-448</td><td>900</td><td rowspan="3">Mothershed, EA, 2004</td></tr><tr><td>R237</td><td>TGC TGC TAT CAT AGC CGC C</td><td>493-475</td><td>900</td></tr><tr><td>Pb196 CY5 (BHQ)</td><td>TGT TTG CCC ACA TGA ATG GCG G</td><td>452-473</td><td>250</td></tr><tr><td rowspan="3">csb/synD</td><td>F737</td><td>GCT ACC CCA TTT CAG ATG ATT TGT</td><td>737-760</td><td>900</td><td rowspan="3">Wang, X, 2011</td></tr><tr><td>R882</td><td>ACC AGC CGA GGG TTT ATT TCT AC</td><td>905-883</td><td>900</td></tr><tr><td>Pb839i CY5 (BHQ)</td><td>AAG AGA TGG GYA ACA ACT ATG TAA TGT CTT TAT TT</td><td>839-873</td><td>250</td></tr><tr><td rowspan="3">csc/synE</td><td>F478</td><td>CCC TGA GTA TGC GAA AAA AAT T</td><td>591-612</td><td>600</td><td rowspan="3">Mothershed, EA, 2004</td></tr><tr><td>R551</td><td>TGC TAA TCC CGC CTG AAT G</td><td>664-646</td><td>600</td></tr><tr><td>Pb495i FAM (BHQ)</td><td>TTT CAA TGC TAA TGA ATA CCA CCG TTT TTT TGC</td><td>612-644</td><td>200</td></tr><tr><td rowspan="3">csy/synF</td><td>F787</td><td>TCC GAG CAG GAA ATT TAT GAG AAT AC</td><td>865-890</td><td>900</td><td rowspan="3">Wang, X, 2011</td></tr><tr><td>R929</td><td>TTG CTA AAA TCA TTC GCT CCA TAT</td><td>1010-987</td><td>900</td></tr><tr><td>Pb1099i HEX (BHQ)</td><td>TAT GGT GTA CGA TAT CCC TAT CCT TGC CTA TAA T</td><td>948-981</td><td>250</td></tr><tr><td rowspan="3">cseE/cap29eF</td><td>cseE-2_Fwd1</td><td>GAGGCTGGCAATGACCAATT</td><td>470-488</td><td>600</td><td rowspan="3">Present study</td></tr><tr><td>cseE-2_Rev1</td><td>CCCAGCATATCGACAACCAA</td><td>548-529</td><td>600</td></tr><tr><td>cseE-2_Pb1 FAM (MGB)</td><td>ATCTTATGTGAACGTGGCGC</td><td>490-509</td><td>200</td></tr><tr><td rowspan="3">cszC/capZD</td><td>cszC-1_Fwd1</td><td>CAG GCC GAA GAG CGT TAT CA</td><td>895-914</td><td>900</td><td rowspan="3">Present study</td></tr><tr><td>cszC-1_Rev1</td><td>CGC CAT TCA GGG CGA TT</td><td>950-934</td><td>900</td></tr><tr><td>cszC-1_Pb1 NED (MGB)</td><td>ACA GCT CTG GCC TTA G</td><td>916-931</td><td>250</td></tr><tr><td rowspan="3">cshC</td><td>cshC-1_Fwd1</td><td>AAG CCC GTT CCA AGA TCA TG</td><td>1491-1510</td><td>600</td><td rowspan="3">Present study</td></tr><tr><td>cshC-1_Rev1</td><td>GCG GTT TGG AGA AAT AAT ATG TGT T</td><td>1563-1539</td><td>600</td></tr><tr><td>cshC-1_Pb1 VIC (MGB)</td><td>AAT GTC AGC CGT AAC TT</td><td>1513-1529</td><td>200</td></tr></tbody></table>

Sequences, positions and concentrations of primers and probes used in the 4 multiplexes, MGB: Minor Groove Binding/ BHQ: Black Whole Quencher, are two different types of Quencher. The probes obtained from CDC had a BHQ while the ones designed for this study the MGB quencher. In bold are base pairs that showed a mismatch with at least one of the alleles present in pubMLST, but not with the majority of them.

https://doi.org/10.1371/journal.pone.0206453.t002

conditions were defined for the qPCR assay using quantitative Ct values. Then, 44 Malian isolates were analysed using all primers and probes in monoplex and multiplex assays in parallel to further evaluate the optimized PCR conditions for multiplex assays and ensure the multiplex assays were comparable to monoplex assays. Later, as more samples were pro- spectively characterized at Oxford University from other MenAfriCar countries, samples from Chad, Ethiopia, Ghana, Niger and Senegal were tested with the optimized multiplex assays (Fig 1). Sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) (see definitions in Table 3) were calculated for primers/probe triplets (forward
<figure><img src="image_4.png" /><figcaption>Fig 1. Summary of the methodology employed in the study and samples tested.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206453.g001

primer, reverse primer and probe for each qPCR) using the results obtained from Oxford University as a gold standard.

# Results

## Literature search and evaluation of primers and probes

Eleven papers describing 112 primers and 16 probes for gel-based or qPCR methods to detect and characterize meningococci were identified (S1 Table). Only papers that described the primer and probe sequences were used. Gel-based PCR assays targeting different genes were identified for *N. meningitidis* detection: 1 assay for crgA, 2 for porA, 2 for ctrA and 1 for sodC. Similarly, qPCR assays for characterization of genogroups were analyzed: 3 for A, 3 for W, 7 for X, 5 for B, 4 for C, 3 for Y, 1 for H, 2 for E, and 1 for Z.

<table><caption>Table 3. Definition of sensitivity, specificity, positive predictive value and negative predictive value.</caption><thead><tr><th colspan="3">Carriage of meningococci (confirmed by the gold standard in Oxford):</th></tr></thead><tbody><tr><td>PCR outcome/carriage</td><td>Carriage present</td><td>No carriage</td></tr><tr><td>Positive</td><td>True Positive (TP)</td><td>False Positive (FP)</td></tr><tr><td>Negative</td><td>False Negative (FN)</td><td>True Negative (TN)</td></tr></tbody></table>

Sensitivity (SS): Proportion of samples positive by the gold standard technique (true carriers) that were positive by PCR. $Sensitivity = TP/(TP+FN)$

Specificity (SP): Proportion of samples negative by the gold standard technique (true non-carriers) that were negative by PCR. $Specificity = TN/(TN+FP)$

Positive Predictive Value (PPV): Probability to have a positive result with the new PCR when it is positive with the gold standard (carriage of meningococci or carriage of the tested strain). $PPV = TP/(TP+FP)$

Negative Predictive Value (NPV): Probability to have a negative result with the new PCR when it is negative with the gold standard (no meningococcal carriage or carriage of the tested strain). $NPV = TN/(TN+FN)$

https://doi.org/10.1371/journal.pone.0206453.t003
For the gel-based PCR assays: the same primers were often chosen for detection and genogrouping in different papers. The most frequent were for characterization of genogroups A, W, Y, B, C [22], and X [33]. For the qPCR assays: primers/probe triplets (forward primer, reverse primer and probe) were retrieved for *porA* (1 triplet), *ctrA* (2 triplets), and *sodC* (1 triplet) for the detection of *N. meningitidis* and for genogroups A (1 triplet), W (2 triplets), X (1 triplet), B (3 triplets), C (3 triplets), Y (2 triplets). In one of the most recent comprehensive assays, two multiplex PCR for genogroups A, W, X and B, C, Y were used successively [23].

Nineteen (14 primers and 5 probes-see Tables 1 and 2) out of the 68 primers and probes sequences evaluated within the PubMLST database had 1 to 3 mismatches, compared with the consensus sequences, which was considered acceptable. The Tm and dimer analyses generally gave acceptable results for gel-based and qPCR, except for H and Z that had to be redesigned to fit the cycling conditions in the gel-based PCR assay (Tables A, B, C, D, E, F, G and H in S2 Table). A final list of 47 primers and 12 probes was evaluated for our detection and characterization of the PCR algorithm (Tables 1 and 2). Primers and probes regions are represented schematically on their respective loci in Panels A, B, C and D in S1 Fig.

**Optimisation of the four multiplex assays (1- porA, sodC, cnl; 2- A, W, X; 3- Y, B, C; 4- H, E, Z)**

Sixteen positive controls from different genogroups were used to optimize both the gel-based and the real-time PCR assays.

For gel-based assay, the optimal MgCl2 concentration for the assay was determined to be 10.4mM and this concentration was retained for the remaining experiments. Each reaction consisted of a 30µl solution with 0.5µl of each primer, 4µl of rIMMx (monoplex) or 6µl of rIMMx (multiplex) and molecular biology grade water (ddH2O).

Table 2 presents the optimal concentrations for the selected primers and probes for real-time PCR (qPCR). The qPCR reaction consisted of a 25µl solution with 12.5 µl of 2X Gene Expression Master Mix and ddH₂O.

**Evaluation of sensitivity, specificity, PPV and NPV compared with our gold standard (S3 Table).**

**Gel-based assay (Table A in S3 Table):**

Sixteen positive controls and the samples from Mali (n = 44) were used to compare monoplex and multiplex PCR. All meningococci identified at the University of Oxford were correctly detected and characterized by monoplex gel-based PCR resulting in a sensitivity of 100% for all genes except 75% for genogroup Y (see S3 Table). The *porA* and *sodC* sensitivities were decreased to 82.14% and 96.43% respectively when working in multiplex whereas it was 100% for *cnl* and all genogroups. The specificity of *porA* assay increased from 77.78% in monoplex to 88.89% in multiplex, whereas specificity of *sodC* decreased from 29.64% in monoplex to 25.93% in multiplex. Specificity was 100% for *cnl* and all genogroups except X and Z with 96.55% and 96.30% respectively in monoplex and 100% for *cnl* and all genogroups except X and B with 96.55% and 96.15% respectively in multiplex (Supplemental Table 3).

**Real-time PCR assay (Table B in S3 Table):**

Starting with 100ng of DNA, eight 2-fold serial dilutions were performed to determine the Ct limit for categorizing a sample as positive. The last dilution, at a concentration of 0.79ng, systematically had a Ct lower or equal to 30. This is roughly equivalent to 40,000 genome copies which is far lower than the expected number of colony forming units after culture. The Ct limit was determined to be 30 cycles for crude DNA isolates. For this assay, a sample with a Ct lower or equal to 30 was considered positive and samples above 30 were considered negative (Fig 2).
<figure><img src="image_5.png" /><figcaption>**Fig 2. qPCR standard curves.** Eight two-fold dilutions of DNA from each positive control were performed and tested in a qPCR assay with appropriate primers and probes. The Ct observed for each dilution is plotted against the DNA concentration demonstrating a decrease in Ct with increased DNA concentration and detection for all genogroups with a threshold of 30. *sodC* was empirically tested later and demonstrated the same threshold (data not shown).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0206453.g002

The sensitivity was 100% for all PCR assays evaluated in this study both in monoplex and multiplex. Specificity was poor for *sodC* with 25.93% in monoplex and 29.63% in multiplex. Specificity for *porA* was slightly decreased from 81.48% to 74.07% in multiplex. All *cnl* and genogroups had a specificity of 100% both in monoplex and multiplex except genogroup X with 96.55% in monoplex. No statistical differences were observed between the new MGB *sodC* probe and the original BHQ *sodC* probe; however, the new MGB probe detected *N. meningitidis* at lower Ct values (data not shown) and was included in the final assay.

**Evaluation of the optimized multiplex assay using the pool of African carriage isolates (Table 4)**

The best performing assay was tested in multiplex on 203 isolates from five other countries of the MenAfriCar study with the gel-based and qPCR assays making a total of 247 samples. Again, the comparison was made with the gold standard results from Oxford: for the detection by gel-based PCR assay, the best sensitivity was obtained with *sodC* (80.92%—with a specificity of 69.47%) and the best specificity with *porA* (93.68% with a sensitivity of 70.39%). Sensitivity for *cnl* and all genogroups ranged from 66.67% (genogroup X) to 100% for genogroup B and E, whereas specificity ranged from 97.99% (*genogroup X*) and 100% (genogroups C, Y, H and Z).

For the qPCR assay, the best sensitivity and specificity was obtained with *porA* (96.05% and 91.58% respectively) showing no advantage of using *sodC* (94.74% and 77.89% respectively). For *cnl* and all genogroups, sensitivity ranged between 90% (genogroup A and Y) and 100%, whereas specificity ranged between 99.33% (genogroup X) and 100% (Table 4).

The benefit of evaluating results from *sodC* and *porA* together, instead of individually from the same multiplex, was analysed by comparing each primers/probe sets taken individually or together. Sensitivity and NPV were improved by using both genes together and considering
<table><caption>Table 4. Statistical analysis results for the multiplex gel-based and qPCR assays.</caption><thead><tr><th colspan="9">Gel-based multiplex</th></tr><tr><th></th><th>SS</th><th>SS 95% CI</th><th>SP</th><th>SP 95% CI</th><th>PPV</th><th>PPV 95% CI</th><th>NPV</th><th>NPV 95% CI</th></tr></thead><tbody><tr><td>porA</td><td>70.4</td><td>62.4, 77.5</td><td>93.7</td><td>86.8, 97.6</td><td>94.7</td><td>88.8, 98.0</td><td>66.4</td><td>57.8, 74.3</td></tr><tr><td>sodC</td><td>80.9</td><td>73.7, 86.8</td><td>69.5</td><td>59.2, 78.5</td><td>80.9</td><td>73.8, 86.8</td><td>69.5</td><td>59.2, 78.5</td></tr><tr><td>Cnl</td><td>85.5</td><td>75.6, 92.5</td><td>98.7</td><td>92.9, 100</td><td>98.5</td><td>91.8, 100</td><td>87.2</td><td>78.3, 93.4</td></tr><tr><td>A</td><td>75.0</td><td>50.9, 91.3</td><td>99.2</td><td>95.9, 100</td><td>93.8</td><td>69.8, 99.8</td><td>96.3</td><td>91.6, 98.8</td></tr><tr><td>W</td><td>81.5</td><td>61.9, 93.7</td><td>98.4</td><td>94.3, 99.8</td><td>91.7</td><td>73.0, 99.0</td><td>96.1</td><td>91.1, 98.7</td></tr><tr><td>X</td><td>66.7</td><td>9.4, 99.1</td><td>98.0</td><td>94.2, 99.6</td><td>40.0</td><td>5.3, 85.3</td><td>99.3</td><td>96.3, 100</td></tr><tr><td>B</td><td>100.0</td><td>29.2, 100</td><td>99.3</td><td>96.3, 100</td><td>75.0</td><td>19.4, 99.4</td><td>100.0</td><td>97.5, 100</td></tr><tr><td>C</td><td>80.0</td><td>28.3, 99.5</td><td>100.0</td><td>97.5, 100</td><td>100.0</td><td>39.8, 100</td><td>99.3</td><td>96.3, 100</td></tr><tr><td>Y</td><td>80.0</td><td>44.4, 97.5</td><td>100.0</td><td>97.4, 100</td><td>100.0</td><td>63.1, 100</td><td>98.6</td><td>95.1, 99.8</td></tr><tr><td>H</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>E</td><td>100.0</td><td>2.5, 100</td><td>99.3</td><td>96.3, 100</td><td>50.0</td><td>1.3, 98.7</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>Z</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>porA/sodC</td><td>88.8</td><td>82.7, 93.3</td><td>67.4</td><td>57.0, 76.6</td><td>81.3</td><td>74.6, 86.9</td><td>79.0</td><td>68.5, 87.3</td></tr><tr><th colspan="9">Real-time multiplex</th></tr><tr><th></th><th>SS</th><th>SS 95% CI</th><th>SP</th><th>SP 95% CI</th><th>PPV</th><th>PPV 95% CI</th><th>NPV</th><th>NPV 95% CI</th></tr><tr><td>porA</td><td>96.1</td><td>91.6, 98.5</td><td>91.6</td><td>84.1, 96.3</td><td>94.8</td><td>90.0, 97.7</td><td>93.5</td><td>86.5, 97.6</td></tr><tr><td>sodC</td><td>94.7</td><td>89.9, 97.7</td><td>77.9</td><td>68.2, 85.8</td><td>87.3</td><td>81.2, 91.9</td><td>90.2</td><td>81.7, 95.7</td></tr><tr><td>Cnl</td><td>94.7</td><td>87.1, 98.5</td><td>100.0</td><td>95.3, 100</td><td>100.0</td><td>95.0, 100</td><td>95.0</td><td>87.7, 98.6</td></tr><tr><td>A</td><td>90.0</td><td>68.3, 98.8</td><td>100.0</td><td>97.2, 100</td><td>100.0</td><td>81.5, 100</td><td>98.5</td><td>94.7, 99.8</td></tr><tr><td>W</td><td>92.6</td><td>75.7, 99.1</td><td>100.0</td><td>97.1, 100</td><td>100.0</td><td>86.3, 100</td><td>98.4</td><td>94.4, 99.8</td></tr><tr><td>X</td><td>100.0</td><td>29.2, 100</td><td>99.3</td><td>96.3, 100</td><td>75.0</td><td>19.4, 99.4</td><td>100.0</td><td>97.5, 100</td></tr><tr><td>B</td><td>100.0</td><td>29.2, 100</td><td>100.0</td><td>97.5, 100</td><td>100.0</td><td>29.2, 100</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>C</td><td>100.0</td><td>47.8, 100</td><td>100.0</td><td>97.5, 100</td><td>100.0</td><td>47.8, 100</td><td>100.0</td><td>97.5, 100</td></tr><tr><td>Y</td><td>90.0</td><td>55.5, 99.7</td><td>100.0</td><td>97.4, 100</td><td>100.0</td><td>66.4, 100</td><td>99.3</td><td>96.2, 100</td></tr><tr><td>H</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>E</td><td>100.0</td><td>2.5, 100</td><td>99.3</td><td>96.4, 100</td><td>50.0</td><td>1.3, 98.7</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>Z</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td><td>NA</td><td>NA</td><td>100.0</td><td>97.6, 100</td></tr><tr><td>porA/sodC</td><td>96.1</td><td>91.6, 98.5</td><td>77.9</td><td>68.2, 85.8</td><td>87.4</td><td>81.4, 92.0</td><td>92.5</td><td>84.4, 97.2</td></tr></tbody></table>

Results of both PCR assays tested only in multiplex on a panel of 247 samples from Chad, Ethiopia, Ghana, Mali, Niger and Senegal. Sensitivity, Specificity, PPV and NPV were assessed for each primers pair and primers/probe in the multiplex assays and presented as percentages (%) with 95% confidence intervals (CI).

https://doi.org/10.1371/journal.pone.0206453.t004

positive any samples that were positive for at least one of the targets for the gel-based PCR assay in comparison with *porA* or *sodC* alone from the same multiplex whereas the specificity and PPV were decreased. For the qPCR assay all statistical results for *sodC* were lower or equal to those for *porA*. Inclusion of *sodC* did not provide additional value (Table 4).

From the panel of 247 samples from all six different countries, 152 were characterized as *N. meningitidis* at the University of Oxford using *rplF* sequencing [29]. Of these, 88.82% (135) isolates were correctly identified as meningococci by the gel-based PCR assay and 96.05% (146) by the qPCR assay. All genogroups, but H and Z, were identified in the panel tested: the most commonly represented was the *cnl* with 43.9%/48.4% (gel based/qPCR). Some genogroups considered as rare in the region were detected, such as B (2%) and E (0.66%) by both assays.
Discussion

Conventional bacteriological methods that are usually routinely used to detect and character-
ize meningococci are complex, long and involve several overnight cultures and chemical tests
which might generate error [15]. Developing a simple and rapid PCR algorithm that detects
and characterizes carriage of meningococci is important, especially in the meningitis belt
where recent vaccinations with conjugate vaccines were carried out and carriage is usually low
[14]. Therefore, it is important to correctly detect *N. meningitidis* in a context where transmis-
sion is high especially in young children and where the period of carriage is short [34].

In this study, after an *in silico* evaluation of the existing primers and probes, we evaluated a multiplex PCR algorithm by comparison with techniques that were validated and published in a reference laboratory at the University of Oxford (24, 30). Our multiplex PCR algorithm was later used for a meningococcal carriage study in The Gambia and improved the detection of carried meningococci by more than twice (239%) [35]. To our knowledge, this is the first time that such an extensive evaluation was made for a carriage study in the meningitis belt, using the most advanced tools and following recommended advanced guidelines for qualitative techniques [32].

To optimise the sensitivity of our algorithm for the detection of meningococci, we targeted both *porA* and *sodC* genes simultaneously. The published *sodC* gel-based PCR primers were re-designed due to too many mismatches against the consensus sequence [18]. Similarly, primers and probe of qPCR targeting *porA* available from the literature [17] yielded too many mismatches and were re-designed entirely. Designing qPCR primers and probes for *cnl* identification was not an easy task as the *cnl* sequence is short and variable, but the design of two complementary forward primers (Table 2), allowed the complete coverage of the *cnl* sequences that were present in pubMLST at that time. Recent comparison of primers/probes sequences used in this study with a larger number of sequences now available in the pubMLST identified additional point mutations in newly published alleles (data not shown). Although most alleles are still covered by the primers and probes developed here, new variations within these primers/probes sequences may lead to decreased sensitivity and we advise a complete new *in silico* evaluation of primers and probes before using these. Since pubMLST is exponentially growing, we recommend our methodology to be applied—or, at a minimum, blasting of primers and probes to make them more robust—when developing or implementing new molecular diagnosis tools for *N. meningitidis* as well as other pathogens. The variability of the targeted genes between meningococcal strains might explain the sensitivity and specificity variations observed between monoplex and multiplex assays and future studies will be needed to investigate this further for the detection of *N. meningitidis*.

A major strength of the developed algorithm is that the first multiplex assay (*porA*, *sodC* and *cnl*) allows the simultaneous detection of meningococci and the presence/absence of their capsule locus. Since *cnl* genotype is most common among meningococcal strains in carriage [11, 14], this assay improves the rapidity of *N. meningitidis* characterization. The sequential use of the following three multiplex assays for genogroup characterization, from the most frequent groups (A, W, and X followed by B, C, and Y) to the carriage uncommon genogroups (E, H and Z) in the meningitis belt, further reduces the time required for meningococcal characterization. The primers/probe sequences used for these genogroup assays could also be modified to fit the variable epidemiology of meningococci by region and over time.

In this study, we have developed an algorithm that: (i) has a high sensitivity for the detec-
tion of meningococci; and, (ii) allows rapid characterization of all genogroups with high speci-
ficity. A strength of gel-based assays is that the assays joined two commonly used genes to
detect *N. meningitidis*: *porA* and *sodC* in multiplex, which improved the sensitivity from
70.39% for *porA* and 80.92% for *sodC* alone to 88.82% (see Table 4) for meningococcal detection. However, no improvement in sensitivity and specificity was observed when adding *sodC* in the multiplex qPCR assay and we advise using only *porA/cnl* for the qPCR. This might be due to the modification of the published *sodC* probe to accommodate a MGB quencher in order to harmonize the quenchers in the first multiplex (*porA*, *sodC*, *cnl*). Indeed, “BHQ” probes are difficult to obtain from oligonucleotide suppliers, because of their unusual design [23]. A subsequent comparison of 44 Malian samples tested with both BHQ and MGB probes showed no differences in detection of isolates. However, the efficacy of the shortened MGB *sodC* probe should be further investigated.

For genogroup characterization, multiplex qPCR that has already been evaluated were validated in gel-based PCR platform. The statistical results varied from one genogroup to the other (Table 4). For the two multiplexes (A, W, and X; B, C, and Y), published primers/probes were used and their statistical results were acceptable for all genogroups. The genogroup X primers set had the lowest sensitivity and PPV in the gel-based assay and may require redesign. Overall, the qPCR assay yielded better results than the gel-based PCR and should be the preferred assay for genogroup characterization.

Our study used the laboratory results obtained from Oxford University as reference standard for the statistical analysis. The sensitivity, specificity, PPV and NPV may be affected by the transport conditions when DNA samples were shipped from MenAfriCar countries and how the samples were processed at Oxford University. This may have affected our comparisons between PCR tests realized in Bamako and analyses in Oxford. In addition, no qPCR assay for the rare genogroups H, E, and Z was available in Oxford, making the direct comparison difficult. The results generated in this study were compared to the *ctrA* sequencing results when available and may explain some of the observed discrepancies; indeed, the comparison of two different techniques could explain the low PPV for group E in both assays.

The quality of the template may have also affected the detection of meningococci as most studies use DNA purified by kits. Boiled cell suspensions were purposely used to make our assays as inexpensive as possible and the good statistics obtained with those crude samples demonstrate the efficacy of the assays. Stringent rules were applied to determine whether a test was positive or not: by qPCR, a CT threshold of 30 was chosen based on preliminary results obtained (Fig 2), whereas the CDC/WHO reference guidelines recommend a threshold of 35 which, with the reagents and conditions used in our study, would have led to a considerable number of false positives.

In conclusion, this algorithm allows the detection of *N. meningitidis* with high sensitivity for samples from 6 different countries of the meningitis belt, it significantly reduces the number of tests required for meningococcal detection and allows very specific characterization of samples for carriage studies in comparison with conventional bacteriological techniques.

**Supporting information**

**S1 Fig. Primers and probes positioned on appropriate genes.** Gel-based (green arrows) and rt (black arrows) PCR primers and probes for each multiplex (A-multiplex 1, B-multiplex 2, C-multiplex 3 and D-multiplex 4); the gene arrow pointing toward the 3' end of the gene. Both forward primers (available in the literature and newly designed) for *porA* are shown on the gene (A). Figure adapted from Harisson OB and al. 2013. (PDF)

**S1 Table. List of papers that were selected for further analyses.**
(XLSX)
**S2 Table. Primer and probe sequences general technical analysis.** Results of the BLAST analyses performed for each primer and probe using the NCBI BLAST tool. All BLAST results shown were done in the nucleotide collection database, against *N. meningitidis* (taxon id 487). The Tm and possible dimers were obtained using the primer express design tool.
(XLSX)

**S3 Table.** A. Statistical results for gel-based PCR assay. Statistical analysis results of the gel-based PCR assay. Sensitivity, Specificity, Positive Predictive Value (PPV) and Negative Predictive Value (NPV) were assessed for each primers pair in the multiplex and monoplex assay pooling results from test on positive controls (n = 16) and the Malian samples (n = 44) and presented as percentages (%). B. Statistical results for the qPCR assay. Statistical analysis results of the qPCR assay. Sensitivity, Specificity, Positive Predictive Value (PPV) and Negative Predictive Value (NPV) were assessed for each primers pair in the multiplex and monoplex assay pooling results from test on positive controls (n = 16) and the Malian samples (n = 44) and presented as percentages (%).
(DOCX)

**S1 Dataset. Database used for the analyses for the Chadian Site.**
(XLSX)

**S2 Dataset. Database used for the analyses for the Ethiopian Site.**
(XLSX)

**S3 Dataset. Database used for the analyses for the Ghanaian Site.**
(XLSX)

**S4 Dataset. Database used for the analyses for the Malian Site.**
(XLSX)

**S5 Dataset. Database used for the analyses for the Niger Site.**
(XLSX)

**S6 Dataset. Database used for the analyses for the Senegalese Site.**
(XLSX)

**Acknowledgments**

K.D. Is a Wellcome Trust Tropical Fellow. O.B.H. and M.C.J.M. were supported by the Wellcome Trust (grant 087622). This research made use of the *Neisseria* Sequence Typing database website (http://pubmlst.org/neisseria/) developed by Keith A. Jolley and Martin C. J. Maiden, hosted by the University of Oxford, and supported by the Wellcome Trust (grant 104992).

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.

^Institutions and individual members of the MenAfriCar consortium

Armauer Hansen Research Institute, Addis Ababa, Ethiopia: Oumer Ali, Abraham Aseffa (PI), Ahmed Bedru, Tsehaynesh Lema, Tesfaye Moti, Yenenesh Tekletsion, Alemayehu Worku, Haimanot Guebre Xabher (deceased), Lawrence Yamuah.

Centre de Recherche Médicale et Sanitaire (CERMES), Niamey, Niger (Member of the International Network of Pasteur Institutes): Rahamatou Moustapha Boukary, Jean-Marc Collard (PI), Ibrahim Dan Dano, Ibrahim Habiboulaye, Bassira Issaka, Jean-François Jusot, Sani Ousmane, Issoufa Rabe.

Centre de Support en Santé International (CSSI), N’Djamena, Chad: Doumagoum Moto Daugla (PI), Jean Pierre Gami, Kadidja Gamougam, Lodoum Mbainadji, Nathan Naibei, Maxime Narbé, Jacques Toralta.
**Centre pour les Vaccins en Développement, Bamako, Mali:** Abdoulaye Berthe, Kanny Diallo, Mahamadou Keita, Uma Onwuchekwa, Samba O Sow (PI), Boubou Tamboura, Awa Traore, Aliou Toure.

**Centers for Disease Control, Atlanta, USA:** Tom Clark, Xin Wang, Leonard Mayer.

**Department of Community Medicine, University of Maiduguri, Maiduguri, Nigeria:** Mary Amodu, Omeiza Beida, Galadima Gadzama, Babatunji Omotara (PI), Zailani Sambo, Shuaibu Yahya.

**Faculty of Infectious Disease, London School of Hygiene & Tropical Medicine, London, UK:** Daniel Chandramohan, Brian M Greenwood (PI), Musa Hassan-King, Olivier Manigart, Maria Nascimento, James M Stuart, Arouna Woukeu.

**Princeton University, USA:** Nicole E Basta

**Public Health England Vaccine Evaluation Unit, Manchester, UK:** Xilian Bai, Jay Lucidarme, Ray Borrow, Helen Findlow.

**Institut de Recherche pour le Développement, Dakar, Senegal:** Serge Alavo, Hubert Bassene, Aldiouma Diallo (PI), Marietou Dieng, Souleymane Doucouré, Jules François Gomis, Assane Ndiaye, Cheikh Sokhna, Jean François Trape.

**Navrongo Health Research Centre, Navrongo, Ghana:** Bugri Akalifa (deceased), Abudulai Forgor, Abraham Hodgson (PI), Isaac Osei, Stephen L Quaye, John Williams, Peter Wontuo.

**University of Bristol, UK:** Thomas Irving.

**University of Cambridge, UK:** Caroline L Trotter

**University of Oxford, UK:** Julia Bennett, Dorothea Hill, Odile Harrison, Martin C Maiden, Lisa Rebbetts, Eleanor Watkins.

MenAfriCar PI and corresponding author: Professor Brian Greenwood brian.green-wood@lshtm.ac.uk

## Author Contributions

**Conceptualization:** Samba O. Sow, Ray Borrow, Brian M. Greenwood, Martin C. J. Maiden, Olivier Manigart.

**Formal analysis:** Kanny Diallo, Mamadou D. Coulibaly, Lisa S. Rebbetts, Olivier Manigart.

**Methodology:** Ray Borrow.

**Resources:** Kadidja Gamougam, Yenenesh K. Tekletsion, Akalifa Bugri, Aliou Toure, Bassira Issaka, Marietou Dieng, Jean-Marc Collard.

**Supervision:** Odile B. Harrison, Jay Lucidarme, Caroline Trotter, Jean-Marc Collard, Samba O. Sow, Xin Wang, Leonard W. Mayer, Ray Borrow, Brian M. Greenwood, Martin C. J. Maiden.

**Writing – original draft:** Caroline Trotter, Brian M. Greenwood, Olivier Manigart.

**Writing – review & editing:** Kanny Diallo, Lisa S. Rebbetts, Odile B. Harrison, Jay Lucidarme, Caroline Trotter, Jean-Marc Collard, Samba O. Sow, Xin Wang, Leonard W. Mayer, Ray Borrow, Brian M. Greenwood, Martin C. J. Maiden, Olivier Manigart.

## References

1. Lapeyssonnie L. [Comparative epidemiologic study of meningococcic cerebrospinal meningitis in temperate regions and in the meningitis belt in Africa. Attempt at synthesis]. *Medecine tropicale*: revue du Corps de sante colonial. 1968; 28(6):709-20. PMID: 5739513.

2. Broome CV. The carrier state: Neisseria meningitidis. *The Journal of antimicrobial chemotherapy*. 1986; 18 Suppl A:25-34. PMID: 3091563.
3. Stephens DS, Greenwood B, Brandtzaeg P. Epidemic meningitis, meningococcaemia, and Neisseria meningitidis. Lancet. 2007; 369(9580):2196-210. https://doi.org/10.1016/S0140-6736(07)61016-2 PMID: 17604802.

4. Harrison OB, Brueggemann AB, Caugant DA, van der Ende A, Frosch M, Gray S, et al. Molecular typing methods for outbreak detection and surveillance of invasive disease caused by Neisseria meningitidis, Haemophilus influenzae and Streptococcus pneumoniae, a review. Microbiology. 2011; 157(Pt 8):2181-95. https://doi.org/10.1099/mic.0.050518-0 PMID: 21622526.

5. Nicolas P. [Meningococcal meningitis epidemics in sub-Saharan Africa and the meningococcal A conjugate vaccine]. Medecine et sante tropicales. 22(3):246-58. https://doi.org/10.1684/mst.2012.0086 PMID: 23186941.

6. Gamougam K, Daugla DM, Toralta J, Ngadoua C, Fermon F, Page AL, et al. Continuing effectiveness of serogroup A meningococcal conjugate vaccine, Chad, 2013. Emerging infectious diseases. 2015; 21 (1):115-8. https://doi.org/10.3201/eid2101.140256 PMID: 25536336; PubMed Central PMCID: PMC4285275.

7. Kristiansen PA, Jorgensen HJ, Caugant DA. Serogroup A meningococcal conjugate vaccines in Africa. Expert review of vaccines. 2015:1-18. https://doi.org/10.1586/14760584.2015.1084232 PMID: 26358167.

8. Maurice J. Vaccine shortage threatens spread of meningitis in Niger. Lancet. 2015; 385(9984):2241. https://doi.org/10.1016/S0140-6736(15)61050-9 PMID: 26088485.

9. Abio A, Neal KR, Beck CR. An epidemiological review of changes in meningococcal biology during the last 100 years. Pathogens and global health. 2013; 107(7):373-80. https://doi.org/10.1179/2047773213Y.0000000119 PMID: 24392681; PubMed Central PMCID: PMC4083158.

10. Caugant DA, Tzanakaki G, Kriz P. Lessons from meningococcal carriage studies. FEMS microbiology reviews. 2007; 31(1):52-63. https://doi.org/10.1111/j.1574-6976.2006.00052.x PMID: 17233635.

11. Claus H, Maiden MC, Maag R, Frosch M, Vogel U. Many carried meningococci lack the genes required for capsule synthesis and transport. Microbiology. 2002; 148(Pt 6):1813-9. https://doi.org/10.1099/00221287-148-6-1813 PMID: 12055301.

12. Kristiansen PA, Diomande F, Wei SC, Ouedraogo R, Sangare L, Sanou I, et al. Baseline meningococcal carriage in Burkina Faso before the introduction of a meningococcal serogroup A conjugate vaccine. Clinical and vaccine immunology: CVI. 2011; 18(3):435-43. https://doi.org/10.1128/CVI.00479-10 PMID: 21228139; PubMed Central PMCID: PMC3067389.

13. Kristiansen PA, Ba AK, Ouedraogo AS, Sanou I, Ouedraogo R, Sangare L, et al. Persistent low carriage of serogroup A Neisseria meningitidis two years after mass vaccination with the meningococcal conjugate vaccine, MenAfriVac. BMC Infect Dis. 2014; 14:663. https://doi.org/10.1186/s12879-014-0663-4 PMID: 25472422; PubMed Central PMCID: PMC4267149.

14. MenAfriCar. The diversity of meningococcal carriage across the African meningitis belt and the impact of vaccination with a group A meningococcal conjugate vaccine. The Journal of infectious diseases. 2015. https://doi.org/10.1093/infdis/jiv211 PMID: 25858956.

15. MenAfriCar. Meningococcal carriage in the African meningitis belt. Tropical medicine & international health: TM & IH. 2013; 18(8):968-78. https://doi.org/10.1111/tmi.12125 PMID: 23682910; PubMed Central PMCID: PMC3950817.

16. Guiver M, Borrow R, Marsh J, Gray SJ, Kaczmarski EB, Howells D, et al. Evaluation of the Applied Bio-systems automated Taqman polymerase chain reaction system for the detection of meningococcal DNA. FEMS immunology and medical microbiology. 2000; 28(2):173-9. PMID: 10799809.

17. Jordens JZ, Heckels JE. A novel porA-based real-time PCR for detection of meningococcal carriage. Journal of medical microbiology. 2005; 54(Pt 5):463-6. https://doi.org/10.1099/jmm.0.45847-0 PMID: 15824424.

18. Dolan Thomas J, Hatcher CP, Satterfield DA, Theodore MJ, Bach MC, Linscott KB, et al. sodC-based real-time PCR for detection of Neisseria meningitidis. PloS one. 2011; 6(5):e19361. Epub 2011/05/17. https://doi.org/10.1371/journal.pone.0019361 PMID: 21573213; PubMed Central PMCID: PMC3088665.

19. Jordens JZ, Williams JN, Jones GR, Heckels JE. Detection of meningococcal carriage by culture and PCR of throat swabs and mouth gargles. Journal of clinical microbiology. 2002; 40(1):75-9. https://doi.org/10.1128/JCM.40.1.75-79.2002 PMID: 11773095; PubMed Central PMCID: PMC120086.

20. Bennett DE, Cafferkey MT. Consecutive use of two multiplex PCR-based assays for simultaneous identification and determination of capsular status of nine common Neisseria meningitidis serogroups associated with invasive disease. Journal of clinical microbiology. 2006; 44(3):1127-31. https://doi.org/10.1128/JCM.44.3.1127-1131.2006 PMID: 16517911; PubMed Central PMCID: PMCPMC1393079.

21. Mothershed EA, Sacchi CT, Whitney AM, Barnett GA, Ajello GW, Schmink S, et al. Use of real-time PCR to resolve slide agglutination discrepancies in serogroup identification of Neisseria meningitidis.
Journal of clinical microbiology. 2004; 42(1):320-8. Epub 2004/01/13. https://doi.org/10.1128/JCM.42.1.320-328.2004 PMID: 14715772; PubMed Central PMCID: PMC321732.

22. Taha MK. Simultaneous approach for nonculture PCR-based identification and serogroup prediction of Neisseria meningitidis. Journal of clinical microbiology. 2000; 38(2):855-7. PMID: 10655397.

23. Wang X, Theodore MJ, Mair R, Trujillo-Lopez E, du Plessis M, Wolter N, et al. Clinical validation of multiplex real-time PCR assays for detection of bacterial meningitis pathogens. Journal of clinical microbiology. 2012; 50(3):702-8. https://doi.org/10.1128/JCM.06087-11 PMID: 22170919; PubMed Central PMCID: PMC3295090.

24. Zhu H, Wang Q, Wen L, Xu J, Shao Z, Chen M, et al. Development of a multiplex PCR assay for detection and genogrouping of Neisseria meningitidis. Journal of clinical microbiology. 50(1):46-51. https://doi.org/10.1128/JCM.00918-11 PMID: 22090406.

25. LaForce FM, Konde K, Viviani S, Preziosi MP. The Meningitis Vaccine Project. Vaccine. 2007; 25 Suppl 1:A97-100. https://doi.org/10.1016/j.vaccine.2007.04.049 PMID: 17521780.

26. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990; 215(3):403-10. https://doi.org/10.1016/S0022-2836(05)80360-2 PMID: 2231712.

27. Jolley K. PubMLST, Neisseria sequence Typing 2010. Available from: http://pubmlst.org/neisseria/.

28. Popovic TA, G; Facklam R. Laboratory manual for the diagnosis of meningitis caused by *Neisseria meningitidis*, *Streptococcus pneumoniae* and *Haemophilus influenzae*. Geneva: World Health Organization; 1999.

29. Bennett JS, Watkins ER, Jolley KA, Harrison OB, Maiden MC. Identifying Neisseria species by use of the 50S ribosomal protein L6 (rplF) gene. Journal of clinical microbiology. 2014; 52(5):1375-81. Epub 2014/02/14. https://doi.org/10.1128/JCM.03529-13 PMID: 24523465; PubMed Central PMCID: PMC3993661.

30. tecnologies ID. PrimerQuest. Available from: https://eu.idtdna.com/PrimerQuest/Home/Index.

31. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods in molecular biology (Clifton, NJ. 2000; 132:365-86. PMID: 10547847.

32. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009; 55(4):611-22. Epub 2009/02/28. https://doi.org/10.1373/clinchem.2008.112797 PMID: 19246619.

33. Yaro S, Drabo A, Ouangraoua S, Kirakoya-Samadoulougou F, Mueller J, Sanou O, et al. *Neisseria meningitis* serogroup X outbreak in Burkina Faso, 2009-2010. *Open Journal of Internal Medicine* 2012. p. 41-29.

34. MenAfriCar C. Household transmission of Neisseria meningitidis in the African meningitis belt: a longitudinal cohort study. Lancet Glob Health. 2016; 4(12):e989-e95. Epub 2016/11/20. https://doi.org/10.1016/S2214-109X(16)30244-3 PMID: 27855873.

35. Manigart O, Okeakpu J, Odutola A, Jarju S, Foster-Nyarko E, Diallo K, et al. Alternative Molecular Methods for Improved Detection of Meningococcal Carriage and Measurement of Bacterial Density. Journal of clinical microbiology. 2016. https://doi.org/10.1128/JCM.01428-16 PMID: 27582517.