Meningococcal Disease-Associated Prophage-Like Elements Are Present in Neisseria gonorrhoeae and Some Commensal Neisseria Species

Manipulation of natural transformation by AbaR-type islands promotes fixation of antibiotic resistance in Acinetobacter baumannii

The opportunistic pathogen Acinetobacter baumannii, carries variants of A. baumannii resistance islands (AbaR)-type genomic islands conferring multidrug resistance. Their pervasiveness in the species has remained enigmatic. The dissemination of AbaRs is intricately linked to their horizontal transfer via natural transformation, a process through which bacteria can import and recombine exogenous DNA, effecting allelic recombination, genetic acquisition, and deletion. In experimental populations of the closely related pathogenic Acinetobacter nosocomialis, we quantified the rates at which these natural transformation events occur between individuals. When integrated into a model of population dynamics, they lead to the swift removal of AbaRs from the population, contrasting with the high prevalence of AbaRs in genomes. Yet, genomic analyses show that nearly all AbaRs specifically disrupt comM, a gene encoding a helicase critical for natural transformation. We found that such disruption impedes gene acquisition, and deletion, while moderately impacting acquisition of single nucleotide polymorphism. A mathematical evolutionary model demonstrates that AbaRs inserted into comM gain a selective advantage over AbaRs inserted in sites that do not inhibit or completely inhibit transformation, in line with the genomic observations. The persistence of AbaRs can be ascribed to their integration into a specific gene, diminishing the likelihood of their removal from the bacterial genome. This integration preserves the acquisition and elimination of alleles, enabling the host bacterium-and thus its AbaR-to adapt to unpredictable environments and persist over the long term. This work underscores how manipulation of natural transformation by mobile genetic elements can drive the prevalence of multidrug resistance.

The potential for bacteriophages and prophage elements in fighting and preventing the gonorrhea

Bacteriophages are the most numerous entities on earth and are found everywhere their bacterial hosts live. As natural bacteria killers, phages are extensively investigated as a potential cure for bacterial infections. Neisseria gonorrhoeae (the gonococcus) is the etiologic agent of a sexually transmitted disease: gonorrhea. The rapid increase of resistance of N. gonorrhoeae to antibiotics urges scientists to look for alternative treatments to combat gonococcal infections. Phage therapy has not been tested as an anti-gonococcal therapy so far. To date, no lytic phage has been discovered against N. gonorrhoeae. Nevertheless, gonococcal genomes contain both dsDNA and ssDNA prophages, and viral particle induction has been documented. In this review, we consider literature data about the attempts of hunting for a bacteriophage specific for gonococci - the gonophage. We also discuss the potential application of prophage elements in the fight against N. gonorrhoeae. Temperate phages may be useful in preventing and treating gonorrhea as a scaffold for anti-gonococcal vaccine development and as a source of lytic enzymes with anti-gonococcal activity.

Genetic variants linked to the phenotypic outcome of invasive disease and carriage of Neisseria meningitidis

Neisseria meningitidis can be a human commensal in the upper respiratory tract but is also capable of causing invasive diseases such as meningococcal meningitis and septicaemia. No specific genetic markers have been detected to distinguish carriage from disease isolates. The aim here was to find genetic traits that could be linked to phenotypic outcomes associated with carriage versus invasive N. meningitidis disease through a bacterial genome-wide association study (GWAS). In this study, invasive N. meningitidis isolates collected in Sweden (n=103) and carriage isolates collected at Örebro University, Sweden (n=213) 2018-2019 were analysed. The GWAS analysis, treeWAS, was applied to single-nucleotide polymorphisms (SNPs), genes and k-mers. One gene and one non-synonymous SNP were associated with invasive disease and seven genes and one non-synonymous SNP were associated with carriage isolates. The gene associated with invasive disease encodes a phage transposase (NEIS1048), and the associated invasive SNP glmU S373C encodes the enzyme N-acetylglucosamine 1-phosphate (GlcNAC 1-P) uridyltransferase. Of the genes associated with carriage isolates, a gene variant of porB encoding PorB class 3, the genes pilE/pilS and tspB have known functions. The SNP associated with carriage was fkbp D33N, encoding a FK506-binding protein (FKBP). K-mers from PilS, tbpB and tspB were found to be associated with carriage, while k-mers from mtrD and tbpA were associated with invasiveness. In the genes fkbp, glmU, PilC and pilE, k-mers were found that were associated with both carriage and invasive isolates, indicating that specific variations within these genes could play a role in invasiveness. The data presented here highlight genetic traits that are significantly associated with invasive or carriage N. meningitidis across the species population. These traits could prove essential to our understanding of the pathogenicity of N. meningitidis and could help to identify future vaccine targets.

Comprehensive Genome Analysis of Neisseria meningitidis from South America Reveals a Distinctive Pathogenicity-Related Prophage Repertoire

Neisseria meningitidis, a bacterium that colonizes in the human nasopharynx, occasionally causes invasive meningococcal disease leading to meningitis or septicemia. Different serogroups and lineages (clonal complexes) are related to the occurrence and epidemiology of N. meningitidis. Despite vaccines for most serogroups, N. meningitidis lineages causing unusual clinical manifestations and a higher fatality rate compared to other lineages have been reported in South America. The present study focused on exploring the diversity of N. meningitidis prophages from South America and their relationship with the epidemiological variables of these strains. We found a high diversity of prophages among the different clonal complexes. By comparing them with previously described N. meningitidis phages and prophages, we revealed groups of prophages sharing similar compositions, which could be useful for prophage comparison in N. meningitidis. Furthermore, we observed a high correlation between the prophage content and epidemiological features, e.g., pathogenicity or clonal complex. Additionally, a distinctive filamentous prophage named here as IMSAR-11 (Invasive Meningococci from South America Related to cc11) was identified. Interestingly, two versions of IMSAR-11, circular and chromosomally integrated, were found. Overall, this study reinforces the importance of the genomic characterization of circulating N. meningitidis lineages to generate new targets for lineage monitoring, diagnosis, or appropriateness of vaccine development. Further studies are necessary to understand the role of these prophages in the persistence, dispersal, and virulence of N. meningitidis in the world.

Global genomic and proteomic analysis indicates co-evolution of Neisseria species and with their human host

Neisseria, a genus from the beta-proteobacteria class, is of potential clinical importance. This genus contains both pathogenic and commensal strains. Gonorrhea and meningitis are two major diseases caused by pathogens belonging to this genus. With the increased use of antimicrobial agents against these pathogens they have evolved the antimicrobial resistance capacity making these diseases nearly untreatable. The set of anti-bacterial resistance genes (resistome) and genes associated with signal processing (secretomes) are crucial for the host-microbial interaction. With the virtue of whole-genome sequences and computational biology, it is now possible to study the genomic and proteomic riddles of Neisseria along with their comprehensive evolutionary and metabolic profiling. We have studied relative synonymous codon usage, amino acid usage, reverse ecology, comparative genomics, evolutionary analysis and pathogen-host (Neisseria-human) interaction through bioinformatics analysis. Our analysis revealed the co-evolution of Neisseria genomes with the human host. Moreover, the co-occurrence of Neisseria and humans has been supported through reverse ecology analysis. A differential pattern of the evolutionary rate of resistomes and secretomes was evident among the pathogenic and commensal strains. Comparative genomics supported the presence of virulent genes in both pathogenic and commensal strains of the select genus. Our analysis also indicated a transition from commensal to pathogenic Neisseria strains through the long run of evolution.

The Neisseria gonorrhoeae Accessory Genome and Its Association with the Core Genome and Antimicrobial Resistance

The bacterial accessory genome provides the genetic flexibility needed to facilitate environment and host adaptation. In Neisseria gonorrhoeae, known accessory elements include plasmids which can transfer and mediate antimicrobial resistance (AMR); however, chromosomal accessory genes could also play a role in AMR. Here, the gonococcal accessory genome was characterized using gene-by-gene approaches and its association with the core genome and AMR were assessed. The gonococcal accessory gene pool consisted of 247 genes, which were mainly genes located on large mobile genetic elements, phage associated genes, or genes encoding putative secretion systems. Accessory elements showed similar synteny across genomes, indicating either a predisposition for particular genomic locations or ancestral inheritance that are conserved during strain expansion. Significant associations were found between the prevalence of accessory elements and core genome multi-locus sequence types (cgMLST), consistent with a structured gonococcal population despite frequent horizontal gene transfer (HGT). Increased prevalence of putative DNA exchange regulators was significantly associated with AMR, which included a putative secretion system, methyltransferases and a toxin-antitoxin system. Although frequent HGT results in high genetic diversity in the gonococcus, we found that this is mediated by a small gene pool. In fact, a highly organized genome composition was identified with a strong association between the accessory and core genome. Increased prevalence of DNA exchange regulators in antimicrobial resistant isolates suggests that genetic material exchange plays a role in the development or maintenance of AMR. These findings enhance our understanding of gonococcal genome architecture and have important implications for gonococcal population biology. IMPORTANCE The emergence of antimicrobial resistance (AMR) against third generation cephalosporins in Neisseria gonorrhoeae is a major public health concern, as these are antibiotics of last resort for the effective treatment of gonorrhea. Although the resistance mechanisms against this class of antibiotics have not been entirely resolved, resistance against other classes of antibiotics, such as tetracyclines, is known to be mediated through plasmids, which are known gonococcal extra-chromosomal accessory elements. A complete assessment of the chromosomal accessory genome content and its role in AMR has not yet been undertaken. Here, we comprehensively characterize the gonococcal accessory genome to better understand genome architecture as well as the evolution and mechanisms of AMR in this species.

Prediction of Prophages and Their Host Ranges in Pathogenic and Commensal Neisseria Species

The genus Neisseria includes two pathogenic species, N. gonorrhoeae and N. meningitidis, and numerous commensal species. Neisseria species frequently exchange DNA with one another, primarily via transformation and homologous recombination and via multiple types of mobile genetic elements (MGEs). Few Neisseria bacteriophages (phages) have been identified, and their impact on bacterial physiology is poorly understood. Furthermore, little is known about the range of species that Neisseria phages can infect. In this study, we used three virus prediction tools to scan 248 genomes of 21 different Neisseria species and identified 1,302 unique predicted prophages. Using comparative genomics, we found that many predictions are dissimilar from prophages and other MGEs previously described to infect Neisseria species. We also identified similar predicted prophages in genomes of different Neisseria species. Additionally, we examined CRISPR-Cas targeting of each Neisseria genome and predicted prophage. While CRISPR targeting of chromosomal DNA appears to be common among several Neisseria species, we found that 20% of the prophages we predicted are targeted significantly more than the rest of the bacterial genome in which they were identified (i.e., backbone). Furthermore, many predicted prophages are targeted by CRISPR spacers encoded by other species. We then used these results to infer additional host species of known Neisseria prophages and predictions that are highly targeted relative to the backbone. Together, our results suggest that we have identified novel Neisseria prophages, several of which may infect multiple Neisseria species. These findings have important implications for understanding horizontal gene transfer between members of this genus. IMPORTANCE Drug-resistant Neisseria gonorrhoeae is a major threat to human health. Commensal Neisseria species are thought to serve as reservoirs of antibiotic resistance and virulence genes for the pathogenic species N. gonorrhoeae and N. meningitidis. Therefore, it is important to understand both the diversity of mobile genetic elements (MGEs) that can mediate horizontal gene transfer within this genus and the breadth of species these MGEs can infect. In particular, few bacteriophages (phages) are known to infect Neisseria species. In this study, we identified a large number of candidate phages integrated in the genomes of commensal and pathogenic Neisseria species, many of which appear to be novel phages. Importantly, we discovered extensive interspecies targeting of predicted phages by Neisseria CRISPR-Cas systems, which may reflect their movement between different species. Uncovering the diversity and host range of phages is essential for understanding how they influence the evolution of their microbial hosts.

Neisseria genes required for persistence identified via in vivo screening of a transposon mutant library

The mechanisms used by human adapted commensal Neisseria to shape and maintain a niche in their host are poorly defined. These organisms are common members of the mucosal microbiota and share many putative host interaction factors with Neisseria meningitidis and Neisseria gonorrhoeae. Evaluating the role of these shared factors during host carriage may provide insight into bacterial mechanisms driving both commensalism and asymptomatic infection across the genus. We identified host interaction factors required for niche development and maintenance through in vivo screening of a transposon mutant library of Neisseria musculi, a commensal of wild-caught mice which persistently and asymptomatically colonizes the oral cavity and gut of CAST/EiJ and A/J mice. Approximately 500 candidate genes involved in long-term host interaction were identified. These included homologs of putative N. meningitidis and N. gonorrhoeae virulence factors which have been shown to modulate host interactions in vitro. Importantly, many candidate genes have no assigned function, illustrating how much remains to be learned about Neisseria persistence. Many genes of unknown function are conserved in human adapted Neisseria species; they are likely to provide a gateway for understanding the mechanisms allowing pathogenic and commensal Neisseria to establish and maintain a niche in their natural hosts. Validation of a subset of candidate genes confirmed a role for a polysaccharide capsule in N. musculi persistence but not colonization. Our findings highlight the potential utility of the Neisseria musculi-mouse model as a tool for studying the pathogenic Neisseria; our work represents a first step towards the identification of novel host interaction factors conserved across the genus.

The Neisseria genus contains many genetically related commensals of animals and humans, and two human pathogens, Neisseria gonorrhoeae and Neisseria meningitidis. The mechanisms allowing commensal Neisseria to maintain a niche in their host is little understood. To identify genes required for persistence, we screened a library of transposon mutants of Neisseria musculi, a commensal of wild-caught mice, in CAST/EiJ mice, which persistently and asymptomatically colonizes. Approximately 500 candidate host interaction genes were identified. A subset of these are homologs of N. meningitidis and N. gonorrhoeae genes known to modulate pathogen-host interactions in vitro. Many candidate genes have no known function, demonstrating how much remains to be learned about N. musculi niche maintenance. As many genes of unknown function are conserved in human adapted Neisseria, they provide a gateway for understanding Neisseria persistence mechanisms in general.

Modelling evolutionary pathways for commensalism and hypervirulence in Neisseria meningitidis

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.

Bioinformatics tools used for whole-genome sequencing analysis of Neisseria gonorrhoeae: a literature review

Whole-genome sequencing (WGS) data are well established for the investigation of gonococcal transmission, antimicrobial resistance prediction, population structure determination and population dynamics. A variety of bioinformatics tools, repositories, services and platforms have been applied to manage and analyze Neisseria gonorrhoeae WGS datasets. This review provides an overview of the various bioinformatics approaches and resources used in 105 published studies (as of 30 April 2021). The challenges in the analysis of N. gonorrhoeae WGS datasets, as well as future bioinformatics requirements, are also discussed.

Simulation study and comparative evaluation of viral contiguous sequence identification tools

BACKGROUND

Viruses, including bacteriophages, are important components of environmental and human associated microbial communities. Viruses can act as extracellular reservoirs of bacterial genes, can mediate microbiome dynamics, and can influence the virulence of clinical pathogens. Various targeted metagenomic analysis techniques detect viral sequences, but these methods often exclude large and genome integrated viruses. In this study, we evaluate and compare the ability of nine state-of-the-art bioinformatic tools, including Vibrant, VirSorter, VirSorter2, VirFinder, DeepVirFinder, MetaPhinder, Kraken 2, Phybrid, and a BLAST search using identified proteins from the Earth Virome Pipeline to identify viral contiguous sequences (contigs) across simulated metagenomes with different read distributions, taxonomic compositions, and complexities.

RESULTS

Of the tools tested in this study, VirSorter achieved the best F1 score while Vibrant had the highest average F1 score at predicting integrated prophages. Though less balanced in its precision and recall, Kraken2 had the highest average precision by a substantial margin. We introduced the machine learning tool, Phybrid, which demonstrated an improvement in average F1 score over tools such as MetaPhinder. The tool utilizes machine learning with both gene content and nucleotide features. The addition of nucleotide features improves the precision and recall compared to the gene content features alone.Viral identification by all tools was not impacted by underlying read distribution but did improve with contig length. Tool performance was inversely related to taxonomic complexity and varied by the phage host. For instance, Rhizobium and Enterococcus phages were identified consistently by the tools; whereas, Neisseria prophage sequences were commonly missed in this study.

CONCLUSION

This study benchmarked the performance of nine state-of-the-art bioinformatic tools to identify viral contigs across different simulation conditions. This study explored the ability of the tools to identify integrated prophage elements traditionally excluded from targeted sequencing approaches. Our comprehensive analysis of viral identification tools to assess their performance in a variety of situations provides valuable insights to viral researchers looking to mine viral elements from publicly available metagenomic data.