Skip directly to search Skip directly to A to Z list Skip directly to page options Skip directly to site content

Volume 24, Number 4—April 2018


Evolution of Sequence Type 4821 Clonal Complex Meningococcal Strains in China from Prequinolone to Quinolone Era, 1972–2013



Technical Appendices

Article Metrics



Views: Views equals page views plus PDF downloads

Metric Details
Qinglan Guo1, Mustapha M. Mustapha1, Mingliang Chen, Di Qu, Xi Zhang, Min ChenComments to Author , Yohei DoiComments to Author , Minggui WangComments to Author , and Lee H. Harrison
Author affiliations: Fudan University Huashan Hospital, Shanghai, China (Q. Guo, M. Wang); University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA (M.M. Mustapha, Y. Doi, L.H. Harrison); Shanghai Municipal Center for Disease Control and Prevention, Shanghai (M. Chen, X. Zhang, M. Chen); Fudan University Shanghai Medical College, Shanghai (D. Qu)

Cite This Article


Highlight and copy the desired format.

EID Guo Q, Mustapha MM, Chen M, Qu D, Zhang X, Chen M, et al. Evolution of Sequence Type 4821 Clonal Complex Meningococcal Strains in China from Prequinolone to Quinolone Era, 1972–2013. Emerg Infect Dis. 2018;24(4):683-690.
AMA Guo Q, Mustapha MM, Chen M, et al. Evolution of Sequence Type 4821 Clonal Complex Meningococcal Strains in China from Prequinolone to Quinolone Era, 1972–2013. Emerging Infectious Diseases. 2018;24(4):683-690. doi:10.3201/eid2404.171744.
APA Guo, Q., Mustapha, M. M., Chen, M., Qu, D., Zhang, X., Chen, M....Harrison, L. H. (2018). Evolution of Sequence Type 4821 Clonal Complex Meningococcal Strains in China from Prequinolone to Quinolone Era, 1972–2013. Emerging Infectious Diseases, 24(4), 683-690.


The expansion of hypervirulent sequence type 4821 clonal complex (CC4821) lineage Neisseria meningitidis bacteria has led to a shift in meningococcal disease epidemiology in China, from serogroup A (MenA) to MenC. Knowledge of the evolution and genetic origin of the emergent MenC strains is limited. In this study, we subjected 76 CC4821 isolates collected across China during 1972–1977 and 2005–2013 to phylogenetic analysis, traditional genotyping, or both. We show that successive recombination events within genes encoding surface antigens and acquisition of quinolone resistance mutations possibly played a role in the emergence of CC4821 as an epidemic clone in China. MenC and MenB CC4821 strains have spread across China and have been detected in several countries in different continents. Capsular switches involving serogroups B and C occurred among epidemic strains, raising concerns regarding possible increases in MenB disease, given that vaccines in use in China do not protect against MenB.

The incidence of meningococcal disease and Neisseria meningitidis strain distribution vary over time, within and between countries and regions (1). Six serogroups (A, B, C W, X, and Y) account for nearly all cases of invasive meningococcal disease (IMD) globally (1). Serogroup C (MenC) cases were rare in China until 2003–2005, when several MenC outbreaks were reported in Anhui Province (24). These outbreaks were caused by a previously unreported hypervirulent clonal complex (CC) 4821 lineage (5). A pharyngeal carriage survey and national public health surveillance during 2004–2005 identified CC4821 among sporadic IMD case-patients and asymptomatic carriers across 11 provinces, demonstrating the wide geographic distribution of CC4821 (2). During 2005–2012, MenC CC4821 became the leading cause of endemic meningococcal disease in China (6). Further analyses of historic isolate collections identified MenB and MenC CC4821 from carriage surveys in the 1970s (2,7). These studies demonstrated that CC4821 had been mostly associated with asymptomatic carriage over several decades before it emerged as a main cause of IMD (2,7,8). Also, our recent analyses of quinolone resistance among historic meningococcal isolate collections in China found a substantial temporal shift toward increased quinolone nonsusceptibility from the prequinolone era (before ≈1985) to the quinolone era (none versus >70%), particularly within hypervirulent CC4821 and CC5 lineages (7). Such findings support the hypothesis that quinolone resistance could have played a role in the emergence of MenC CC4821 outbreaks in China.

Meningococci have a dynamic genome that evolves rapidly through point mutations and frequent recombination. Such genetic changes give rise to strains with novel capsular or other major surface antigens that evade existing population immunity (9). A study by Zhu et al. found extensive genomic variation among 22 CC4821 invasive and carriage isolates from 12 provinces in China during 2005–2011 (8). In that study, CC4821 isolates belonged to 2 distinct phylogenetic groups, and results indicated that group 1, containing the epidemic reference strain 053442, might be more invasive than group 2 and that MenB and MenC coexisted within both groups 1 and 2 (8).

Our study describes phylogenetic relationships within a collection of historic and current isolates from Shanghai in the prequinolone and quinolone eras and explores how the isolates fit into the larger genetic profile of CC4821 from China (8). We aimed to shed light on the genomic factors underlying the abrupt transition of this lineage from a minority strain to a leading cause of endemic disease and outbreaks.

Materials and Methods

Strain Selection and Molecular Typing

A total of 374 meningococcal isolates were collected from IMD case-patients, close contacts, and asymptomatic carriers during pharyngeal carriage surveys and routine, laboratory-based public health surveillance conducted during 1965–1985 and 2005–2013 (7). The interruption during 1986–2004 was because of the decreased incidence of meningococcal cases and the increased ability to identify N. meningitidis in hospitals, necessitating fewer isolates to be referred (7). In all, 52 CC4821 isolates were identified and underwent molecular characterization using traditional PCR sequencing of the porA, porB, fetA, fHbp, nadA, nhba, and gyrA genes and pulsed-field gel electrophoresis as described previously (7,10,11). We selected 8 of these 52 isolates, representing strains from different periods, serogroups, or pulsed-field gel electrophoresis groups, for whole-genome sequencing (WGS) and in-depth phylogenetic analysis. We downloaded assembled contiguous genome sequences (contigs) for 24 additional genome sequences from previous studies of CC4821 in China from GenBank (NM11003, accession no. NZ_ANBU00000000) (5,8). Therefore, a total of 32 whole-genome sequences underwent core genome phylogenetic analyses for this study (Technical Appendix [PDF - 368 KB - 20 pages] Table 1).

Genome Sequencing and Assembly

We performed single molecule real-time sequencing (PacBio; Pacific Biosciences, Menlo Park, CA, USA) on 4 CC4821 isolates (NM040, NM062, NM205, and NM323) and sequenced the remaining 4 isolates (NM001, NM050, NM193, and NM313) by using Illumina HiSeq paired-end sequencing (Illumina, San Diego, CA, USA). We assembled PacBio genomes by using HGAP 4.0 ( and Illumina genomes by using SPAdes 3.7 (12). We annotated assembled contigs by using the Prokka v1.11 (13) pipeline and submitted them to the PubMLST Neisseria genome database ( with ID numbers 41414–41421, where allelic numbers were assigned to all identified genes (14).

Phylogenetic Analyses

We aligned assembled genomes (n = 32) and produced a core genome phylogenetic tree with 1,000 rapid bootstrap replicates by using RAxML and Mauve 2.3 as previously described (1517). Serogroup A reference genome Z2491 was included as an outgroup (18). We aligned gene sequences selected for focused analyses with MEGA 5.2 ( and constructed maximum-likelihood phylogenetic trees under the HKY model of evolution (19). We included comparison sequences from a global collection of 133 genomes from the PubMLST Neisseria genome database representing all major invasive disease lineages in some of these phylogenetic analyses as references.


We assessed recombination by using ClonalFrameML and Gubbins (20,21). We then mapped the recombination to the CC4821 reference genome 053442 (Technical Appendix [PDF - 368 KB - 20 pages]).

Gene Content

We assessed gene content by using the Roary 3.6 ortholog clustering program (22), which identifies presence or absence of orthologous gene sequences using a cutoff of 90% sequence identity (–i = 90). We defined core genes as genes present in >90% of the genomes. When comparing the gene content of 2 groups of genomes, we defined a gene as specific to that group if it was present in >90% of the genomes in the group and in <20% of the genomes in the comparison group. We downloaded gene functional annotation from the COG database (23) and compared genes containing recombinant sequences with nonrecombinant ones based on the major COG classes (cellular processes and signaling, information storage and processing, metabolism, poorly characterized) using uncorrected χ2 tests.


Out of 52 Shanghai isolates, 18 (34%) were from 1972–1977 (3 IMD isolates and 15 asymptomatic carriage isolates), and the remaining 34 (23 IMD isolates, 6 isolates from close contacts, and 5 carriage isolates) were isolated during 2005–2013. Most (56%) Shanghai isolates were serogroup C, 42% were serogroup B, and 1 was nongroupable (Technical Appendix [PDF - 368 KB - 20 pages] Table 2).

Phylogenetic Analyses of CC4821 Isolates in China

Figure 1

Thumbnail of Core genome phylogenetic tree showing relationships between Neisseria meningitidis serogroups C and B CC4821 strains, China, 1972–2011. The strains cluster within 1 of 2 distinct phylogenetic groups, group 1 and group 2. Within group 1 is an antigenically distinct clonal group (epidemic clone) containing outbreak-associated strains. Tree is rooted using serogroup A reference strain (Z2491) as an outgroup. Maximum-likelihood phylogenetic trees of aligned core genome sequences were ge

Figure 1. Core genome phylogenetic tree showing relationships between Neisseria meningitidis serogroups C and B CC4821 strains, China, 1972–2011. The strains cluster within 1 of 2 distinct phylogenetic groups, group 1 and group...

We conducted comparative genome analyses for 8 isolates from Shanghai and 24 publicly available genomes from across China (Figure 1). Among the 8 newly sequenced isolates, 4 were from asymptomatic carriers from 1972–1977 (NM193, NM205, NM313, and NM323), 3 were IMD isolates from 2005–2011 (NM001, NM062, and NM040), and 1 was from an asymptomatic contact of a meningitis patient (NM050). The remaining 24 genome sequences included invasive and carriage CC4821 isolates from 12 provinces across China from 2004–2011, as characterized in previous studies (4,8). Overall, MenC represented 66% (21/32) of isolates that we analyzed, and the remaining 34% (11/32) belonged to MenB.

The core genome phylogenetic tree classified CC4821 into 2 distinct groups (Figure 1). Group 1 consisted of several isolates that were very closely related; this group contained 15 isolates from 2004 to 2011, and most were from IMD case-patients. In addition, the epidemic reference strain 053442 clustered with a small group of highly similar group 1 isolates. This subgroup is defined as the epidemic clone based on core genome phylogenetic analyses and antigen gene profile (described hereafter). A second more diverse phylogenetic group (group 2) consisted of isolates from Shanghai during 1972–1977 and more recent ones collected from 9 provinces during 2005–2011. In contrast to group 1, only 29% of group 2 isolates were from IMD case-patients (Figure 1). MenB was found in both groups and was interspersed with MenC genomes. Two of 11 MenB isolates belonged to group 1, the remainder to group 2 (Figure 1; Technical Appendix [PDF - 368 KB - 20 pages] Table 1).

Characterization of Antigen Gene Content

Figure 2

Thumbnail of Genetic diversity within capsule and major antigen encoding genes among 32 Neisseria meningitidis clonal complex 4821 isolates, China, 1972–2011. Gene sequences from 32 clonal complex 4821 isolates were compared with the epidemic reference strain 053442 (the topmost isolate on the juxtaposed core genome phylogenetic tree). Scale bar represents total substitutions per site. SNP, single-nucleotide polymorphism.

Figure 2. Genetic diversity within capsule and major antigen encoding genes among 32 Neisseria meningitidis clonal complex 4821 isolates, China, 1972–2011. Gene sequences from 32 clonal complex 4821 isolates were compared with the...

We examined 32 isolates that had undergone WGS to identify the major antigen gene (PorA VR1 and VR2, porB, FetA, nhba, and fHbp) alleles that corresponded to the epidemic strain group 1 or group 2, as determined by the core genome phylogenetic analysis (Figures 1, 2; Technical Appendix [PDF - 368 KB - 20 pages] Table 1, Figures 1–5). nadA was missing in all CC4821 study genomes. Group 1 genomes had diverse fHbp and porA alleles. All 15 group 1 isolates contained the porB 3–48 and FetA F3–3 alleles and nhba allele 124. The epidemic clone contained a few highly related porA alleles that all encoded unique, conserved porA variable regions (PorA VR1 and VR2: P1.7–2, 14). Group 2 isolates had the most antigen gene diversity, containing 5–12 different alleles for each antigen-encoding gene at the nucleotide level and no clear predominance of any single allelic profile. We observed little overlap between the antigen gene allelic profiles in groups 1 and 2. Only 2 of 13 group 1 alleles were also found in group 2. None of the antigen gene alleles found in the epidemic clone was present in group 2, suggesting that the epidemic clone had a nonoverlapping repertoire of antigens compared with historic and current group 2 isolates (Figure 2; Technical Appendix [PDF - 368 KB - 20 pages] Table 1, Figures 1–5). Group 2 was predominantly associated with PorA P1.20 variants (82%); PorB 3–229 (29%); FetA F3–9 (18%) and FetA F5–2 (18%); FHbp peptide 16 (65%); and nhba 553 (41%) (Technical Appendix [PDF - 368 KB - 20 pages] Table 1).

Recombination Events Separating Group 1 and Group 2 Isolates

We identified extensive recombination within CC4821 lineage, as detected by ClonalFrameML (20). The estimated rate of recombination relative to mutation (R/ϴ) for the entire CC4821 core genome was 1.37 (95% CI 1.40–1.34), and the relative impact of recombination to mutation (r/m) was 20.83, indicating that ≈21 nt were acquired by recombination for every point mutation within the core genome of CC4821. Sublineage specific recombination rates (R/ϴ) were 1.79 (95% CI 1.53–2.05) for group 1, 1.47 (95% CI 1.15–1.79) for group 2, and 0.5 (95% CI 0.26–0.75) for the epidemic clone.

A total of 46 recombination fragments containing 120 genes were mapped to the node that separated group 1 from group 2, indicating that these recombination events could be linked to the divergence of group 1 and 2 sublineages (Technical Appendix [PDF - 368 KB - 20 pages] Table 3). Sequence alignment and phylogenetic analysis confirmed that the 46 genomic loci represented regions of marked sequence divergence, presumably caused by homologous recombination within a common ancestral genome. These 120 recombinant genes belonged to diverse functional categories; metabolism was the overrepresented functional assignment (51/120 [42.5%]), compared with 28% of nonrecombinant genes (p<0.001 by χ2 test). Proportions of other functional groups were similar between recombinant and nonrecombinant genes.

Recombination Events Unique to the Epidemic Clone

The epidemic clone diverged from group 1 through several additional recombination events, affecting 24 genes across 7 genomic loci (Table). These genomic loci included capsule translocation genes (ctrE and ctrF); major outer membrane protein (OMP) antigen genes porA, porB, and fHbp; pts genes involved in carbohydrate transport and metabolism; and the tkl locus encoding DNA polymerase and several metabolic enzymes.

Gene Content

We identified a total of 3,292 unique genes, of which 1,730 core genes were shared by most CC4821 genomes. Eleven genes were present in >90% of group 1 isolates but missing in most group 2 isolates. Genes that were predominantly found in group 1 included lbpB encoding lactoferrin binding protein B, nhaP encoding Na+/K+ antiporter, and genes encoding several putative enzymes and 4 hypothetical proteins whose functions are unknown (Technical Appendix [PDF - 368 KB - 20 pages] Table 4). These genes were spread across several locations on the reference genome 053442, suggesting that they were acquired separately rather than in a single event. The epidemic clone had no noticeable gene gain or loss compared with other group 1 isolates.

Capsular Gene cps Cluster

The cps cluster of N. meningitidis consists of 6 regions (D-A-C-E-D′-B) required for capsule biosynthesis (region A), transport (region C), and translocation (region B) (24). The group 1 isolates were characterized by a novel capsule region A (ctrG4-cssE1-csc1-cssC3-cssB1-cssA3), associated with region C (ctrA5-ctrB1-ctrC4-ctrD1) and region E (tex-orf1-orf2) Technical Appendix [PDF - 368 KB - 20 pages] Table 1). Nine of 15 regions A-C-E identified among group 1 isolates were identical, and another 4 were almost identical (1-bp variation in csc, cssA, or ctrA over 14,489 bp). Only 2 diverged because of the allelic replacement containing the serogroup-specific polysialyltransferase gene (csb), resulting in capsular switches involving serogroups B and C (NM341215 and NM440902) (Technical Appendix [PDF - 368 KB - 20 pages] Table 2). Group 2 had substantial genetic diversity within the cps cluster, with no 2 isolates possessing an identical cps gene allelic profile (Figure 2).

Characterization of CC4821 Isolates from Shanghai

Our results classified 52 CC4821 isolates from Shanghai into epidemic clone or groups 1 and 2 on the basis of associated porA, porB, fetA, fHbp, and nhba antigen genes defined by WGS (Figure 1; online Technical Appendix Table 2). A substantial proportion of strains in Shanghai (23/52 [44%]) belonged to group 1, with 0–1 antigen gene differences (Figure 1; Technical Appendix [PDF - 368 KB - 20 pages] Table 2, Figures 1–5). Shanghai group 1 isolates were exclusively isolated during 2005–2013; were MenC (19/23), MenB (3/23), or nongroupable (1/23); and were isolated from IMD case-patients (18/23), close contacts of IMD case-patients (4/23), or an asymptomatic carrier (1/23). Within group 1, a total of 16 isolates (16/23 [70%]) contained PorA P1.7–2,14 and were consistent with the epidemic clone by all 5 antigen genes. Twenty-nine (29/52 [56%]) Shanghai isolates belonged to group 2 and differed from group 1 by 4–5 out of 5 antigen loci. Shanghai group 2 was dominated by carriage isolates (21/29 [72%]) from 1972–1977 and 2005–2013. MenB equaled or outnumbered MenC among historic (9/18 [50%]) and recent (10/11 [91%]) group 2 CC4821 isolates. MenB strains in group 1 possessed fHbp 498 or 22 encoding peptides 1.80 or 2.22 (subvariants belonging to FHbp variant group 1 or 2), as most (17/19 [89%]) of MenC epidemic clones did (Technical Appendix [PDF - 368 KB - 20 pages] Tables 1, 2).

All group 1 isolates contained fluoroquinolone-resistant gyrA allele 71 (corresponding to allele R1 in our previous study [7]), which was generated by a nonsynonymous mutation of fluoroquinolone-susceptible gyrA allele 12 (allele S1), creating an amino acid substitution of T91I. The gyrA allele 12/S1 was carried by 76% (22/29) of the group 2 strains (Technical Appendix [PDF - 368 KB - 20 pages] Table 2). In contrast, only 4 of 29 group 2 isolates from Shanghai and another 6 publicly available group 2 genomes contained resistant gyrA alleles, which were different from each other. Group 2 isolates with genotypic resistance to quinolones were genetically highly diverse as evidenced by diverse sequence types, different antigen gene profiles, and expression of both MenB and MenC capsules. However, all quinolone-resistant isolates were from 2005–2013 (Figure 1; Technical Appendix [PDF - 368 KB - 20 pages] 1, 2).


This study presents detailed genomic analyses of the current major endemic meningococcal disease lineages in China, serogroups C and B CC4821. Phylogenetic analyses have suggested that the epidemic clone (corresponding to the epidemic ChinaCC4821-R1-C/B clone described in our previous work [7)]) responsible for most recent CC4821 disease cases was nested within a distinct phylogenetic group (group 1), consistent with recent clonal expansion of a genetically distinct strain. In contrast, group 2 was temporally and genetically more diverse and accounted for a smaller proportion of IMD cases, as evidenced by the preponderance of asymptomatic carriage isolates within this group (8). Group 1 isolates were first identified during MenC outbreaks during 2003–2005, whereas group 2 isolates were mostly associated with carriage from the 1970s to 2013 (7).

This work adds to the description of CC4821 by a previous study (8) by demonstrating that, within group 1, a genetically distinct clone exists that shares the porA antigenic formula P1.7–2,14 and comprises strains associated with outbreaks and a large proportion of epidemic disease cases. This finding is in keeping with surveillance studies reporting MenC strains during 2005–2012 indicating that 55% of 238 confirmed meningococcal disease cases and 84% of 131 MenC strains belong to CC4821 with porA P1.7–2,14 (6), suggesting that the expansion of CC4821 was caused by clonal expansion of this antigenic type. Genomic analyses demonstrated that the epidemic CC4821 lineage had undergone 2 crucial genetic events compared with historic asymptomatic carriage isolates. First, CC4821 diverged into 2 major sublineages through extensive recombination events predominantly affecting genes involved in metabolic functions. Such extensive allelic exchanges might have enhanced the transmission fitness or the invasive potential within group 1. Second, a virulent, antigenically unique, epidemic strain emerged from within group 1 in a second set of more focused recombination events affecting major antigen genes, porA, porB, fHbp, capsule genes ctrE and ctrF, and several metabolic genes associated with oxidative phosphorylation and glycolytic processes. These genetic changes possibly account for the rapid dissemination and increased invasive potential of the epidemic clone.

This work also adds to the evidence that the emergence and persistence of virulent meningococcal strains occurs through the introduction of a novel antigenic variant in an immunologically naive population (9,15,2527). In addition, even though a few major antigen gene repertoires mediate the microevolution of virulent meningococcal strains, a larger and more assorted number of metabolic genes might be involved in the divergence of sublineages. Additional research is needed to elucidate the intricate interactions between various metabolic pathways in the fitness and virulence potential of meningococci.

Selection pressure of fluoroquinolones might also affect the evolution and adaption of meningococcal strains, as evidenced by various resistant gyrA alleles recovered in many meningococcal lineages and singletons only in the quinolone era in China (7,28). All of the group 1 CC4821 strains contain gyrA allele 71/R1, in contrast to various resistant gyrA alleles in group 2, suggesting acquisition of this trait by group 1 at an early stage in evolution. The gyrA allele 71/R1 derived from gyrA allele 12/S1, which was possessed by most of the group 2 strains both in prequinolone and quinolone eras, indicating their common origin. The precise role of quinolone resistance in the emergence of the epidemic clone requires further study.

The relative rates and impact of recombination within the CC4821 lineage were considerably higher than previous estimates that examined recombination across 7 housekeeping genes (2931). We identified 4,026 recombination events and 21 recombinant single-nucleotide polymorphisms (SNPs) for every point mutation within CC4821. Estimates of meningococcal recombination rates from whole-genome sequence data are limited. A study of MenA CC5 strains in Africa (27), an epidemic lineage notable for relatively low genetic diversity, found 34 recombinant sequences and 12 recombinant SNPs for every point mutation. A study of recombination rates within conserved multilocus sequence type loci among carriage isolates from the Czech Republic in 1993 found 6.2–16.8 SNPs caused by recombination for each point mutation (29). These data suggest that although homologous recombination is a shared mechanism for meningococcal strain emergence and persistence (29,3133), the frequency and extent of recombination likely differs substantially between lineages.

High rates of recombination have led to multiple distinct capsular switch strains expressing group B and C capsule, as described in this study, and serogroup W (8,34). Both MenB and MenC were phylogenetically diverse and interspersed within both groups 1 and 2, suggesting multiple distinct capsular switches between these serogroups rather than clonal expansion of a single capsular switch strain (35,36). Moreover, recombination events that led to no apparent change in capsular phenotype were also prevalent.

Wide geographic and temporal spread of MenB strains in China is of concern given that vaccines currently in use do not protect against MenB disease. Furthermore, marked heterogeneity exists among both MenB and MenC strains within the gene encoding FHbp, a key component of MenB vaccines.

At the time of its emergence, CC4821 was confined to China and Taiwan (37). However, review of the PubMLST Neisseria genome database suggests a recent increase in the geographic spread of CC4821. CC4821 isolates were reported from a small number of asymptomatic carriers in Brazil (2014) (38,39) and Australia (2012–2017), as well as isolated disease cases and carriers from France (2009–2011), the United Kingdom (2011–2014), the United States (2007–2016), and India (2017) (14). A recent case report of a quinolone-resistant sequence type 4821 strain from Japan in a patient with no history of foreign travel also suggests local transmission of CC4821 in Japan (40). Global CC4821 carriage and IMD isolates on the PubMLST Neisseria database were antigenically diverse, containing PorA types associated with both group 1 and 2 strains, and expressed both serogroup C and B capsules. This finding suggests low-level dissemination of CC4821 strains with diverse virulence and antigenic types as opposed to clonal spread of a single epidemic strain. This pattern is in contrast to the pandemic spread of highly clonal serogroup A epidemic strains from China from the 1960s through the 1990s (25). Genome-based surveillance of these global CC4821 strains is needed to monitor the global spread of this clonal lineage.

Our study is limited by lack of isolates before the 1970s and during 1986–2004; isolates from those periods might have provided a clearer picture of the multiple evolutionary steps that led to the epidemic CC4821 clone. Also, the high frequency of recombination within capsular genes makes it difficult to accurately determine the direction of capsular switch.

In summary, we have presented detailed genomic analysis of a major hypervirulent meningococcal lineage associated with MenC and MenB in China and identified key genomic factors that might have led to the emergence and persistence of MenC in China. The potential emergence of MenB is of public health concern. Strengthened laboratory surveillance for disease cases and carefully planned carriage surveys are needed to monitor global trends, detect outbreaks, and inform immunization policies.

Dr. Guo is an associate professor in the Institute of Antibiotics, Huashan Hospital, Fudan University. Her research interests include mechanisms of antimicrobial resistance in clinical isolates.



We thank Jane Marsh for providing thoughtful comments on an early draft of this manuscript.

This study made use of the PubMLST Neisseria genome database (, developed by Keith Jolley and hosted by the University of Oxford. The development of this site has been funded by the Wellcome Trust and European Union. This study also made use of the Meningitis Research Foundation Meningococcus Genome Library (, developed by Public Health England, the Wellcome Trust Sanger Institute, and the University of Oxford as a collaboration. The project is funded by the Meningitis Research Foundation.

This work was supported by grants from the National Natural Science Foundation of China (81673479, 81120108024, and 81601801); the Shanghai Pujiang Program (16PJD010); the Shanghai Rising-Star Program (17QA1403100), the Shanghai Medical Health Plan for Outstanding Young Talents (2017YQ039); and the 4th Three-Year Action Plan for Public Health of Shanghai (GWTD2015S01).



  1. Halperin SA, Bettinger JA, Greenwood B, Harrison LH, Jelfs J, Ladhani SN, et al. The changing and dynamic epidemiology of meningococcal disease. Vaccine. 2012;30(Suppl 2):B2636. DOIPubMed
  2. Zhang X, Shao Z, Yang E, Xu L, Xu X, Li M, et al. Molecular characterization of serogroup C Neisseria meningitidis isolated in China. J Med Microbiol. 2007;56:12249. DOIPubMed
  3. Ni JD, Jin YH, Dai B, Wang XP, Liu DQ, Chen X, et al. Recent epidemiological changes in meningococcal disease may be due to the displacement of serogroup A by serogroup C in Hefei City, China. Postgrad Med J. 2008;84:8792. DOIPubMed
  4. Shao Z, Li W, Ren J, Liang X, Xu L, Diao B, et al. Identification of a new Neisseria meningitidis serogroup C clone from Anhui province, China. Lancet. 2006;367:41923. DOIPubMed
  5. Peng J, Yang L, Yang F, Yang J, Yan Y, Nie H, et al. Characterization of ST-4821 complex, a unique Neisseria meningitidis clone. Genomics. 2008;91:7887. DOIPubMed
  6. Zhou H, Shan X, Sun X, Xu L, Gao Y, Li M, et al. Clonal characteristics of invasive Neisseria meningitidis following initiation of an A + C vaccination program in China, 2005-2012. J Infect. 2015;70:3743. DOIPubMed
  7. Chen M, Guo Q, Wang Y, Zou Y, Wang G, Zhang X, et al. Shifts in the antibiotic susceptibility, serogroups, and clonal complexes of Neisseria meningitidis in Shanghai, China: a time trend analysis of the pre-quinolone and quinolone eras. PLoS Med. 2015;12:e1001838, discussion e1001838. DOIPubMed
  8. Zhu B, Xu Z, Du P, Xu L, Sun X, Gao Y, et al. Sequence type 4821 clonal complex serogroup B Neisseria meningitidis in China, 1978–2013. Emerg Infect Dis. 2015;21:92532. DOIPubMed
  9. Harrison LH, Jolley KA, Shutt KA, Marsh JW, O’Leary M, Sanza LT, et al.; Maryland Emerging Infections Program. Antigenic shift and increased incidence of meningococcal disease. J Infect Dis. 2006;193:126674. DOIPubMed
  10. Thompson EA, Feavers IM, Maiden MC. Antigenic diversity of meningococcal enterobactin receptor FetA, a vaccine component. Microbiology. 2003;149:184958. DOIPubMed
  11. Urwin R, Fox AJ, Musilek M, Kriz P, Maiden MC. Heterogeneity of the PorB protein in serotype 22 Neisseria meningitidis. J Clin Microbiol. 1998;36:36802.PubMed
  12. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:45577. DOIPubMed
  13. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:20689. DOIPubMed
  14. Jolley KA, Maiden MC. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595. DOIPubMed
  15. Mustapha MM, Marsh JW, Krauland MG, Fernandez JO, de Lemos APS, Dunning Hotopp JC, et al. Genomic epidemiology of hypervirulent serogroup W, ST-11 Neisseria meningitidis. EBioMedicine. 2015;2:144755. DOIPubMed
  16. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:13123. DOIPubMed
  17. Darling AE, Mau B, Perna NT. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5:e11147. DOIPubMed
  18. Parkhill J, Achtman M, James KD, Bentley SD, Churcher C, Klee SR, et al. Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491. Nature. 2000;404:5026. DOIPubMed
  19. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:27319. DOIPubMed
  20. Didelot X, Wilson DJ. ClonalFrameML: efficient inference of recombination in whole bacterial genomes. PLOS Comput Biol. 2015;11:e1004041. DOIPubMed
  21. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43:e15. DOIPubMed
  22. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:36913. DOIPubMed
  23. Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:336. DOIPubMed
  24. Harrison OB, Claus H, Jiang Y, Bennett JS, Bratcher HB, Jolley KA, et al. Description and nomenclature of Neisseria meningitidis capsule locus. Emerg Infect Dis. 2013;19:56673. DOIPubMed
  25. Zhu P, van der Ende A, Falush D, Brieske N, Morelli G, Linz B, et al. Fit genotypes and escape variants of subgroup III Neisseria meningitidis during three pandemics of epidemic meningitis. Proc Natl Acad Sci U S A. 2001;98:52349. DOIPubMed
  26. Taha MK, Giorgini D, Ducos-Galand M, Alonso JM. Continuing diversification of Neisseria meningitidis W135 as a primary cause of meningococcal disease after emergence of the serogroup in 2000. J Clin Microbiol. 2004;42:415863. DOIPubMed
  27. Lamelas A, Harris SR, Röltgen K, Dangy JP, Hauser J, Kingsley RA, et al. Emergence of a new epidemic Neisseria meningitidis serogroup A Clone in the African meningitis belt: high-resolution picture of genomic changes that mediate immune evasion. MBio. 2014;5:e0197414. DOIPubMed
  28. Zhu B, Fan Y, Xu Z, Xu L, Du P, Gao Y, et al. Genetic diversity and clonal characteristics of ciprofloxacin-resistant meningococcal strains in China. J Med Microbiol. 2014;63:14118. DOIPubMed
  29. Jolley KA, Wilson DJ, Kriz P, McVean G, Maiden MC. The influence of mutation, recombination, population history, and selection on patterns of genetic diversity in Neisseria meningitidis. Mol Biol Evol. 2005;22:5629. DOIPubMed
  30. Didelot X, Maiden MC. Impact of recombination on bacterial evolution. Trends Microbiol. 2010;18:31522. DOIPubMed
  31. Vos M, Didelot X. A comparison of homologous recombination rates in bacteria and archaea. ISME J. 2009;3:199208. DOIPubMed
  32. Budroni S, Siena E, Dunning Hotopp JC, Seib KL, Serruto D, Nofroni C, et al. Neisseria meningitidis is structured in clades associated with restriction modification systems that modulate homologous recombination. Proc Natl Acad Sci U S A. 2011;108:44949. DOIPubMed
  33. Kong Y, Ma JH, Warren K, Tsang RS, Low DE, Jamieson FB, et al. Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis. Genome Biol Evol. 2013;5:161127. DOIPubMed
  34. He B, Jia Z, Zhou H, Wang Y, Jiang X, Ma H, et al. CC4821 serogroup W meningococcal disease in China. Int J Infect Dis. 2014;29:1134. DOIPubMed
  35. Lucidarme J, Lekshmi A, Parikh SR, Bray JE, Hill DM, Bratcher HB, et al. Frequent capsule switching in ‘ultra-virulent’ meningococci - Are we ready for a serogroup B ST-11 complex outbreak? J Infect. 2017;75:95103. DOIPubMed
  36. Mustapha MM, Marsh JW, Krauland MG, Fernandez JO, de Lemos AP, Dunning Hotopp JC, et al. Genomic investigation reveals highly conserved, mosaic, recombination events associated with capsular switching among invasive Neisseria meningitidis serogroup W sequence type (ST)-11 strains. Genome Biol Evol. 2016;8:206575. DOIPubMed
  37. Chiou CS, Liao JC, Liao TL, Li CC, Chou CY, Chang HL, et al. Molecular epidemiology and emergence of worldwide epidemic clones of Neisseria meningitidis in Taiwan. BMC Infect Dis. 2006;6:25. DOIPubMed
  38. Nunes AM, Ribeiro GS, Ferreira IE, Moura AR, Felzemburgh RD, de Lemos AP, et al. Meningococcal carriage among adolescents after mass meningococcal c conjugate vaccination campaigns in Salvador, Brazil. PLoS One. 2016;11:e0166475. DOIPubMed
  39. Moura ARSS, Kretz CB, Ferreira IE, Nunes AMPB, de Moraes JC, Reis MG, et al. Molecular characterization of Neisseria meningitidis isolates recovered from 11-19-year-old meningococcal carriers in Salvador, Brazil. PLoS One. 2017;12:e0185038. DOIPubMed
  40. Kawasaki Y, Matsubara K, Takahashi H, Morita M, Ohnishi M, Hori M, et al. Invasive meningococcal disease due to ciprofloxacin-resistant Neisseria meningitidis sequence type 4821: the first case in Japan. J Infect Chemother. 2017 Dec 7:S1341–321X(17)30274-X [Epub ahead of print].




Technical Appendix


Cite This Article

DOI: 10.3201/eid2404.171744

1These authors contributed equally to this article.

Table of Contents – Volume 24, Number 4—April 2018


Please use the form below to submit correspondence to the authors or contact them at the following address:

Min Chen, Shanghai Municipal Center for Disease Control and Prevention, 1380 W ZhongShan Rd, Shanghai, China; ; Yohei Doi, Division of Infectious Diseases, University of Pittsburgh School of Medicine, S829 Scaife Hall, 3550 Terrace St, Pittsburgh, PA 15261, USA; ; Minggui Wang, Institute of Antibiotics, Huashan Hospital, Fudan University, 12 M Wulumuqi Rd, Shanghai 200040, China

character(s) remaining.

Comment submitted successfully, thank you for your feedback.