Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

Genetic characterization and estimated 4CMenB vaccine strain coverage of 284 Neisseria meningitidis isolates causing invasive meningococcal disease in Argentina in 2010-2014

Meningococcal (Neisseria meningitidis) serogroup B (MenB) strain antigens are diverse and a limited number of strains can be evaluated using the human serum bactericidal antibody (hSBA) assay. The genetic Meningococcal Antigen Typing System (gMATS) was developed to predict the likelihood of coverage for large numbers of isolates by the 4CMenB vaccine, which includes antigens Neisseria adhesin A (NadA), Neisserial Heparin-Binding Antigen (NHBA), factor H-binding protein (fHbp), and Porin A (PorA). In this study, we characterized by whole-genome analyses 284 invasive MenB isolates collected from 2010 to 2014 by the Argentinian National Laboratories Network (52-61 isolates per year). Strain coverage was estimated by gMATS on all isolates and by hSBA assay on 74 randomly selected isolates, representative of the whole panel. The four most common clonal complexes (CCs), accounting for 81.3% of isolates, were CC-865 (75 isolates, 26.4%), CC-32 (59, 20.8%), CC-35 (59, 20.8%), and CC-41/44 (38, 13.4%). Vaccine antigen genotyping showed diversity. The most prevalent variants/peptides were fHbp variant 2, NHBA peptides 24, 21, and 2, and PorA variable region 2 profiles 16-36 and 14. The nadA gene was present in 66 (23.2%) isolates. Estimated strain coverage by hSBA assay showed 78.4% of isolates were killed by pooled adolescent sera, and 51.4% and 64.9% (based on two different thresholds) were killed by pooled infant sera. Estimated coverage by gMATS (61.3%; prediction interval: 55.5%, 66.7%) was consistent with the infant hSBA assay results. Continued genomic surveillance is needed to evaluate the persistence of major MenB CCs in Argentina.

Sweet complexity: O-linked protein glycosylation in pathogenic Neisseria

The genus Neisseria, which colonizes mucosal surfaces, includes both commensal and pathogenic species that are exclusive to humans. The two pathogenic Neisseria species are closely related but cause quite different diseases, meningococcal sepsis and meningitis (Neisseria meningitidis) and sexually transmitted gonorrhea (Neisseria gonorrhoeae). Although obvious differences in bacterial niches and mechanisms for transmission exists, pathogenic Neisseria have high levels of conservation at the levels of nucleotide sequences, gene content and synteny. Species of Neisseria express broad-spectrum O-linked protein glycosylation where the glycoproteins are largely transmembrane proteins or lipoproteins localized on the cell surface or in the periplasm. There are diverse functions among the identified glycoproteins, for example type IV biogenesis proteins, proteins involved in antimicrobial resistance, as well as surface proteins that have been suggested as vaccine candidates. The most abundant glycoprotein, PilE, is the major subunit of pili which are an important colonization factor. The glycans attached can vary extensively due to phase variation of protein glycosylation (pgl) genes and polymorphic pgl gene content. The exact roles of glycosylation in Neisseria remains to be determined, but increasing evidence suggests that glycan variability can be a strategy to evade the human immune system. In addition, pathogenic and commensal Neisseria appear to have significant glycosylation differences. Here, the current knowledge and implications of protein glycosylation genes, glycan diversity, glycoproteins and immunogenicity in pathogenic Neisseria are summarized and discussed.

The Evolutionary History of a DNA Methylase Reveals Frequent Horizontal Transfer and Within-Gene Recombination

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages, we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. Additionally, we found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long-distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

Meningococcal disease in North America: Updates from the Global Meningococcal Initiative

This review summarizes the recent Global Meningococcal Initiative (GMI) regional meeting, which explored meningococcal disease in North America. Invasive meningococcal disease (IMD) cases are documented through both passive and active surveillance networks. IMD appears to be decreasing in many areas, such as the Dominican Republic (2016: 18 cases; 2021: 2 cases) and Panama (2008: 1 case/100,000; 2021: <0.1 cases/100,000); however, there is notable regional and temporal variation. Outbreaks persist in at-risk subpopulations, such as people experiencing homelessness in the US and migrants in Mexico. The recent emergence of β-lactamase-positive and ciprofloxacin-resistant meningococci in the US is a major concern. While vaccination practices vary across North America, vaccine uptake remains relatively high. Monovalent and multivalent conjugate vaccines (which many countries in North America primarily use) can provide herd protection. However, there is no evidence that group B vaccines reduce meningococcal carriage. The coronavirus pandemic illustrates that following public health crises, enhanced surveillance of disease epidemiology and catch-up vaccine schedules is key. Whole genome sequencing is a key epidemiological tool for identifying IMD strain emergence and the evaluation of vaccine strain coverage. The Global Roadmap on Defeating Meningitis by 2030 remains a focus of the GMI.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Transmission Dynamics and Microevolution of Neisseria meningitidis During Carriage and Invasive Disease in High School Students in Georgia and Maryland, 2006-2007

BACKGROUND

The mechanisms by which Neisseria meningitidis cause persistent human carriage and transition from carriage to invasive disease have not been fully elucidated.

METHODS

Georgia and Maryland high school students were sampled for pharyngeal carriage of N. meningitidis during the 2006-2007 school year. A total of 321 isolates from 188 carriers and all 67 invasive disease isolates collected during the same time and from the same geographic region underwent whole-genome sequencing. Core-genome multilocus sequence typing was used to compare allelic profiles, and direct read mapping was used to study strain evolution.

RESULTS

Among 188 N. meningitidis culture-positive students, 98 (52.1%) were N. meningitidis culture positive at 2 or 3 samplings. Most students who were positive at >1 sampling (98%) had persistence of a single strain. More than a third of students carried isolates that were highly genetically related to isolates from other students in the same school, and occasional transmission within the same county was also evident. The major pilin subunit gene, pilE, was the most variable gene, and no carrier had identical pilE sequences at different time points.

CONCLUSION

We found strong evidence of local meningococcal transmission at both the school and county levels. Allelic variation within genes encoding bacterial surface structures, particularly pilE, was common.

Diversification in immunogenicity genes caused by selective pressures in invasive meningococci

We studied population genomics of 486 Neisseria meningitidis isolates causing meningitis in the Netherlands during the period 1979-2003 and 2006-2013 using whole-genome sequencing to evaluate the impact of a hyperendemic period of serogroup B invasive disease. The majority of serogroup B isolates belonged to ST-41/44 (41 %) and ST-32 complex (16 %). Comparing the time periods, before and after the decline of serogroup B invasive disease, there was a decrease of ST-41/44 complex sequences (P=0.002). We observed the expansion of a sub-lineage within ST-41/44 complex sequences being associated with isolation from the 1979-2003 time period (P=0.014). Isolates belonging to this sub-lineage expansion within ST-41/44 complex were marked by four antigen allele variants. Presence of these allele variants was associated with isolation from the 1979-2003 time period after correction for multiple testing (Wald test, P=0.0043 for FetA 1-5; P=0.0035 for FHbp 14; P=0.012 for PorA 7-2.4 and P=0.0031 for NHBA two peptide allele). These sequences were associated with 4CMenB vaccine coverage (Fisher's exact test, P<0.001). Outside of the sub-lineage expansion, isolates with markedly lower levels of predicted vaccine coverage clustered in phylogenetic groups showing a trend towards isolation in the 2006-2013 time period (P=0.08). In conclusion, we show the emergence and decline of a sub-lineage expansion within ST-41/44 complex isolates concurrent with a hyperendemic period in meningococcal meningitis. The expansion was marked by specific antigen peptide allele combinations. We observed preliminary evidence for decreasing 4CMenB vaccine coverage in the post-hyperendemic period.

Patterns of inter- and intrasubspecific homologous recombination inform eco-evolutionary dynamics of Xylella fastidiosa

High rates of homologous recombination (HR) in the bacterial plant pathogen Xylella fastidiosa have been previously detected. This study aimed to determine the extent and explore the ecological significance of HR in the genomes of recombinants experimentally generated by natural transformation and wild-type isolates. Both sets of strains displayed widespread HR and similar average size of recombined fragments consisting of random events (2-10 kb) of inter- and intrasubspecific recombination. A significantly higher proportion and greater lengths (>10 kb, maximum 31.5 kb) of recombined fragments were observed in subsp. morus and in strains isolated in Europe from intercepted coffee plants shipped from the Americas. Such highly recombinant strains pose a serious risk of emergence of novel variants, as genetically distinct and formerly geographically isolated genotypes are brought in close proximity by global trade. Recently recombined regions in wild-type strains included genes involved in regulation and signaling, host colonization, nutrient acquisition, and host evasion, all fundamental traits for X. fastidiosa ecology. Identification of four recombinant loci shared between wild-type and experimentally generated recombinants suggests potential hotspots of recombination in this naturally competent pathogen. These findings provide insights into evolutionary forces possibly affecting the adaptive potential to colonize the host environments of X. fastidiosa.

The comparative population genetics of Neisseria meningitidis and Neisseria gonorrhoeae

Neisseria meningitidis and N. gonorrhoeae are closely related pathogenic bacteria. To compare their population genetics, we compiled a dataset of 1,145 genes found across 20 N. meningitidis and 15 N. gonorrhoeae genomes. We find that N. meningitidis is seven-times more diverse than N. gonorrhoeae in their combined core genome. Both species have acquired the majority of their diversity by recombination with divergent strains, however, we find that N. meningitidis has acquired more of its diversity by recombination than N. gonorrhoeae. We find that linkage disequilibrium (LD) declines rapidly across the genomes of both species. Several observations suggest that N. meningitidis has a higher effective population size than N. gonorrhoeae; it is more diverse, the ratio of non-synonymous to synonymous polymorphism is lower, and LD declines more rapidly to a lower asymptote in N. meningitidis. The two species share a modest amount of variation, half of which seems to have been acquired by lateral gene transfer and half from their common ancestor. We investigate whether diversity varies across the genome of each species and find that it does. Much of this variation is due to different levels of lateral gene transfer. However, we also find some evidence that the effective population size varies across the genome. We test for adaptive evolution in the core genome using a McDonald–Kreitman test and by considering the diversity around non-synonymous sites that are fixed for different alleles in the two species. We find some evidence for adaptive evolution using both approaches.

The Patchy Distribution of Restriction⁻Modification System Genes and the Conservation of Orphan Methyltransferases in Halobacteria

Restriction⁻modification (RM) systems in bacteria are implicated in multiple biological roles ranging from defense against parasitic genetic elements, to selfish addiction cassettes, and barriers to gene transfer and lineage homogenization. In bacteria, DNA-methylation without cognate restriction also plays important roles in DNA replication, mismatch repair, protein expression, and in biasing DNA uptake. Little is known about archaeal RM systems and DNA methylation. To elucidate further understanding for the role of RM systems and DNA methylation in Archaea, we undertook a survey of the presence of RM system genes and related genes, including orphan DNA methylases, in the halophilic archaeal class Halobacteria. Our results reveal that some orphan DNA methyltransferase genes were highly conserved among lineages indicating an important functional constraint, whereas RM systems demonstrated patchy patterns of presence and absence. This irregular distribution is due to frequent horizontal gene transfer and gene loss, a finding suggesting that the evolution and life cycle of RM systems may be best described as that of a selfish genetic element. A putative target motif (CTAG) of one of the orphan methylases was underrepresented in all of the analyzed genomes, whereas another motif (GATC) was overrepresented in most of the haloarchaeal genomes, particularly in those that encoded the cognate orphan methylase.

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

The expansion of hypervirulent sequence type 4821 clonal complex (CC4821) lineage Neisseria meningitidis bacteria has led to a shift in meningococcal disease epidemiology in China, from serogroup A (MenA) to MenC. Knowledge of the evolution and genetic origin of the emergent MenC strains is limited. In this study, we subjected 76 CC4821 isolates collected across China during 1972–1977 and 2005–2013 to phylogenetic analysis, traditional genotyping, or both. We show that successive recombination events within genes encoding surface antigens and acquisition of quinolone resistance mutations possibly played a role in the emergence of CC4821 as an epidemic clone in China. MenC and MenB CC4821 strains have spread across China and have been detected in several countries in different continents. Capsular switches involving serogroups B and C occurred among epidemic strains, raising concerns regarding possible increases in MenB disease, given that vaccines in use in China do not protect against MenB.

A statistical method to identify recombination in bacterial genomes based on SNP incompatibility

BACKGROUND

Phylogeny estimation for bacteria is likely to reflect their true evolutionary histories only if they are highly clonal. However, recombination events could occur during evolution for some species. The reconstruction of phylogenetic trees from an alignment without considering recombination could be misleading, since the relationships among strains in some parts of the genome might be different than in others. Using a single, global tree can create the appearance of homoplasy in recombined regions. Hence, the identification of recombination breakpoints is essential to better understand the evolutionary relationships of isolates among a bacterial population.

RESULTS

Previously, we have developed a method (called ACR) to detect potential breakpoints in an alignment by evaluating compatibility of polymorphic sites in a sliding window. To assess the statistical significance of candidate breakpoints, we propose an extension of the algorithm (ptACR) that applies a permutation test to generate a null distribution for comparing the average local compatibility. The performance of ptACR is evaluated on both simulated and empirical datasets. ptACR is shown to have similar sensitivity (true positive rate) but a lower false positive rate and higher F1 score compared to basic ACR. When used to analyze a collection of clinical isolates of Staphylococcus aureus, ptACR finds clear evidence of recombination events in this bacterial pathogen, and is able to identify statistically significant boundaries of chromosomal regions with distinct phylogenies.

CONCLUSIONS

ptACR is an accurate and efficient method for identifying genomic regions affected by recombination in bacterial genomes.

Moraxella catarrhalis Restriction-Modification Systems Are Associated with Phylogenetic Lineage and Disease

Moraxella catarrhalis is a human-adapted pathogen, and a major cause of otitis media (OM) and exacerbations of chronic obstructive pulmonary disease. The species is comprised of two main phylogenetic lineages, RB1 and RB2/3. Restriction–modification (R-M) systems are among the few lineage-associated genes identified in other bacterial genera and have multiple functions including defense against foreign invading DNA, maintenance of speciation, and epigenetic regulation of gene expression. Here, we define the repertoire of R-M systems in 51 publicly available M. catarrhalis genomes and report their distribution among M. catarrhalis phylogenetic lineages. An association with phylogenetic lineage (RB1 or RB2/3) was observed for six R-M systems, which may contribute to the evolution of the lineages by restricting DNA transformation. In addition, we observed a relationship between a mutually exclusive Type I R-M system and a Type III R-M system at a single locus conserved throughout a geographically and clinically diverse set of M. catarrhalis isolates. The Type III R-M system at this locus contains the phase-variable Type III DNA methyltransferase, modM, which controls a phasevarion (phase-variable regulon). We observed an association between modM presence and OM-associated middle ear isolates, indicating a potential role for ModM-mediated epigenetic regulation in OM pathobiology.

Genotypic and Phenotypic Characterization of the O-Linked Protein Glycosylation System Reveals High Glycan Diversity in Paired Meningococcal Carriage Isolates

Species within the genus Neisseria display significant glycan diversity associated with the O-linked protein glycosylation (pgl) systems due to phase variation and polymorphic genes and gene content. The aim of this study was to examine in detail the pgl genotype and glycosylation phenotype in meningococcal isolates and the changes occurring during short-term asymptomatic carriage. Paired meningococcal isolates derived from 50 asymptomatic meningococcal carriers, taken about 2 months apart, were analyzed with whole-genome sequencing. The O-linked protein glycosylation genes were characterized in detail using the Genome Comparator tool at the https://pubmlst.org/ database. Immunoblotting with glycan-specific antibodies (Abs) was used to investigate the protein glycosylation phenotype. All major pgl locus polymorphisms identified in Neisseria meningitidis to date were present in our isolate collection, with the variable presence of pglG and pglH, both in combination with either pglB or pglB2 We identified significant changes and diversity in the pgl genotype and/or glycan phenotype in 96% of the paired isolates. There was also a high degree of glycan microheterogeneity, in which different variants of glycan structures were found at a given glycoprotein. The main mechanism responsible for the observed differences was phase-variable expression of the involved glycosyltransferases and the O-acetyltransferase. To our knowledge, this is the first characterization of the pgl genotype and glycosylation phenotype in a larger strain collection. This report thus provides important insight into glycan diversity in N. meningitidis and into the phase variability changes that influence the expressed glycoform repertoire during meningococcal carriage.IMPORTANCE Bacterial meningitis is a serious global health problem, and one of the major causative organisms is Neisseria meningitidis, which is also a common commensal in the upper respiratory tract of healthy humans. In bacteria, numerous loci involved in biosynthesis of surface-exposed antigenic structures that are involved in the interaction between bacteria and host are frequently subjected to homologous recombination and phase variation. These mechanisms are well described in Neisseria, and phase variation provides the ability to change these structures reversibly in response to the environment. Protein glycosylation systems are becoming widely identified in bacteria, and yet little is known about the mechanisms and evolutionary forces influencing glycan composition during carriage and disease.

Investigation of potential targets of Porphyromonas CRISPRs among the genomes of Porphyromonas species

The oral bacterial species Porphyromonas gingivalis, a periodontal pathogen, has plastic genomes that may be driven by homologous recombination with exogenous deoxyribonucleic acid (DNA) that is incorporated by natural transformation and conjugation. However, bacteriophages and plasmids, both of which are main resources of exogenous DNA, do not exist in the known P. gingivalis genomes. This could be associated with an adaptive immunity system conferred by clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated (cas) genes in P. gingivalis as well as innate immune systems such as a restriction-modification system. In a previous study, few immune targets were predicted for P. gingivalis CRISPR/Cas. In this paper, we analyzed 51 P. gingivalis genomes, which were newly sequenced, and publicly available genomes of 13 P. gingivalis and 46 other Porphyromonas species. We detected 6 CRISPR/Cas types (classified by sequence similarity of repeat) in P. gingivalis and 12 other types in the remaining species. The Porphyromonas CRISPR spacers with potential targets in the genus Porphyromonas were approximately 23 times more abundant than those with potential targets in other genus taxa (1,720/6,896 spacers vs. 74/6,896 spacers). Porphyromonas CRISPR/Cas may be involved in genome plasticity by exhibiting selective interference against intra- and interspecies nucleic acids.

Dynamics and impact of homologous recombination on the evolution of Legionella pneumophila

Legionella pneumophila is an environmental bacterium and the causative agent of Legionnaires' disease. Previous genomic studies have shown that recombination accounts for a high proportion (>96%) of diversity within several major disease-associated sequence types (STs) of L. pneumophila. This suggests that recombination represents a potentially important force shaping adaptation and virulence. Despite this, little is known about the biological effects of recombination in L. pneumophila, particularly with regards to homologous recombination (whereby genes are replaced with alternative allelic variants). Using newly available population genomic data, we have disentangled events arising from homologous and non-homologous recombination in six major disease-associated STs of L. pneumophila (subsp. pneumophila), and subsequently performed a detailed characterisation of the dynamics and impact of homologous recombination. We identified genomic "hotspots" of homologous recombination that include regions containing outer membrane proteins, the lipopolysaccharide (LPS) region and Dot/Icm effectors, which provide interesting clues to the selection pressures faced by L. pneumophila. Inference of the origin of the recombined regions showed that isolates have most frequently imported DNA from isolates belonging to their own clade, but also occasionally from other major clades of the same subspecies. This supports the hypothesis that the possibility for horizontal exchange of new adaptations between major clades of the subspecies may have been a critical factor in the recent emergence of several clinically important STs from diverse genomic backgrounds. However, acquisition of recombined regions from another subspecies, L. pneumophila subsp. fraseri, was rarely observed, suggesting the existence of a recombination barrier and/or the possibility of ongoing speciation between the two subspecies. Finally, we suggest that multi-fragment recombination may occur in L. pneumophila, whereby multiple non-contiguous segments that originate from the same molecule of donor DNA are imported into a recipient genome during a single episode of recombination.

Emergence of a new Neisseria meningitidis clonal complex 11 lineage 11.2 clade as an effective urogenital pathogen

Neisseria meningitidis (Nm) clonal complex 11 (cc11) lineage is a hypervirulent pathogen responsible for outbreaks of invasive meningococcal disease, including among men who have sex with men, and is increasingly associated with urogenital infections. Recently, clusters of Nm urethritis have emerged primarily among heterosexual males in the United States. We determined that nonencapsulated meningococcal isolates from an ongoing Nm urethritis outbreak among epidemiologically unrelated men in Columbus, Ohio, are linked to increased Nm urethritis cases in multiple US cities, including Atlanta and Indianapolis, and that they form a unique clade (the US Nm urethritis clade, US_NmUC). The isolates belonged to the cc11 lineage 11.2/ET-15 with fine type of PorA P1.5-1, 10-8; FetA F3-6; PorB 2-2 and express a unique FHbp allele. A common molecular fingerprint of US_NmUC isolates was an IS1301 element in the intergenic region separating the capsule ctr-css operons and adjacent deletion of cssA/B/C and a part of csc, encoding the serogroup C capsule polymerase. This resulted in the loss of encapsulation and intrinsic lipooligosaccharide sialylation that may promote adherence to mucosal surfaces. Furthermore, we detected an IS1301-mediated inversion of an ∼20-kb sequence near the cps locus. Surprisingly, these isolates had acquired by gene conversion the complete gonococcal denitrification norB-aniA gene cassette, and strains grow well anaerobically. The cc11 US_NmUC isolates causing urethritis clusters in the United States may have adapted to a urogenital environment by loss of capsule and gene conversion of the Neisseria gonorrheae norB-aniA cassette promoting anaerobic growth.

Gene flow, recombination, and positive selection in Stenotrophomonas maltophilia: mechanisms underlying the diversity of the widespread opportunistic pathogen

Stenotrophomonas maltophilia is a global multidrug-resistant human opportunistic pathogen in clinical environments. Stenotrophomonas maltophilia is also ubiquitous in aqueous environments, soil, and plants. Various molecular typing methods have revealed that S. maltophilia exhibits high levels of phenotypic and genotypic diversity. However, information regarding the genomic diversity within S. maltophilia and the corresponding genetic mechanisms resulting in said diversity remain scarce. The genome sequences of 17 S. maltophilia strains were selected to investigate the mechanisms contributing to genetic diversity at the genome level. The core and large pan-genomes of the species were first estimated, resulting in a large, open pan-genome. A species phylogeny was also reconstructed based on 344 orthologous genes with one copy per genome, and the contribution of four evolutionary mechanisms to the species genome diversity was quantified: 15%-35% of the genes showed evidence for recombination, 0%-25% of the genes in one genome were likely gained, 0%-44% of the genes in some genomes were likely lost, and less than 0.3% of the genes in a genome were under positive selection pressures. We observed that, among the four main mechanisms, homologous recombination plays a key role in maintaining diversity in S. maltophilia. In this study, we provide an overview of evolution in S. maltophilia to provide a better understanding of its evolutionary dynamics and its relationship with genome diversity.

Genomic Investigation Reveals Highly Conserved, Mosaic, Recombination Events Associated with Capsular Switching among Invasive Neisseria meningitidis Serogroup W Sequence Type (ST)-11 Strains

Neisseria meningitidis is an important cause of meningococcal disease globally. Sequence type (ST)-11 clonal complex (cc11) is a hypervirulent meningococcal lineage historically associated with serogroup C capsule and is believed to have acquired the W capsule through a C to W capsular switching event. We studied the sequence of capsule gene cluster (cps) and adjoining genomic regions of 524 invasive W cc11 strains isolated globally. We identified recombination breakpoints corresponding to two distinct recombination events within W cc11: A 8.4-kb recombinant region likely acquired from W cc22 including the sialic acid/glycosyl-transferase gene, csw resulted in a C→W change in capsular phenotype and a 13.7-kb recombinant segment likely acquired from Y cc23 lineage includes 4.5 kb of cps genes and 8.2 kb downstream of the cps cluster resulting in allelic changes in capsule translocation genes. A vast majority of W cc11 strains (497/524, 94.8%) retain both recombination events as evidenced by sharing identical or very closely related capsular allelic profiles. These data suggest that the W cc11 capsular switch involved two separate recombination events and that current global W cc11 meningococcal disease is caused by strains bearing this mosaic capsular switch.

Distribution of the type III DNA methyltransferases modA, modB and modD among Neisseria meningitidis genotypes: implications for gene regulation and virulence

Neisseria meningitidis is a human-specific bacterium that varies in invasive potential. All meningococci are carried in the nasopharynx, and most genotypes are very infrequently associated with invasive meningococcal disease; however, those belonging to the ‘hyperinvasive lineages’ are more frequently associated with sepsis or meningitis. Genome content is highly conserved between carriage and disease isolates, and differential gene expression has been proposed as a major determinant of the hyperinvasive phenotype. Three phase variable DNA methyltransferases (ModA, ModB and ModD), which mediate epigenetic regulation of distinct phase variable regulons (phasevarions), have been identified in N. meningitidis. Each mod gene has distinct alleles, defined by their Mod DNA recognition domain, and these target and methylate different DNA sequences, thereby regulating distinct gene sets. Here 211 meningococcal carriage and >1,400 disease isolates were surveyed for the distribution of meningococcal mod alleles. While modA11-12 and modB1-2 were found in most isolates, rarer alleles (e.g., modA15, modB4, modD1-6) were specific to particular genotypes as defined by clonal complex. This suggests that phase variable Mod proteins may be associated with distinct phenotypes and hence invasive potential of N. meningitidis strains.

Global epidemiology of capsular group W meningococcal disease (1970-2015): Multifocal emergence and persistence of hypervirulent sequence type (ST)-11 clonal complex

Following an outbreak in Mecca Saudi Arabia in 2000, meningococcal strains expressing capsular group W (W) emerged as a major cause of invasive meningococcal disease (IMD) worldwide. The Saudi Arabian outbreak strain (Hajj clone) belonging to the ST-11 clonal complex (cc11) is similar to W cc11 causing occasional sporadic disease before 2000. Since 2000, W cc11 has caused large meningococcal disease epidemics in the African meningitis belt and endemic disease in South America, Europe and China. Traditional molecular epidemiologic typing suggested that a majority of current W cc11 burden represented global spread of the Hajj clone. However, recent whole genome sequencing (WGS) analyses revealed significant genetic heterogeneity among global W cc11 strains. While continued spread of the Hajj clone occurs in the Middle East, the meningitis belt and South Africa have co-circulation of the Hajj clone and other unrelated W cc11 strains. Notably, South America, the UK, and France share a genetically distinct W cc11 strain. Other W lineages persist in low numbers in Europe, North America and the meningitis belt. In summary, WGS is helping to unravel the complex genomic epidemiology of group W meningococcal strains. Wider application of WGS and strengthening of global IMD surveillance is necessary to monitor the continued evolution of group W lineages.

Plasmid pEC156, a Naturally Occurring Escherichia coli Genetic Element That Carries Genes of the EcoVIII Restriction-Modification System, Is Mobilizable among Enterobacteria

Type II restriction-modification systems are ubiquitous in prokaryotes. Some of them are present in naturally occurring plasmids, which may facilitate the spread of these systems in bacterial populations by horizontal gene transfer. However, little is known about the routes of their dissemination. As a model to study this, we have chosen an Escherichia coli natural plasmid pEC156 that carries the EcoVIII restriction modification system. The presence of this system as well as the cis-acting cer site involved in resolution of plasmid multimers determines the stable maintenance of pEC156 not only in Escherichia coli but also in other enterobacteria. We have shown that due to the presence of oriT-type F and oriT-type R64 loci it is possible to mobilize pEC156 by conjugative plasmids (F and R64, respectively). The highest mobilization frequency was observed when pEC156-derivatives were transferred between Escherichia coli strains, Enterobacter cloacae and Citrobacter freundii representing coliform bacteria. We found that a pEC156-derivative with a functional EcoVIII restriction-modification system was mobilized in enterobacteria at a frequency lower than a plasmid lacking this system. In addition, we found that bacteria that possess the EcoVIII restriction-modification system can efficiently release plasmid content to the environment. We have shown that E. coli cells can be naturally transformed with pEC156-derivatives, however, with low efficiency. The transformation protocol employed neither involved chemical agents (e.g. CaCl₂) nor temperature shift which could induce plasmid DNA uptake.

DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates

The Gram-negative bacterium Neisseria meningitidis features extensive genetic variability. To present, proposed virulence genotypes are also detected in isolates from asymptomatic carriers, indicating more complex mechanisms underlying variable colonization modes of N. meningitidis. We applied the Single Molecule, Real-Time (SMRT) sequencing method from Pacific Biosciences to assess the genome-wide DNA modification profiles of two genetically related N. meningitidis strains, both of serogroup A. The resulting DNA methylomes revealed clear divergences, represented by the detection of shared and of strain-specific DNA methylation target motifs. The positional distribution of these methylated target sites within the genomic sequences displayed clear biases, which suggest a functional role of DNA methylation related to the regulation of genes. DNA methylation in N. meningitidis has a likely underestimated potential for variability, as evidenced by a careful analysis of the ORF status of a panel of confirmed and predicted DNA methyltransferase genes in an extended collection of N. meningitidis strains of serogroup A. Based on high coverage short sequence reads, we find phase variability as a major contributor to the variability in DNA methylation. Taking into account the phase variable loci, the inferred functional status of DNA methyltransferase genes matched the observed methylation profiles. Towards an elucidation of presently incompletely characterized functional consequences of DNA methylation in N. meningitidis, we reveal a prominent colocalization of methylated bases with Single Nucleotide Polymorphisms (SNPs) detected within our genomic sequence collection. As a novel observation we report increased mutability also at 6mA methylated nucleotides, complementing mutational hotspots previously described at 5mC methylated nucleotides. These findings suggest a more diverse role of DNA methylation and Restriction-Modification (RM) systems in the evolution of prokaryotic genomes.

Genomic Epidemiology of Hypervirulent Serogroup W, ST-11 Neisseria meningitidis

Neisseria meningitidis is a leading bacterial cause of sepsis and meningitis globally with dynamic strain distribution over time. Beginning with an epidemic among Hajj pilgrims in 2000, serogroup W (W) sequence type (ST) 11 emerged as a leading cause of epidemic meningitis in the African ‘meningitis belt’ and endemic cases in South America, Europe, Middle East and China. Previous genotyping studies were unable to reliably discriminate sporadic W ST-11 strains in circulation since 1970 from the Hajj outbreak strain (Hajj clone). It is also unclear what proportion of more recent W ST-11 disease clusters are caused by direct descendants of the Hajj clone. Whole genome sequences of 270 meningococcal strains isolated from patients with invasive meningococcal disease globally from 1970 to 2013 were compared using whole genome phylogenetic and major antigen-encoding gene sequence analyses. We found that all W ST-11 strains were descendants of an ancestral strain that had undergone unique capsular switching events. The Hajj clone and its descendants were distinct from other W ST-11 strains in that they shared a common antigen gene profile and had undergone recombination involving virulence genes encoding factor H binding protein, nitric oxide reductase, and nitrite reductase. These data demonstrate that recent acquisition of a distinct antigen-encoding gene profile and variations in meningococcal virulence genes was associated with the emergence of the Hajj clone. Importantly, W ST-11 strains unrelated to the Hajj outbreak contribute a significant proportion of W ST-11 cases globally. This study helps illuminate genomic factors associated with meningococcal strain emergence and evolution.

•

Genomic characterization of serogroup W ST-11 of Neisseria meningitidis.

.

•

Epidemic W ST-11 strain (Hajj clone) emerged through recombination affecting virulence genes.

•

Both the Hajj clone and W ST-11 strains unrelated to the Hajj outbreak have persisted globally.

Neisseria meningitidis, a bacterial cause of frequently fatal brain (meningitis) and blood stream (sepsis) infections, has variable strain distribution over time. Serogroup W sequence type 11 (W ST-11) lineage is associated on one hand with strains causing only rare (sporadic) disease cases, and the Hajj clone – a major global cause of epidemic and endemic meningococcal disease. In this study we analyzed complete genome sequences of a global collection of 270 W ST-11 isolates causing meningococcal disease from 1970-2013. The Hajj clone acquired novel gene sequences within genes involved in nitrogen metabolism (nitrogen oxide reductase, nitrite reductase) and evasion of human immune response (factor H binding protein). These genes may be the cause of increased virulence of the Hajj clone and can be used to trace continuing spread of the clone. These results shed light on mechanisms of meningococcal strain emergence.

The Contribution of Genetic Recombination to CRISPR Array Evolution

CRISPR (clustered regularly interspaced short palindromic repeats) is a microbial immune system against foreign DNA. Recognition sequences (spacers) encoded within the CRISPR array mediate the immune reaction in a sequence-specific manner. The known mechanisms for the evolution of CRISPR arrays include spacer acquisition from foreign DNA elements at the time of invasion and array erosion through spacer deletion. Here, we consider the contribution of genetic recombination between homologous CRISPR arrays to the evolution of spacer repertoire. Acquisition of spacers from exogenic arrays via recombination may confer the recipient with immunity against unencountered antagonists. For this purpose, we develop a novel method for the detection of recombination in CRISPR arrays by modeling the spacer order in arrays from multiple strains from the same species. Because the evolutionary signal of spacer recombination may be similar to that of pervasive spacer deletions or independent spacer acquisition, our method entails a robustness analysis of the recombination inference by a statistical comparison to resampled and perturbed data sets. We analyze CRISPR data sets from four bacterial species: two Gammaproteobacteria species harboring CRISPR type I and two Streptococcus species harboring CRISPR type II loci. We find that CRISPR array evolution in Escherichia coli and Streptococcus agalactiae can be explained solely by vertical inheritance and differential spacer deletion. In Pseudomonas aeruginosa, we find an excess of single spacers potentially incorporated into the CRISPR locus during independent acquisition events. In Streptococcus thermophilus, evidence for spacer acquisition by recombination is present in 5 out of 70 strains. Genetic recombination has been proposed to accelerate adaptation by combining beneficial mutations that arose in independent lineages. However, for most species under study, we find that CRISPR evolution is shaped mainly by spacer acquisition and loss rather than recombination. Since the evolution of spacer content is characterized by a rapid turnover, it is likely that recombination is not beneficial for improving phage resistance in the strains under study, or that it cannot be detected in the resolution of intraspecies comparisons.

Not seeing the forest for the trees: size of the minimum spanning trees (MSTs) forest and branch significance in MST-based phylogenetic analysis

Trees, including minimum spanning trees (MSTs), are commonly used in phylogenetic studies. But, for the research community, it may be unclear that the presented tree is just a hypothesis, chosen from among many possible alternatives. In this scenario, it is important to quantify our confidence in both the trees and the branches/edges included in such trees. In this paper, we address this problem for MSTs by introducing a new edge betweenness metric for undirected and weighted graphs. This spanning edge betweenness metric is defined as the fraction of equivalent MSTs where a given edge is present. The metric provides a per edge statistic that is similar to that of the bootstrap approach frequently used in phylogenetics to support the grouping of taxa. We provide methods for the exact computation of this metric based on the well known Kirchhoff’s matrix tree theorem. Moreover, we implement and make available a module for the PHYLOViZ software and evaluate the proposed metric concerning both effectiveness and computational performance. Analysis of trees generated using multilocus sequence typing data (MLST) and the goeBURST algorithm revealed that the space of possible MSTs in real data sets is extremely large. Selection of the edge to be represented using bootstrap could lead to unreliable results since alternative edges are present in the same fraction of equivalent MSTs. The choice of the MST to be presented, results from criteria implemented in the algorithm that must be based in biologically plausible models.

Chaperone-usher fimbriae in a diverse selection of Gallibacterium genomes

Background

Fimbriae are bacterial cell surface organelles involved in the pathogenesis of many bacterial species, including Gallibacterium anatis, in which a F17-like fimbriae of the chaperone-usher (CU) family was recently shown to be an important virulence factor and vaccine candidate. To reveal the distribution and variability of CU fimbriae 22 genomes of the avian host-restricted bacteria Gallibacterium spp. were investigated. Fimbrial clusters were classified using phylogeny-based and conserved domain (CD) distribution-based approaches. To characterize the fimbriae in depth evolutionary analysis and in vitro expression of the most prevalent fimbrial clusters was performed.

Results

Overall 48 CU fimbriae were identified in the genomes of the examined Gallibacterium isolates. All fimbriae were assigned to γ4 clade of the CU fimbriae of Gram-negative bacteria and were organized in four-gene clusters encoding a putative major fimbrial subunit, a chaperone, an usher and a fimbrial adhesin. Five fimbrial clusters (Flf-Flf4) and eight conserved domain groups were defined to accommodate the identified fimbriae. Although, the number of different fimbrial clusters in individual Gallibacterium genomes was low, there was substantial amino acid sequence variability in the major fimbrial subunit and the adhesin proteins. The distribution of CDs among fimbrial clusters, analysis of their flanking regions, and evolutionary comparison of the strains revealed that Gallibacterium fimbrial clusters likely underwent evolutionary divergence resulting in highly host adapted and antigenically variable fimbriae. In vitro, only the fimbrial subunit FlfA was expressed in most Gallibacterium strains encoding this protein. The absence or scarce expression of the two other common fimbrial subunits (Flf1A and Flf3A) indicates that their expression may require other in vitro or in vivo conditions.

Conclusions

This is the first approach establishing a systematic fimbria classification system within Gallibacterium spp., which indicates a species-wide distribution of γ₄ CU fimbriae among a diverse collection of Gallibacterium isolates. The expression of only one out of up to three fimbriae present in the individual genomes in vitro suggests that fimbriae expression in Gallibacterium is highly regulated. This information is important for future attempts to understand the role of Gallibacterium fimbriae in pathogenesis, and may prove useful for improved control of Gallibacterium infections in chickens.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1093) contains supplementary material, which is available to authorized users.

Emergence of a new epidemic Neisseria meningitidis serogroup A Clone in the African meningitis belt: high-resolution picture of genomic changes that mediate immune evasion

In the African “meningitis belt,” outbreaks of meningococcal meningitis occur in cycles, representing a model for the role of host-pathogen interactions in epidemic processes. The periodicity of the epidemics is not well understood, nor is it currently possible to predict them. In our longitudinal colonization and disease surveys, we have observed waves of clonal replacement with the same serogroup, suggesting that immunity to noncapsular antigens plays a significant role in natural herd immunity. Here, through comparative genomic analysis of 100 meningococcal isolates, we provide a high-resolution view of the evolutionary changes that occurred during clonal replacement of a hypervirulent meningococcal clone (ST-7) by a descendant clone (ST-2859). We show that the majority of genetic changes are due to homologous recombination of laterally acquired DNA, with more than 20% of these events involving acquisition of DNA from other species. Signals of adaptation to evade herd immunity were indicated by genomic hot spots of recombination. Most striking is the high frequency of changes involving the pgl locus, which determines the glycosylation patterns of major protein antigens. High-frequency changes were also observed for genes involved in the regulation of pilus expression and the synthesis of Maf3 adhesins, highlighting the importance of these surface features in host-pathogen interaction and immune evasion.

While established meningococcal capsule polysaccharide vaccines are protective through the induction of anticapsular antibodies, findings of our longitudinal studies in the African meningitis belt have indicated that immunity to noncapsular antigens plays a significant role in natural herd immunity. Our results show that meningococci evade herd immunity through the rapid homologous replacement of just a few key genomic loci that affect noncapsular cell surface components. Identification of recombination hot spots thus represents an eminent approach to gain insight into targets of protective natural immune responses. Moreover, our results highlight the role of the dynamics of the protein glycosylation repertoire in immune evasion by Neisseria meningitidis. These results have major implications for the design of next-generation protein-based subunit vaccines.

Defining the estimated core genome of bacterial populations using a Bayesian decision model

The bacterial core genome is of intense interest and the volume of whole genome sequence data in the public domain available to investigate it has increased dramatically. The aim of our study was to develop a model to estimate the bacterial core genome from next-generation whole genome sequencing data and use this model to identify novel genes associated with important biological functions. Five bacterial datasets were analysed, comprising 2096 genomes in total. We developed a Bayesian decision model to estimate the number of core genes, calculated pairwise evolutionary distances (p-distances) based on nucleotide sequence diversity, and plotted the median p-distance for each core gene relative to its genome location. We designed visually-informative genome diagrams to depict areas of interest in genomes. Case studies demonstrated how the model could identify areas for further study, e.g. 25% of the core genes with higher sequence diversity in the Campylobacter jejuni and Neisseria meningitidis genomes encoded hypothetical proteins. The core gene with the highest p-distance value in C. jejuni was annotated in the reference genome as a putative hydrolase, but further work revealed that it shared sequence homology with beta-lactamase/metallo-beta-lactamases (enzymes that provide resistance to a range of broad-spectrum antibiotics) and thioredoxin reductase genes (which reduce oxidative stress and are essential for DNA replication) in other C. jejuni genomes. Our Bayesian model of estimating the core genome is principled, easy to use and can be applied to large genome datasets. This study also highlighted the lack of knowledge currently available for many core genes in bacterial genomes of significant global public health importance.

Whole genome sequencing has revolutionised the study of pathogenic microorganisms. It has also become so affordable that hundreds of samples can reasonably be sequenced in an individual project, creating a wealth of data. Estimating the bacterial core genome – traditionally defined as those genes present in all genomes – is an important initial step in population genomics analyses. We developed a simple statistical model to estimate the number of core genes in a bacterial genome dataset, calculated pairwise evolutionary distances (p-distances) based on differences among nucleotide sequences, and plotted the median p-distance for each core gene relative to its genome location. Low p-distance values indicate highly-conserved genes; high values suggest genes under selection and/or undergoing recombination. The genome diagrams depict areas of interest in genomes that can be explored in further detail. Using our method, we analysed five bacterial species comprising a total of 2096 genomes. This revealed new information related to antibiotic resistance and virulence for two bacterial species and demonstrated that the function of many core genes in bacteria is still unknown. Our model provides a highly-accessible, publicly-available tool to use on the vast quantities of genome sequence data now available.

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.

Horizontal transfer and gene conversion as an important driving force in shaping the landscape of mitochondrial introns

Group I introns are highly dynamic and mobile, featuring extensive presence-absence variation and widespread horizontal transfer. Group I introns can invade intron-lacking alleles via intron homing powered by their own encoded homing endonuclease gene (HEG) after horizontal transfer or via reverse splicing through an RNA intermediate. After successful invasion, the intron and HEG are subject to degeneration and sequential loss. It remains unclear whether these mechanisms can fully address the high dynamics and mobility of group I introns. Here, we found that HEGs undergo a fast gain-and-loss turnover comparable with introns in the yeast mitochondrial 21S-rRNA gene, which is unexpected, as the intron and HEG are generally believed to move together as a unit. We further observed extensively mosaic sequences in both the introns and HEGs, and evidence of gene conversion between HEG-containing and HEG-lacking introns. Our findings suggest horizontal transfer and gene conversion can accelerate HEG/intron degeneration and loss, or rescue and propagate HEG/introns, and ultimately result in high HEG/intron turnover rate. Given that up to 25% of the yeast mitochondrial genome is composed of introns and most mitochondrial introns are group I introns, horizontal transfer and gene conversion could have served as an important mechanism in introducing mitochondrial intron diversity, promoting intron mobility and consequently shaping mitochondrial genome architecture.

Natural competence and the evolution of DNA uptake specificity

Many bacteria are naturally competent, able to actively transport environmental DNA fragments across their cell envelope and into their cytoplasm. Because incoming DNA fragments can recombine with and replace homologous segments of the chromosome, competence provides cells with a potent mechanism of horizontal gene transfer as well as access to the nutrients in extracellular DNA. This review starts with an introductory overview of competence and continues with a detailed consideration of the DNA uptake specificity of competent proteobacteria in the Pasteurellaceae and Neisseriaceae. Species in these distantly related families exhibit strong preferences for genomic DNA from close relatives, a self-specificity arising from the combined effects of biases in the uptake machinery and genomic overrepresentation of the sequences this machinery prefers. Other competent species tested lack obvious uptake bias or uptake sequences, suggesting that strong convergent evolutionary forces have acted on these two families. Recent results show that uptake sequences have multiple "dialects," with clades within each family preferring distinct sequence variants and having corresponding variants enriched in their genomes. Although the genomic consensus uptake sequences are 12 and 29 to 34 bp, uptake assays have found that only central cores of 3 to 4 bp, conserved across dialects, are crucial for uptake. The other bases, which differ between dialects, make weaker individual contributions but have important cooperative interactions. Together, these results make predictions about the mechanism of DNA uptake across the outer membrane, supporting a model for the evolutionary accumulation and stability of uptake sequences and suggesting that uptake biases may be more widespread than currently thought.