Whole-genome sequencing of non-typeable Haemophilus influenzae isolated from a tertiary care hospital in Surabaya, Indonesia.

BACKGROUND: Haemophilus influenzae causes life-threatening invasive diseases such as septicaemia and meningitis. Reports on circulating H. influenzae causing invasive disease in lower-middle income settings, including Indonesia, are lacking. This study describes the serotype distributions and whole-genome sequence (WGS) data of H. influenzae isolated from hospitalized patients at Soetomo Hospital, Surabaya, Indonesia. METHODS: H. influenzae isolates were isolated from blood and pleural fluid specimens and identified using culture-based and molecular methods, followed by serotyping and WGS using RT‒PCR and Illumina MiSeq, respectively. Sequencing reads were assembled, and further analyses were undertaken to determine the genomic content and reconstruct the phylogeny. A second dataset consisting of publicly available H. influenzae genomes was curated to conduct phylogenetic analyses of isolates in this study in the context of globally circulating isolates. RESULTS: Ten H. influenzae isolates from hospitalized patients were collected, and septicaemia was the most common diagnosis (n=8). RT‒PCR and WGS were performed to determine whether all the isolates were nontypeable H. influenzae (NTHi). There were four newly identified STs distributed across the two main clusters. A total of 91 out of 126 virulence factor (VF)-related genes in Haemophilus sp. were detected in at least one isolate. Further evaluation incorporating a global collection of H. influenzae genomes confirmed the diverse population structure of NTHi in this study. CONCLUSION: This study showed that all H. influenzae recovered from invasive disease patients were nonvaccine-preventable NTHi isolates. The inclusion of WGS revealed four novel STs and the possession of key VF-associated genes.

Genome-wide association studies reveal the role of polymorphisms affecting factor H binding protein expression in host invasion by Neisseria meningitidis.

Many invasive bacterial diseases are caused by organisms that are ordinarily harmless components of the human microbiome. Effective interventions against these microbes require an understanding of the processes whereby symbiotic or commensal relationships transition into pathology. Here, we describe bacterial genome-wide association studies (GWAS) of Neisseria meningitidis, a common commensal of the human respiratory tract that is nevertheless a leading cause of meningitis and sepsis. An initial GWAS discovered bacterial genetic variants, including single nucleotide polymorphisms (SNPs), associated with invasive meningococcal disease (IMD) versus carriage in several loci across the meningococcal genome, encoding antigens and other extracellular components, confirming the polygenic nature of the invasive phenotype. In particular, there was a significant peak of association around the fHbp locus, encoding factor H binding protein (fHbp), which promotes bacterial immune evasion of human complement by recruiting complement factor H (CFH) to the meningococcal surface. The association around fHbp with IMD was confirmed by a validation GWAS, and we found that the SNPs identified in the validation affected the 5' region of fHbp mRNA, altering secondary RNA structures, thereby increasing fHbp expression and enhancing bacterial escape from complement-mediated killing. This finding is consistent with the known link between complement deficiencies and CFH variation with human susceptibility to IMD. These observations demonstrate the importance of human and bacterial genetic variation across the fHbp:CFH interface in determining IMD susceptibility, the transition from carriage to disease.

Genome-Wide Association Studies Identify an Association of Transferrin Binding Protein B Variation and Invasive Serogroup Y Meningococcal Disease in Older Adults.

BACKGROUND: Neisseria meningitidis serogroup Y, especially ST-23 clonal complex (Y:cc23), represents a larger proportion of invasive meningococcal disease (IMD) in older adults compared to younger individuals. This study explored the meningococcal genetic variation underlying this association. METHODS: Maximum-likelihood phylogenies and the pangenome were analyzed using whole-genome sequence (WGS) data from 200 Y:cc23 isolates in the Neisseria PubMLST database. Genome-wide association studies (GWAS) were performed on WGS data from 250 Y:cc23 isolates from individuals with IMD aged ≥65 years versus < 65 years. RESULTS: Y:cc23 meningococcal variants did not cluster by age group or disease phenotype in phylogenetic analyses. Pangenome comparisons found no differences in presence or absence of genes in IMD isolates from the different age groups. GWAS identified differences in nucleotide polymorphisms within the transferrin-binding protein B (tbpB) gene in isolates from individuals ≥65 years of age. TbpB structure modelling suggests these may impact binding of human transferrin. CONCLUSIONS: These data suggest differential iron scavenging capacity amongst Y:cc23 meningococci isolated from older compared to younger patients. Iron acquisition is essential for many bacterial pathogens including the meningococcus. These polymorphisms may facilitate colonization, thereby increasing the risk of disease in vulnerable older people with altered nasopharyngeal microbiomes and nutritional status.

Genetic determinants of genus-level glycan diversity in a bacterial protein glycosylation system.

The human pathogens N. gonorrhoeae and N. meningitidis display robust intra- and interstrain glycan diversity associated with their O-linked protein glycosylation (pgl) systems. In an effort to better understand the evolution and function of protein glycosylation operating there, we aimed to determine if other human-restricted, Neisseria species similarly glycosylate proteins and if so, to assess the levels of glycoform diversity. Comparative genomics revealed the conservation of a subset of genes minimally required for O-linked protein glycosylation glycan and established those pgl genes as core genome constituents of the genus. In conjunction with mass spectrometric-based glycan phenotyping, we found that extant glycoform repertoires in N. gonorrhoeae, N. meningitidis and the closely related species N. polysaccharea and N. lactamica reflect the functional replacement of a progenitor glycan biosynthetic pathway. This replacement involved loss of pgl gene components of the primordial pathway coincident with the acquisition of two exogenous glycosyltransferase genes. Critical to this discovery was the identification of a ubiquitous but previously unrecognized glycosyltransferase gene (pglP) that has uniquely undergone parallel but independent pseudogenization in N. gonorrhoeae and N. meningitidis. We suggest that the pseudogenization events are driven by processes of compositional epistasis leading to gene decay. Additionally, we documented instances where inter-species recombination influences pgl gene status and creates discordant genetic interactions due ostensibly to the multi-locus nature of pgl gene networks. In summary, these findings provide a novel perspective on the evolution of protein glycosylation systems and identify phylogenetically informative, genetic differences associated with Neisseria species.

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

Expansion of quinolone-resistant Neisseria meningitidis clone ChinaCC4821-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 ChinaCC4821-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.

Association of Neisseria gonorrhoeae Plasmids With Distinct Lineages and The Economic Status of Their Country of Origin.

Plasmids are vehicles for horizontal gene transfer between bacteria, and in Neisseria gonorrhoeae plasmids can mediate high-level antimicrobial resistance (AMR). Using genomic and phylogenetic analyses, we show that plasmids are widespread in a collection of 3724 gonococcal isolates from 56 countries, and characterized the conjugative, ß-lactamase and cryptic plasmids. We found that variants of the conjugative plasmid (which can mediate tetracycline resistance) and the ß-lactamase plasmid expressing TEM-135 are associated with distinct gonococcal lineages. Furthermore, AMR plasmids are significantly more prevalent in gonococci from less wealthy countries, highlighting the need for further studies. More than 94% of gonococci possess the cryptic plasmid, with its absence correlated with the presence of a novel chromosomal type IV secretion system. Our results reveal the extent of plasmid-mediated AMR in the gonococcus, particularly in less wealthy countries, where diagnostic and therapeutic options can be limited, and highlight the risk of their global spread.

Neisseria gonorrhoeae Population Genomics: Use of the Gonococcal Core Genome to Improve Surveillance of Antimicrobial Resistance.

BACKGROUND: Gonorrhea, caused by the bacterium Neisseria gonorrhoeae, is a globally prevalent sexually transmitted infection. The dynamics of gonococcal population biology have been poorly defined due to a lack of resolution in strain typing methods. METHODS: In this study, we assess how the core genome can be used to improve our understanding of gonococcal population structure compared with current typing schemes. RESULTS: A total of 1668 loci were identified as core to the gonococcal genome. These were organized into a core genome multilocus sequence typing scheme (N gonorrhoeae cgMLST v1.0). A clustering algorithm using a threshold of 400 allelic differences between isolates resolved gonococci into discrete and stable core genome groups, some of which persisted for multiple decades. These groups were associated with antimicrobial genotypes and non-overlapping NG-STAR and NG-MAST sequence types. The MLST-STs were more widely distributed among core genome groups. CONCLUSIONS: Clustering with cgMLST identified globally distributed, persistent, gonococcal lineages improving understanding of the population biology of gonococci and revealing its population structure. These findings have implications for the emergence of antimicrobial resistance in gonococci and how this is associated with lineages, some of which are more predisposed to developing antimicrobial resistance than others.

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.

Identification of Novel Neisseria gonorrhoeae Lineages Harboring Resistance Plasmids in Coastal Kenya.

Background: Africa has the highest incidence of gonorrhea in the world. However, little is known about gonococcal populations in this continent or mechanisms of antimicrobial resistance (AMR). Methods: Whole-genome sequence data were analyzed from 103 Neisseria gonorrhoeae isolates from 73 patients, mainly men who have sex with men, from coastal Kenya. We annotated loci, defined the core genome, defined mechanisms of AMR, and performed phylogenetic analysis. For patients with multiple episodes of gonorrhea, we determined whether infections occurred with related strains. Results: We identified 3 clusters of isolates that are phylogenetically distinct from isolates found elsewhere. Plasmids were virtually ubiquitous: pTetM and pblaTEM were found in 97%, and 55% of isolates, respectively. This was associated with high doxycycline use for undiagnosed sexually transmitted infections. Twenty-three percent of multiple episodes of gonorrhea in the same individual were caused by a related strain, suggesting inadequate treatment or reinfection. Conclusions: The prevalence of plasmid-mediated AMR in Kenyan gonococci contrasts with that in wealthy countries, where AMR is largely chromosomally mediated. Antimicrobials have a profound effect on the maintenance of lineages harboring plasmids. Doxycycline can select for tetracycline and penicillin resistance, through plasmid cooperation. Understanding the mechanisms of AMR in high-risk groups is required to inform treatment strategies.

Genomic analyses of the Chlamydia trachomatis core genome show an association between chromosomal genome, plasmid type and disease.

BACKGROUND: Chlamydia trachomatis (Ct) plasmid has been shown to encode genes essential for infection. We evaluated the population structure of Ct using whole-genome sequence data (WGS). In particular, the relationship between the Ct genome, plasmid and disease was investigated. RESULTS: WGS data from 157 Ct isolates deposited in the Chlamydiales pubMLST database ( http://pubMLST.org/chlamydiales/ ) were annotated with 902 genes including the core and accessory genome. Plasmid associated genes were annotated and a plasmid MLST scheme was defined allowing plasmid sequence types to be determined. Plasmid allelic variation was investigated. Phylogenetic relationships were examined using the Genome Comparator tool available in pubMLST. Phylogenetic analyses identified four distinct Ct core genome clusters and six plasmid clusters, with a strong association between the chromosomal genotype and plasmid. This in turn was linked to ompA genovars and disease phenotype. Horizontal genetic transfer of plasmids was observed for three urogenital-associated isolates, which possessed plasmids more commonly found in isolates resulting from ocular infections. The pgp3 gene was identified as the most polymorphic plasmid gene and pgp4 was the most conserved. CONCLUSION: A strong association between chromosomal genome, plasmid type and disease was observed, consistent with previous studies. This suggests co-evolution of the Ct chromosome and their plasmids, but we confirmed that plasmid transfer can occur between isolates. These data provide a better understanding of the genetic diversity occurring across the Ct genome in association with the plasmid content.

Hierarchical genomic analysis of carried and invasive serogroup A Neisseria meningitidis during the 2011 epidemic in Chad.

BACKGROUND: Serogroup A Neisseria meningitidis (NmA) was the cause of the 2011 meningitis epidemics in Chad. This bacterium, often carried asymptomatically, is considered to be an "accidental pathogen"; however, the transition from carriage to disease phenotype remains poorly understood. This study examined the role genetic diversity might play in this transition by comparing genomes from geographically and temporally matched invasive and carried NmA isolates. RESULTS: All 23 NmA isolates belonged to the ST-5 clonal complex (cc5). Ribosomal MLST comparison with other publically available NmA:cc5 showed that isolates were closely related, although those from Chad formed two distinct branches and did not cluster with other NmA, based on their MLST profile, geographical and temporal location. Whole genome MLST (wgMLST) comparison identified 242 variable genes among all Chadian isolates and clustered them into three distinct phylogenetic groups (Clusters 1, 2, and 3): no systematic clustering by disease or carriage source was observed. There was a significant difference (p = 0.0070) between the mean age of the individuals from which isolates from Cluster 1 and Cluster 2 were obtained, irrespective of whether the person was a case or a carrier. CONCLUSIONS: Whole genome sequencing provided high-resolution characterization of the genetic diversity of these closely related NmA isolates. The invasive meningococcal isolates obtained during the epidemic were not homogeneous; rather, a variety of closely related but distinct clones were circulating in the human population with some clones preferentially colonizing specific age groups, reflecting a potential age-related niche adaptation. Systematic genetic differences were not identified between carriage and disease isolates consistent with invasive meningococcal disease being a multi-factorial event resulting from changes in host-pathogen interactions along with the bacterium.

Molecular characterization of invasive capsule null Neisseria meningitidis in South Africa.

BACKGROUND: The meningococcal capsule is an important virulence determinant. Unencapsulated meningococci lacking capsule biosynthesis genes and containing the capsule null locus (cnl) are predominantly non-pathogenic. Rare cases of invasive meningococcal disease caused by cnl isolates belonging to sequence types (ST) and clonal complexes (cc) ST-845 (cc845), ST-198 (cc198), ST-192 (cc192) and ST-53 (cc53) have been documented. The clinical significance of these isolates however remains unclear. We identified four invasive cnl meningococci through laboratory-based surveillance in South Africa from 2003 through 2013, which we aimed to characterize using whole genome data. RESULTS: One isolate [NG: P1.7-2,30: F1-2: ST-53 (cc53)] contained cnl allele 12, and caused empyema in an adult male with bronchiectasis from tuberculosis, diabetes mellitus and a smoking history. Three isolates were NG: P1.18-11,42-2: FΔ: ST-192 (cc192) and contained cnl allele 2. One patient was an adolescent male with meningitis. The remaining two isolates were from recurrent disease episodes (8 months apart) in a male child with deficiency of the sixth complement component, and with the exception of two single nucleotide polymorphisms, contained identical core genomes. The ST-53 (cc53) isolate possessed alleles for NHBA peptide 191 and fHbp variant 2; whilst the ST-192 (cc192) isolates contained NHBA peptide 704 and fHbp variant 3. All four isolates lacked nadA. Comparison of the South African genomes to 61 additional cnl genomes on the PubMLST Neisseria database ( http://pubmlst.org/neisseria/ ), determined that most putative virulence genes could be found in both invasive and carriage phenotypes. CONCLUSIONS: Although rare, invasive disease by cnl meningococci may be associated with host immunodeficiency and such patients may benefit from protein-based meningococcal vaccines.

Genomic analysis of urogenital and rectal Neisseria meningitidis isolates reveals encapsulated hyperinvasive meningococci and coincident multidrug-resistant gonococci.

OBJECTIVE: Invasive meningococcal disease (IMD) outbreaks in men who have sex with men (MSM) have been associated with meningococcal colonisation of the urethra and rectum, but little is known about this colonisation or co-colonisation with the closely related gonococcus. Whole genome sequencing (WGS) was employed to explore these phenomena. METHODS: Meningococci isolated from the urogenital tract and rectum (n=23) and coincident gonococci (n=14) were analysed by WGS along with contemporary meningococci from IMD (n=11). All isolates were obtained from hospital admissions in Brighton, UK, 2011-2013. Assembled WGS were deposited in the PubMLST/neisseria database (http://pubmlst.org/neisseria) and compared at genomic loci common to gonococci or meningococci. RESULTS: As expected, most meningococci from IMD were encapsulated and belonged to hyperinvasive lineages. So too were meningococci found in the urogenital tract and rectum, contrasting to those asymptomatically carried in the nasopharynx where such meningococci are rare. Five hyperinvasive meningococcal lineages and four distinct gonococcal genotypes were recovered, including multiresistant ST-1901 (NG MAST-1407) gonococci. CONCLUSIONS: These data were consistent with a predisposition for potentially virulent encapsulated hyperinvasive meningococci to colonise the urethra and rectum, which suggests their involvement in MSM IMD outbreaks. The coincidence of multiresistant gonococci raises wider public health concerns.

Genomic Analysis of Serogroup Y Neisseria meningitidis Isolates Reveals Extensive Similarities Between Carriage-Associated and Disease-Associated Organisms.

BACKGROUND: Neisseria meningitidis is a frequent colonizer of the human nasopharynx, with asymptomatic carriage providing the reservoir for invasive, disease-causing strains. Serogroup Y (MenY) strains are a major cause of meningococcal disease. High-resolution genetic analyses of carriage and disease isolates can establish epidemiological relationships and identify potential virulence factors. METHODS: Whole-genome sequence data were obtained for 99 MenY carriage isolates recovered in the United Kingdom during 1997-2010. Sequences were compared to those of 73 MenY invasive isolates recovered during 2010-2011, using a gene-by-gene approach. RESULTS: Comparisons across 1605 core genes resolved 91% of isolates into one of 8 clusters containing closely related disease and carriage isolates. Six clusters contained carried meningococci isolated during 1997-2001, suggesting temporal stability. One cluster of isolates, predominately sharing the designation Y: P1.5-1,10-1: F4-1: ST-1655 (cc23), was resolved into one subcluster with 86% carriage isolates and a second with 90% invasive isolates. These subclusters were defined by specific allelic differences in 5 core genes encoding glycerate kinase (glxK), valine-pyruvate transaminase (avtA), superoxide dismutase (sodB), and 2 hypothetical proteins. CONCLUSIONS: High-resolution genetic analyses detected long-term temporal stability and temporally overlapping carriage and disease populations for MenY clones but also evidence of a disease-associated clone.

Population and Functional Genomics of Neisseria Revealed with Gene-by-Gene Approaches.

Rapid low-cost whole-genome sequencing (WGS) is revolutionizing microbiology; however, complementary advances in accessible, reproducible, and rapid analysis techniques are required to realize the potential of these data. Here, investigations of the genus Neisseria illustrated the gene-by-gene conceptual approach to the organization and analysis of WGS data. Using the gene and its link to phenotype as a starting point, the BIGSdb database, which powers the PubMLST databases, enables the assembly of large open-access collections of annotated genomes that provide insight into the evolution of the Neisseria, the epidemiology of meningococcal and gonococcal disease, and mechanisms of Neisseria pathogenicity.

Isolation and characterization of Neisseria musculi sp. nov., from the wild house mouse.

Members of the genus Neisseria have been isolated from or detected in a wide range of animals, from non-human primates and felids to a rodent, the guinea pig. By means of selective culture, biochemical testing, Gram staining and PCR screening for the Neisseria-specific internal transcribed spacer region of the rRNA operon, we isolated four strains of the genus Neisseria from the oral cavity of the wild house mouse, Mus musculus subsp. domesticus. The isolates are highly related and form a separate clade in the genus, as judged by tree analyses using either multi-locus sequence typing of ribosomal genes or core genes. One isolate, provisionally named Neisseria musculi sp. nov. (type strain AP2031T=DSM 101846T=CCUG 68283T=LMG 29261T), was studied further. Strain AP2031T/N. musculi grew well in vitro. It was naturally competent, taking up DNA in a DNA uptake sequence and pilT-dependent manner, and was amenable to genetic manipulation. These and other genomic attributes of N. musculi sp. nov. make it an ideal candidate for use in developing a mouse model for studying Neisseria-host interactions.

Bacteriocin-mediated competition in cystic fibrosis lung infections.

Bacteriocins are toxins produced by bacteria to kill competitors of the same species. Theory and laboratory experiments suggest that bacteriocin production and immunity play a key role in the competitive dynamics of bacterial strains. The extent to which this is the case in natural populations,especially human pathogens, remains to be tested. We examined the role of bacteriocins in competition using Pseudomonas aeruginosa strains infecting lungs of humans with cystic fibrosis (CF). We assessed the ability of different strains to kill each other using phenotypic assays, and sequenced their genomes to determine what bacteriocins (pyocins) they carry. We found that(i) isolates from later infection stages inhibited earlier infecting strains less,but were more inhibited by pyocins produced by earlier infecting strains and carried fewer pyocin types; (ii) this difference between early and late infections appears to be caused by a difference in pyocin diversity between competing genotypes and not by loss of pyocin genes within a lineage overtime; (iii) pyocin inhibition does not explain why certain strains outcompete others within lung infections; (iv) strains frequently carry the pyocin-killing gene, but not the immunity gene, suggesting resistance occurs via other unknown mechanisms. Our results show that, in contrast to patterns observed in experimental studies, pyocin production does not appear to have a major influence on strain competition during CF lung infections.

Identifying Neisseria species by use of the 50S ribosomal protein L6 (rplF) gene.

The comparison of 16S rRNA gene sequences is widely used to differentiate bacteria; however, this gene can lack resolution among closely related but distinct members of the same genus. This is a problem in clinical situations in those genera, such as Neisseria, where some species are associated with disease while others are not. Here, we identified and validated an alternative genetic target common to all Neisseria species which can be readily sequenced to provide an assay that rapidly and accurately discriminates among members of the genus. Ribosomal multilocus sequence typing (rMLST) using ribosomal protein genes has been shown to unambiguously identify these bacteria. The PubMLST Neisseria database (http://pubmlst.org/neisseria/) was queried to extract the 53 ribosomal protein gene sequences from 44 genomes from diverse species. Phylogenies reconstructed from these genes were examined, and a single 413-bp fragment of the 50S ribosomal protein L6 (rplF) gene was identified which produced a phylogeny that was congruent with the phylogeny reconstructed from concatenated ribosomal protein genes. Primers that enabled the amplification and direct sequencing of the rplF gene fragment were designed to validate the assay in vitro and in silico. Allele sequences were defined for the gene fragment, associated with particular species names, and stored on the PubMLST Neisseria database, providing a curated electronic resource. This approach provides an alternative to 16S rRNA gene sequencing, which can be readily replicated for other organisms for which more resolution is required, and it has potential applications in high-resolution metagenomic studies.

Characterizing the diversity and commensal origins of penA mosaicism in the genus Neisseria.

Mosaic penA alleles formed through horizontal gene transfer (HGT) have been instrumental to the rising incidence of ceftriaxone-resistant gonococcal infections. Although interspecies HGT of regions of the penA gene between Neisseria gonorrhoeae and commensal Neisseria species has been described, knowledge concerning which species are the most common contributors to mosaic penA alleles is limited, with most studies examining only a small number of alleles. Here, we investigated the origins of recombinant penA alleles through in silico analyses that incorporated 1700 penA alleles from 35 513 Neisseria isolates, comprising 15 different Neisseria species. We identified Neisseria subflava and Neisseria cinerea as the most common source of recombinant sequences in N. gonorrhoeae penA. This contrasted with Neisseria meningitidis penA, for which the primary source of recombinant DNA was other meningococci, followed by Neisseria lactamica. Additionally, we described the distribution of polymorphisms implicated in antimicrobial resistance in penA, and found that these are present across the genus. These results provide insight into resistance-related changes in the penA gene across human-associated Neisseria species, illustrating the importance of genomic surveillance of not only the pathogenic Neisseria, but also of the oral niche-associated commensals from which these pathogens are sourcing key genetic variation.

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae.

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.