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

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Metagenomic data from cerebrospinal fluid permits tracing the origin and spread of Neisseria meningitidis CC4821 in China

Metagenomic next-generation sequencing (mNGS) is useful for difficult to cultivate pathogens. Here, we use cerebrospinal fluid mNGS to diagnose invasive meningococcal disease. The complete genome sequences of Neisseria meningitidis were assembled using N. meningitidis of ST4821-serotype C isolated from four patients. To investigate the phylogeny, 165 CC4821 N. meningitidis genomes from 1972 to 2017 were also included. The core genome accumulated variation at a rate of 4.84×10⁻⁸ substitutions/nucleotide site/year. CC4821 differentiated into four sub-lineages during evolution (A, B, C, and D). While evolving from sub-lineage A (early stage) to sub-lineage D (late stage), the ST and CC4821 serotype converged into the ST4821-serotype C clone. Most strains of sub-lineage D were isolated from invasive meningococcal disease, with increasing resistance to quinolones. Phylogeographic analysis suggests that CC4821 has spread across 14 countries. Thus, the selective pressure of quinolones may cause CC4821 to converge evolutionarily, making it more invasive and facilitating its spread.

Metagenomic data from cerebrospinal fluid was used to genotype Neisseria meningitidis in patients with invasive meningococcal disease and trace the origin of the pathobiont, providing a phylogeographic analysis of the strain’s evolution in China and its global spread.

Serogroup Y Clonal Complex 23 Meningococcus in China Acquiring Penicillin Resistance from Commensal Neisseria lactamica Species

Invasive meningococcal disease (IMD) due to serogroup Y Neisseria meningitidis (NmY) is rare in China; recently, an invasive NmY isolate, Nm512, was discovered in Shanghai with decreased susceptibility to penicillin (PenNS). Here, we investigated the epidemiology of NmY isolates in Shanghai and explored the potential commensal Neisseria lactamica donor of the PenNS NmY isolate. A total of 491 N. meningitidis and 724 commensal Neisseria spp. isolates were collected. Eleven NmY isolates were discovered from IMD (n = 1) and carriers (n = 10), including two PenNS isolates with five-key-mutation-harboring (F504L-A510V-I515V-H541N-I566V) penA genes. Five of the eight ST-175 complex (CC175) isolates had a genotype [Y:P1.5-1,2-2:F5-8:ST-175(CC175)] identical to that of the predominant invasive clone found in South Africa. Only one invasive NmY CC23 isolate (Nm512) was discovered; this isolate carried a novel PenNSpenA832 allele, which was identified in commensal N. lactamica isolates locally. Recombination analysis and transformation of the penA allele highlighted that N. meningitidis Nm512 may acquire resistance from its commensal donor; this was supported by the similar distribution of transformation-required DNA uptake sequence variants and the highly cognate receptor ComP between N. meningitidis and N. lactamica. In 2,309 NmY CC23 genomes from the PubMLST database, isolates with key-mutation-harboring penA genes comprised 12% and have been increasing since the 1990s, accompanied by recruitment of the blaROB-1 and/or quinolone resistance allele. Moreover, penA22 was predominant among genomes without key mutations in penA. These results strongly suggest that Nm512 is a descendant of the penA22-harboring CC23 isolate from Europe and acquired its penicillin resistance locally from commensal N. lactamica species by natural transformation.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Evolution of Sequence Type 4821 Clonal Complex Hyperinvasive and Quinolone-Resistant Meningococci

Expansion of quinolone-resistant Neisseria meningitidis clone China^{CC4821-R1-C/B} from sequence type (ST) 4821 clonal complex (CC4821) caused a serogroup shift from serogroup A to serogroup C invasive meningococcal disease (IMD) in China. To determine the relationship among globally distributed CC4821 meningococci, we analyzed whole-genome sequence data from 173 CC4821 meningococci isolated from 4 continents during 1972–2019. These meningococci clustered into 4 sublineages (1–4); sublineage 1 primarily comprised of IMD isolates (41/50, 82%). Most isolates from outside China (40/49, 81.6%) formed a distinct sublineage, the Europe–USA cluster, with the typical strain designation B:P1.17-6,23:F3-36:ST-3200(CC4821), harboring mutations in penicillin-binding protein 2. These data show that the quinolone-resistant clone China^{CC4821-R1-C/B} has expanded to other countries. The increasing distribution worldwide of serogroup B CC4821 raises the concern that CC4821 has the potential to cause a pandemic that would be challenging to control, despite indirect evidence that the Trumenba vaccine might afford some protection.

Meningococcal vaccines in China

Meningococcal meningitis caused by Neisseria meningitidis is a reportable infectious disease in China, due to the high incidence of meningitis in the era before the availability of vaccines. The disease incidence was markedly reduced after meningococcal vaccination was introduced in the 1980s. Currently, there are polysaccharide, conjugate, and combined vaccine formulations against meningococcal meningitis in the Chinese market, almost all of which are produced by domestic manufacturers. It is necessary to further enhance national meningococcal surveillance to improve the level of prevention and control of meningococcus. However, the immune efficacy and persistence of immunity of vaccines should be monitored. More importantly, additional investments should be made to develop serogroup B meningococcal vaccines.

Genomic analysis of the meningococcal ST-4821 complex-Western clade, potential sexual transmission and predicted antibiotic susceptibility and vaccine coverage

INTRODUCTION

The ST-4821 complex (cc4821) is a leading cause of serogroup C and serogroup B invasive meningococcal disease in China where diverse strains in two phylogenetic groups (groups 1 and 2) have acquired fluoroquinolone resistance. cc4821 was recently prevalent among carriage isolates in men who have sex with men in New York City (USA). Genome-level population studies have thus far been limited to Chinese isolates. The aim of the present study was to build upon these with an extended panel of international cc4821 isolates.

METHODS

Genomes of isolates from Asia (1972 to 2017), Europe (2011 to 2018), North America (2007), and South America (2014) were sequenced or obtained from the PubMLST Neisseria database. Core genome comparisons were performed in PubMLST.

RESULTS

Four lineages were identified. Western isolates formed a distinct, mainly serogroup B sublineage with alleles associated with fluoroquinolone susceptibility (MIC <0.03 mg/L) and reduced penicillin susceptibility (MIC 0.094 to 1 mg/L). A third of these were from anogenital sites in men who have sex with men and had unique denitrification gene alleles. Generally 4CMenB vaccine strain coverage was reliant on strain-specific NHBA peptides.

DISCUSSION

The previously identified cc4821 group 2 was resolved into three separate lineages. Clustering of western isolates was surprising given the overall diversity of cc4821. Possible association of this cluster with the anogenital niche is worthy of monitoring given concerns surrounding antibiotic resistance and potential subcapsular vaccine escape.

Toward a Global Genomic Epidemiology of Meningococcal Disease

Whole-genome sequencing (WGS) is invaluable for studying the epidemiology of meningococcal disease. Here we provide a perspective on the use of WGS for meningococcal molecular surveillance and outbreak investigation, where it helps to characterize pathogens, predict pathogen traits, identify emerging pathogens, and investigate pathogen transmission during outbreaks. Standardization of WGS workflows has facilitated their implementation by clinical and public health laboratories (PHLs), but further development is required for metagenomic shotgun sequencing and targeted sequencing to be widely available for culture-free characterization of bacterial meningitis pathogens. Internet-accessible servers are being established to support bioinformatics analysis, data management, and data sharing among PHLs. However, establishing WGS capacity requires investments in laboratory infrastructure and technical knowledge, which is particularly challenging in resource-limited regions, including the African meningitis belt. Strategic WGS implementation is necessary to monitor the molecular epidemiology of meningococcal disease in these regions and construct a global view of meningococcal disease epidemiology.

Evolutionary analysis of gyrA gene from Neisseria meningitidis bacterial strains of clonal complex 4821 collected in China between 1978 and 2016

BACKGROUND

Neisseria meningitidis (N.meningitidis) bacteria belonging to clonal complex 4821 (CC4821) have been mainly reported in China and have been characterized by a high resistance rate to ciprofloxacin (CIP). The aim of this study was to assess the evolution of the DNA gyrase A (gyrA) gene from N.meningitidis CC4821 strains collected in China between 1978 and 2016. The complete sequence of gyrA gene from 77 strains are reported in this study and analyzed in the context of publicly available sequences from N. meningitidis of other CCs as well as other Neisseria species.

RESULTS

The phylogenetic analysis of CC4821 gyrA gene reveals at least 5 distinct genetic clusters. These clusters are not CC4821-specific showing that gyrA evolution is independent of CC4821 evolution. Some clusters contain sequences from other Neisseria species. Recombination within N.meningitidis strains and between Neisseria species was identified in SimPlot analysis. Finally, amino acid substitutions within GyrA protein were analyzed. Only one position, 91 (83 in E.coli gyrA gene), was linked to CIP resistance. Thirty-one additional putative resistance markers were identified, as amino acid substitutions were only found in resistant strains.

CONCLUSIONS

The evolution of gyrA gene of CC4821 N.meningitidis strains is not dependent on CC4821 evolution or on CIP resistance phenotype. Only amino acid 91 is linked to CIP resistance phenotype. Finally, recombination inter- and intra-species is likely to result in the acquisition of various resistance markers, 31 of them being putatively mapped in the present study. Analyzing the evolution of gyrA gene within CC4821 strains is critical to monitor the CIP resistance phenotype and the acquisition of new resistance markers. Such studies are necessary for the control of the meningococcal disease and the development of new drugs targeting DNA gyrase.

Meningococcal Quinolone Resistance Originated from Several Commensal Neisseria Species

Quinolone resistance is increasing in Neisseria meningitidis, with its prevalence in China being high (>70%), but its origin remains unknown. The aim of this study was to investigate the donors of mutation-harboring gyrA alleles in N. meningitidis A total of 198 N. meningitidis isolates and 293 commensal Neisseria isolates were collected between 2005 and 2018 in Shanghai, China. The MICs of ciprofloxacin were determined using the agar dilution method. The resistance-associated genes gyrA and parC were sequenced for all isolates, while a few isolates were sequenced on the Illumina platform. The prevalences of quinolone resistance in the N. meningitidis and commensal Neisseria isolates were 67.7% (134/198) and 99.3% (291/293), respectively. All 134 quinolone-resistant N. meningitidis isolates possessed mutations in T91 (n = 123) and/or D95 (n = 12) of GyrA, with 7 isolates also harboring ParC mutations and exhibiting higher MICs. Phylogenetic analysis of the gyrA sequence identified six clusters. Among the 71 mutation-harboring gyrA alleles found in 221 N. meningitidis isolates and genomes (n = 221), 12 alleles (n = 103, 46.6%) were included in the N. meningitidis cluster, while 20 alleles (n = 56) were included in the N. lactamica cluster, 27 alleles (n = 49) were included in the N. cinerea cluster, and 9 alleles (n = 10) were included in the N. subflava cluster. Genomic analyses identified the exact N. lactamica donors of seven mutation-harboring gyrA alleles (gyrA92, gyrA97, gyrA98, gyrA114, gyrA116, gyrA151, and gyrA230) and the N. subflava donor isolate of gyrA171, with the sizes of the recombinant fragments ranging from 634 to 7,499 bp. Transformation of gyrA fragments from these donor strains into a meningococcal isolate increased its ciprofloxacin MIC from 0.004 μg/ml to 0.125 or 0.19 μg/ml and to 0.5 μg/ml with further transformation of an additional ParC mutation. Over half of the quinolone-resistant N. meningitidis isolates acquired resistance by horizontal gene transfer from three commensal Neisseria species. Quinolone resistance in N. meningitidis increases in a stepwise manner.