Extensive genomic variation within clonal complexes of Neisseria meningitidis

Meningococcal disease is a widely distributed complex disease affecting all age categories. It can cause severe meningitis and septicemia, especially in unvaccinated infants and young children. The causative agent, Neisseria meningitidis (Nm), can be phenotypically and genetically differentiated into serogroups and sequence types (STs) and has a highly dynamic population structure. To obtain a deeper understanding of the epidemiology of Nm, we sequenced seven Nm genomes. Large-scale genomic analysis was conducted with these 7 Nm genomes, 27 additional Nm genomes from GenBank, and 4 other Neisseria genomes. We observed extensive homologous recombination in all gene functional categories among different Nm genomes. Homologous recombination is so frequent that it has resulted in numerous chimeric open reading frames, including genes in the capsule biosynthesis cluster and loci targeted by commercial vaccines. Our results reveal that, despite widespread use, evolutionary relationships inferred from the standard seven-gene multilocus sequence typing (MLST) method could not predict virulence gene content or strain phenotype. In fact, up to 28% of the virulence-associated genes could differ between strains of identical STs. Consistent with previous studies, we found that allelic recombination is also associated with alterations in antibiotic susceptibility. Overall, these findings emphasize the extensive genomic plasticity of Nm and the limitations of standard molecular methods to quantify this genotypic and phenotypic diversity.