Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 Oct;7(10):000662.
doi: 10.1099/mgen.0.000662.

Modelling evolutionary pathways for commensalism and hypervirulence in Neisseria meningitidis

Affiliations

Modelling evolutionary pathways for commensalism and hypervirulence in Neisseria meningitidis

Christopher A Mullally et al. Microb Genom. 2021 Oct.

Abstract

Neisseria meningitidis, the meningococcus, resides exclusively in humans and causes invasive meningococcal disease (IMD). The population of N. meningitidis is structured into stable clonal complexes by limited horizontal recombination in this naturally transformable species. N. meningitidis is an opportunistic pathogen, with some clonal complexes, such as cc53, effectively acting as commensal colonizers, while other genetic lineages, such as cc11, are rarely colonizers but are over-represented in IMD and are termed hypervirulent. This study examined theoretical evolutionary pathways for pathogenic and commensal lineages by examining the prevalence of horizontally acquired genomic islands (GIs) and loss-of-function (LOF) mutations. Using a collection of 4850 genomes from the BIGSdb database, we identified 82 GIs in the pan-genome of 11 lineages (10 hypervirulent and one commensal lineage). A new computational tool, Phaser, was used to identify frameshift mutations, which were examined for statistically significant association with genetic lineage. Phaser identified a total of 144 frameshift loci of which 105 were shown to have a statistically significant non-random distribution in phase status. The 82 GIs, but not the LOF loci, were associated with genetic lineage and invasiveness using the disease carriage ratio metric. These observations have been integrated into a new model that infers the early events of the evolution of the human adapted meningococcus. These pathways are enriched for GIs that are involved in modulating attachment to the host, growth rate, iron uptake and toxin expression which are proposed to increase competition within the meningococcal population for the limited environmental niche of the human nasopharynx. We surmise that competition for the host mucosal surface with the nasopharyngeal microbiome has led to the selection of isolates with traits that enable access to cell types (non-phagocytic and phagocytic) in the submucosal tissues leading to an increased risk for IMD.

Keywords: clonal complexes; commensalism; competition; pathogenicity.

PubMed Disclaimer

Conflict of interest statement

The authors declare that there are no conflicts of interest.

Figures

Fig. 1.
Fig. 1.
Location of 82 genetic islands (GIs) in 11 representative genomes from different clonal complexes. GIs that were associated with at least one clonal complex are marked numerically on the outer ring from 1 to 63. Where there are multiple islands at one locus the islands are labelled alphabetically as islands A, B, C and D. These flexible loci bring the total to 82 GIs. The closed genome FAM18 was used as the reference genome and is shown as the innermost ring. The capsule island is labelled in green and was not assigned an island number. The image was generated using BRIG [24]. The innermost black trace shows G+C content.
Fig. 2.
Fig. 2.
Prevalence of 82 genomic islands within 11 meningococcal clonal complexes. The prevalence of each island within each lineage is shown in blue. The function of each island is shown in a different colour. The presence of IS elements, repeats flanking the island and atypical G+C content is also shown. If a cell is white, the island or feature is absent. The GIs were grouped according to functionality and ranked within each functional group by prevalence. The clonal complexes were grouped based on hierarchical clustering of the presence and absence of the GIs, and a corresponding dendrogram was generated. GG-I is shown in orange and GG-II is shown in green. Outgroup cc53 can be seen on the left hand side of the tree.
Fig. 3.
Fig. 3.
Disease/carriage ratio of each pathogenic clonal complex plotted against all 82 GIs of the accessory genome (a) and hypervirulence-associated genes (HVAGs) consisting of 35 GIs and 26 genes (b). The disease/carriage ratio was an in silico calculation (see Methods and Table 1). GG-I is shaded in orange while GG-II is shaded in green.
Fig. 4.
Fig. 4.
Two-way hierarchical clustering showing the association of hypervirulence-associated genes (HVAGs) with all hypervirulent lineages. Clonal complexes are hierarchically clustered based on the presence of the HVAGs and a corresponding two-way dendrogram was generated. Where a GI has been assigned a number, it has been listed next to the genomic island. Genomic islands comprising hypothetical genes have been denoted as HI 1 to 5.
Fig. 5.
Fig. 5.
A decision tree model illustrating the putative relationships of modern clonal complexes with the acquisition and loss of genomic islands in ancestral N. meningitidis . The combined genomic island dataset from Figs 2 and 4 was input into MrBayes to draw a Bayesian tree (see Methods for parameters). Each node shows the Bayesian posterior probability (0 to 1.0) of the split into the different clusters of GIs that are associated with a clonal complex. A node with posterior probability ≥0.8 is considered to be highly significant. At each node, the GIs that contributed to the formation of each group are shown. The far left hand node is the putative ancestral split between the commensal cc53 (blue text) and the early specialists (yellow text) that evolved to become genogroups I (red text) and II (green text). The GIs that are characteristic of cc11 are indicated in purple text. GIs characteristic of the ancestral node are coloured according to lineage. Absence of GIs in any given lineage relative to the origin node is shown in a solid black box. This model cannot distinguish between whether these GIs were not acquired at the ancestral node by the founding progenitor of the lineage or that the GI was acquired and subsequently lost during the evolution of the clonal complex. Created using Biorender.com. 1 = NEIS1012/2064/2583/2380/2404/2854/2551/2910/2789/2406/2536/ 2537. 2 = NEIS2778/2463/0955/0881/1861/1866/1715/0223/0445. 3 = NEIS1866/1715/0223/0445. 4 = NEIS/2778/2463/0445.
Fig. 6.
Fig. 6.
Heatmap of phase state of LOF loci by SSR (a) and indels (b). Clonal complexes were hierarchically clustered based on similarity of the phase state of each LOF loci and a corresponding dendrogram was generated. If a gene is always intact, it is shown in dark blue. If it is always interrupted, it is shown in white. If the gene is absent in the clonal complex, it is shaded grey. The phase state for a gene was considered significantly different between each clonal complex and cc53 if the corrected p-value was <0.05.
Fig. 7.
Fig. 7.
Putative biological functions of ancestral meningococci: commensals, early specialists and modern IMD lineage cc11. (a) Proposed stages in the evolution of the meningococcus. The early meningococcus, the commensal lineage (cc53, blue), an early specialist population (yellow node), GG-I (red), GG-II (green) and cc-11 (mauve). These colour codes correspond to the nodes in Fig. 5. Known virulence determinants are represented in the bottom panel by icons which appear in relationship to the cartoon representation of the meningococcus. Pilin glycan is represented by two different colour codes corresponding to: Hex(2-4)-GATDH O-linked glycan (blue) and O-acetyl-Hex (1-2) diNAcBac glycan (green). A yellow circle represents O-acetylation of the glycan. Type IV pili are shown as two classes: class I (blue) and II (purple). (b) Proposed interactions of cc53, early specialists and cc11 with host cells derived from the known functions of GIs enriched in each population (Fig. 5). Abbreviations: DC = dendritic cells, Mφ = macrophages, Nφ = neutrophils. Figure created using Biorender.com.

Similar articles

Cited by

References

    1. Stephens DS, Greenwood B, Brandtzaeg P. Epidemic meningitis, meningococcaemia, and Neisseria meningitidis . Lancet. 2007;369:2196–2210. doi: 10.1016/S0140-6736(07)61016-2. - DOI - PubMed
    1. Johswich K. Innate immune recognition and inflammation in Neisseria meningitidis infection. Pathog Dis. 2017;75 doi: 10.1093/femspd/ftx022. - DOI - PubMed
    1. Christensen H, May M, Bowen L, Hickman M, Trotter CL. Meningococcal carriage by age: a systematic review and meta-analysis. Lancet Infect Dis. 2010;10:853–861. doi: 10.1016/S1473-3099(10)70251-6. - DOI - PubMed
    1. Peterson ME, Mile R, Li Y, Nair H, Kyaw MH. Meningococcal carriage in high-risk settings: A systematic review. Int J Infect Dis. 2018;73:109–117. doi: 10.1016/j.ijid.2018.05.022. - DOI - PubMed
    1. Bidmos FA, Neal KR, Oldfield NJ, Turner DP, Ala’Aldeen DA, et al. Persistence, replacement, and rapid clonal expansion of meningococcal carriage isolates in a 2008 university student cohort. J Clin Microbiol. 2011;49:506–512. doi: 10.1128/JCM.01322-10. - DOI - PMC - PubMed

Publication types

Substances