Genealogical typing of Neisseria meningitidis

Multilocus sequence analysis reveals genetic diversity in Staphylococcus aureus isolate of goat with mastitis persistent after treatment with enrofloxacin

Staphylococcus aureus is one of the main bacterial agents responsible for cases of mastitis in ruminants, playing an important role in the persistence and chronicity of diseases treated with antimicrobials. Using the multilocus sequence typing technique, network approaches and study of the population diversity of microorganisms, we performed analyzes of S. aureus (ES-GPM) isolated from goats with persistent mastitis (GPM). The most strains of ES-GPM were categorically different phylogenetically from the others and could be divided into two lineages: one with a majority belonging to ES-GPM and the other to varied strains. These two lineages were separated by 27 nuclear polymorphisms. The 43 strains comprised 22 clonal complexes (CCs), of which the ES-GPM strains were present in CC133, CC5 and a new complex formed by the sequence type 4966. The genetic diversity of some alleles showed be greater diversity and polymorphism than others, such as of the aroE and yqiL genes less than glpF gene. In addition, the sequences ES-GPM to the arc gene and glpF alleles showed the greatest number of mutations for ES-GPM in relation to non-ES-GPM. Therefore, this study identified genetic polymorphisms characteristic of S. aureus isolated from milk of goats diagnosed with persistent mastitis after the failed treatment with the antibiotic enrofloxacin. This study may help in the future to identify and discriminate this agent in cases of mastitis, and with that, the most appropriate antibiotic treatment can be performed in advance of the appearance of persistent mastitis caused by the agent, reducing the chances of premature culling and animal suffering.

Establishment of the European meningococcal strain collection genome library (EMSC-GL) for the 2011 to 2012 epidemiological year

Invasive meningococcal disease surveillance in Europe combines isolate characterisation and epidemiological data to support public health intervention. A representative European Meningococcal Strain Collection (EMSC) of IMD isolates was obtained, and whole genome sequenced to characterise 799 EMSC isolates from the epidemiological year July 2011-June 2012. To establish a genome library (GL), the isolate information was deposited in the pubMLST.org/neisseria database. Genomes were curated and annotated at 2,429 meningococcal loci, including those defining clonal complex, capsule, antigens, and antimicrobial resistance. Most genomes contained genes encoding B (n = 525; 65.7%) or C (n = 163; 20.4%) capsules; isolates were genetically highly diverse, with >20 genomic lineages, five of which comprising 60.7% (n = 485) of isolates. There were >350 antigenic fine-types: 307 were present once, the most frequent (P1.7-2,4:F5-1) comprised 8% (n = 64) of isolates. Each genome was characterised for Bexsero Antigen Sequence Typing (BAST): 25.5% (n = 204) of isolates contained alleles encoding the fHbp and/or the PorA VR1 vaccine component, but most genomes (n = 513; 64.2%) did not contain the NadA component. EMSC-GL will support an integrated surveillance of disease-associated genotypes in Europe, enabling the monitoring of hyperinvasive lineages, outbreak identification, and supporting vaccine programme implementation.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

How clonal are Neisseria species? The epidemic clonality model revisited

The three species Neisseria meningitidis, Neisseria gonorrheae, and Neisseria lactamica are often regarded as highly recombining bacteria. N. meningitidis has been considered a paradigmatic case of the "semiclonal model" or of "epidemic clonality," demonstrating occasional bouts of clonal propagation in an otherwise recombining species. In this model, occasional clonality generates linkage disequilibrium in the short term. In the long run, however, the effects of clonality are countered by recombination. We show that many data are at odds with this proposal and that N. meningitidis fits the criteria that we have proposed for predominant clonal evolution (PCE). We point out that (i) the proposed way to distinguish epidemic clonality from PCE may be faulty and (ii) the evidence of deep phylogenies by microarrays and whole-genome sequencing is at odds with the predictions of the semiclonal model. Last, we revisit the species status of N. meningitidis, N. gonorrheae, and N. lactamica in the light of the PCE model.

ClonalFrameML: efficient inference of recombination in whole bacterial genomes

Recombination is an important evolutionary force in bacteria, but it remains challenging to reconstruct the imports that occurred in the ancestry of a genomic sample. Here we present ClonalFrameML, which uses maximum likelihood inference to simultaneously detect recombination in bacterial genomes and account for it in phylogenetic reconstruction. ClonalFrameML can analyse hundreds of genomes in a matter of hours, and we demonstrate its usefulness on simulated and real datasets. We find evidence for recombination hotspots associated with mobile elements in Clostridium difficile ST6 and a previously undescribed 310kb chromosomal replacement in Staphylococcus aureus ST582. ClonalFrameML is freely available at http://clonalframeml.googlecode.com/.

A gene-by-gene population genomics platform: de novo assembly, annotation and genealogical analysis of 108 representative Neisseria meningitidis genomes

BACKGROUND

Highly parallel, 'second generation' sequencing technologies have rapidly expanded the number of bacterial whole genome sequences available for study, permitting the emergence of the discipline of population genomics. Most of these data are publically available as unassembled short-read sequence files that require extensive processing before they can be used for analysis. The provision of data in a uniform format, which can be easily assessed for quality, linked to provenance and phenotype and used for analysis, is therefore necessary.

RESULTS

The performance of de novo short-read assembly followed by automatic annotation using the pubMLST.org Neisseria database was assessed and evaluated for 108 diverse, representative, and well-characterised Neisseria meningitidis isolates. High-quality sequences were obtained for >99% of known meningococcal genes among the de novo assembled genomes and four resequenced genomes and less than 1% of reassembled genes had sequence discrepancies or misassembled sequences. A core genome of 1600 loci, present in at least 95% of the population, was determined using the Genome Comparator tool. Genealogical relationships compatible with, but at a higher resolution than, those identified by multilocus sequence typing were obtained with core genome comparisons and ribosomal protein gene analysis which revealed a genomic structure for a number of previously described phenotypes. This unified system for cataloguing Neisseria genetic variation in the genome was implemented and used for multiple analyses and the data are publically available in the PubMLST Neisseria database.

CONCLUSIONS

The de novo assembly, combined with automated gene-by-gene annotation, generates high quality draft genomes in which the majority of protein-encoding genes are present with high accuracy. The approach catalogues diversity efficiently, permits analyses of a single genome or multiple genome comparisons, and is a practical approach to interpreting WGS data for large bacterial population samples. The method generates novel insights into the biology of the meningococcus and improves our understanding of the whole population structure, not just disease causing lineages.

Characterization of pseudogenes in members of the order Frankineae

Pseudogenes are defined as non-functional relatives of genes whose protein-coding abilities are lost and are no longer expressed within cells. They are an outcome of accumulation of mutations within a gene whose end product is not essential for survival. Proper investigation of the procedure of pseudogenization is relevant for estimating occurrence of duplications in genomes. Frankineae houses an interesting group of microorganisms, carving a niche in the microbial world. This study was undertaken with the objective of determining the abundance of pseudogenes, understanding strength of purifying selection, investigating evidence of pseudogene expression, and analysing their molecular nature, their origin, evolution and deterioration patterns amongst domain families. Investigation revealed the occurrence of 956 core pFAM families sharing common characteristics indicating co-evolution. WD40, Rve_3, DDE_Tnp_IS240 and phage integrase core domains are larger families, having more pseudogenes, signifying a probability of harmful foreign genes being disabled within transposable elements. High selective pressure depicted that gene families rapidly duplicating and evolving undoubtedly facilitated creation of a number of pseudogenes in Frankineae. Codon usage analysis between protein-coding genes and pseudogenes indicated a wide degree of variation with respect to different factors. Moreover, the majority of pseudogenes were under the effect of purifying selection. Frankineae pseudogenes were under stronger selective constraints, indicating that they were functional for a very long time and became pseudogenes abruptly. The origin and deterioration of pseudogenes has been attributed to selection and mutational pressure acting upon sequences for adapting to stressed soil environments.

Ribosomal proteins as biomarkers for bacterial identification by mass spectrometry in the clinical microbiology laboratory

Whole-cell matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) is a rapid method for identification of microorganisms that is increasingly used in microbiology laboratories. This identification is based on the comparison of the tested isolate mass spectrum with reference databases. Using Neisseria meningitidis as a model organism, we showed that in one of the available databases, the Andromas database, 10 of the 13 species-specific biomarkers correspond to ribosomal proteins. Remarkably, one biomarker, ribosomal protein L32, was subject to inter-strain variability. The analysis of the ribosomal protein patterns of 100 isolates for which whole genome sequences were available, confirmed the presence of inter-strain variability in the molecular weight of 29 ribosomal proteins, thus establishing a correlation between the sequence type (ST) and/or clonal complex (CC) of each strain and its ribosomal protein pattern. Since the molecular weight of three of the variable ribosomal proteins (L30, L31 and L32) was included in the spectral window observed by MALDI-TOF MS in clinical microbiology, i.e., 3640-12000 m/z, we were able by analyzing the molecular weight of these three ribosomal proteins to classify each strain in one of six subgroups, each of these subgroups corresponding to specific STs and/or CCs. Their detection by MALDI-TOF allows therefore a quick typing of N. meningitidis isolates.

Description of an unusual Neisseria meningitidis isolate containing and expressing Neisseria gonorrhoeae-Specific 16S rRNA gene sequences

An apparently rare Neisseria meningitidis isolate containing one copy of a Neisseria gonorrhoeae 16S rRNA gene is described herein. This isolate was identified as N. meningitidis by biochemical identification methods but generated a positive signal with Gen-Probe Aptima assays for the detection of Neisseria gonorrhoeae. Direct 16S rRNA gene sequencing of the purified isolate revealed mixed bases in signature regions that allow for discrimination between N. meningitidis and N. gonorrhoeae. The mixed bases were resolved by sequencing individually PCR-amplified single copies of the genomic 16S rRNA gene. A total of 121 discrete sequences were obtained; 92 (76%) were N. meningitidis sequences, and 29 (24%) were N. gonorrhoeae sequences. Based on the ratio of species-specific sequences, the N. meningitidis strain seems to have replaced one of its four intrinsic 16S rRNA genes with the gonococcal gene. Fluorescence in situ hybridization (FISH) probes specific for meningococcal and gonococcal rRNA were used to demonstrate the expression of the rRNA genes. Interestingly, the clinical isolate described here expresses both N. meningitidis and N. gonorrhoeae 16S rRNA genes, as shown by positive FISH signals with both probes. This explains why the probes for N. gonorrhoeae in the Gen-Probe Aptima assays cross-react with this N. meningitidis isolate. The N. meningitidis isolate described must have obtained N. gonorrhoeae-specific DNA through interspecies recombination.

Microevolutionary analysis of Clostridium difficile genomes to investigate transmission

BACKGROUND

The control of Clostridium difficile infection is a major international healthcare priority, hindered by a limited understanding of transmission epidemiology for these bacteria. However, transmission studies of bacterial pathogens are rapidly being transformed by the advent of next generation sequencing.

RESULTS

Here we sequence whole C. difficile genomes from 486 cases arising over four years in Oxfordshire. We show that we can estimate the times back to common ancestors of bacterial lineages with sufficient resolution to distinguish whether direct transmission is plausible or not. Time depths were inferred using a within-host evolutionary rate that we estimated at 1.4 mutations per genome per year based on serially isolated genomes. The subset of plausible transmissions was found to be highly associated with pairs of patients sharing time and space in hospital. Conversely, the large majority of pairs of genomes matched by conventional typing and isolated from patients within a month of each other were too distantly related to be direct transmissions.

CONCLUSIONS

Our results confirm that nosocomial transmission between symptomatic C. difficile cases contributes far less to current rates of infection than has been widely assumed, which clarifies the importance of future research into other transmission routes, such as from asymptomatic carriers. With the costs of DNA sequencing rapidly falling and its use becoming more and more widespread, genomics will revolutionize our understanding of the transmission of bacterial pathogens.

Microevolution of symbiotic Bradyrhizobium populations associated with soybeans in east North America

Microevolution and origins of Bradyrhizobium populations associated with soybeans at two field sites (A and B, 280 km apart in Canada) with contrasting histories of inoculation was investigated using probabilistic analyses of six core (housekeeping) gene sequences. These analyses supported division of 220 isolates in five lineages corresponding either to B. japonicum groups 1 and 1a or to one of three novel lineages within the genus Bradyrhizobium. None of the isolates from site A and about 20% from site B (the only site with a recent inoculation history) were attributed to inoculation sources. The data suggest that most isolates were of indigenous origin based on sequence analysis of 148 isolates of soybean-nodulating bacteria from native legumes (Amphicarpaea bracteata and Desmodium canadense). Isolates from D. canadense clustered with B. japonicum group 1, whereas those from A. bracteata were placed in two novel lineages encountered at soybean field sites. One of these novel lineages predominated at soybean sites and exhibited a significant clonal expansion likely reflecting selection by the plant host. Homologous recombination events detected in the 35 sequence types from soybean sites had an effect on genetic diversification that was approximately equal to mutation. Interlineage transfer of core genes was infrequent and mostly attributable to gyrB that had a history of frequent recombination. Symbiotic gene sequences (nodC and nifH) of isolates from soybean sites and native legumes clustered in two lineages corresponding to B. japonicum and B. elkani with the inheritance of these genes appearing predominantly by vertical transmission. The data suggest that soybean-nodulating bacteria associated with native legumes represent a novel source of ecologically adapted bacteria for soybean inoculation.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Resolution of a meningococcal disease outbreak from whole-genome sequence data with rapid Web-based analysis methods

The increase in the capacity and reduction in cost of whole-genome sequencing methods present the imminent prospect of such data being used routinely in real time for investigations of bacterial disease outbreaks. For this to be realized, however, it is necessary that generic, portable, and robust analysis frameworks be available, which can be readily interpreted and used in real time by microbiologists, clinicians, and public health epidemiologists. We have achieved this with a set of analysis tools integrated into the PubMLST.org website, which can in principle be used for the analysis of any pathogen. The approach is demonstrated with genomic data from isolates obtained during a well-characterized meningococcal disease outbreak at the University of Southampton, United Kingdom, that occurred in 1997. Whole-genome sequence data were collected, de novo assembled, and deposited into the PubMLST Neisseria BIGSdb database, which automatically annotated the sequences. This enabled the immediate and backwards-compatible classification of the isolates with a number of schemes, including the following: conventional, extended, and ribosomal multilocus sequence typing (MLST, eMLST, and rMLST); antigen gene sequence typing (AGST); analysis based on genes conferring antibiotic susceptibility. The isolates were also compared to a reference isolate belonging to the same clonal complex (ST-11) at 1,975 loci. Visualization of the data with the NeighborNet algorithm, implemented in SplitsTree 4 within the PubMLST website, permitted complete resolution of the outbreak and related isolates, demonstrating that multiple closely related but distinct strains were simultaneously present in asymptomatic carriage and disease, with two causing disease and one responsible for the outbreak itself.

Evolutionary and genomic insights into meningococcal biology

Epidemic disease caused by Neisseria meningitidis, the meningococcus, has been recognized for two centuries, but remains incompletely controlled and understood. There have been dramatic reductions in serogroup A and C meningococcal disease following the introduction of protein-polysaccharide conjugate vaccines, but there is currently no comprehensive vaccine against serogroup B meningococci. Genetic analyses of meningococcal populations have provided many insights into the biology, evolution and pathogenesis of this important pathogen. The meningococcus, and its close relative the gonococcus, are the only pathogenic members of the genus Neisseria, and the invasive propensity of meningococci varies widely, with approximately a dozen 'hyperinvasive lineages' responsible for most disease. Despite this, attempts to identify a 'pathogenome', a subset of genes associated with the invasive phenotypes, have failed; however, genome-wide studies of representative meningococcal isolates using high-throughput sequencing are beginning to provide details on the relationship of invasive phenotype and genotype in this fascinating organism and how this relationship has evolved.

Impact of homologous and non-homologous recombination in the genomic evolution of Escherichia coli

Background

Escherichia coli is an important species of bacteria that can live as a harmless inhabitant of the guts of many animals, as a pathogen causing life-threatening conditions or freely in the non-host environment. This diversity of lifestyles has made it a particular focus of interest for studies of genetic variation, mainly with the aim to understand how a commensal can become a deadly pathogen. Many whole genomes of E. coli have been fully sequenced in the past few years, which offer helpful data to help understand how this important species evolved.

Results

We compared 27 whole genomes encompassing four phylogroups of Escherichia coli (A, B1, B2 and E). From the core-genome we established the clonal relationships between the isolates as well as the role played by homologous recombination during their evolution from a common ancestor. We found strong evidence for sexual isolation between three lineages (A+B1, B2, E), which could be explained by the ecological structuring of E. coli and may represent on-going speciation. We identified three hotspots of homologous recombination, one of which had not been previously described and contains the aroC gene, involved in the essential shikimate metabolic pathway. We also described the role played by non-homologous recombination in the pan-genome, and showed that this process was highly heterogeneous. Our analyses revealed in particular that the genomes of three enterohaemorrhagic (EHEC) strains within phylogroup B1 have converged from originally separate backgrounds as a result of both homologous and non-homologous recombination.

Conclusions

Recombination is an important force shaping the genomic evolution and diversification of E. coli, both by replacing fragments of genes with an homologous sequence and also by introducing new genes. In this study, several non-random patterns of these events were identified which correlated with important changes in the lifestyle of the bacteria, and therefore provide additional evidence to explain the relationship between genomic variation and ecological adaptation.

Whole genome sequencing to investigate the emergence of clonal complex 23 Neisseria meningitidis serogroup Y disease in the United States

In the United States, serogroup Y, ST-23 clonal complex Neisseria meningitidis was responsible for an increase in meningococcal disease incidence during the 1990s. This increase was accompanied by antigenic shift of three outer membrane proteins, with a decrease in the population that predominated in the early 1990s as a different population emerged later in that decade. To understand factors that may have been responsible for the emergence of serogroup Y disease, we used whole genome pyrosequencing to investigate genetic differences between isolates from early and late N. meningitidis populations, obtained from meningococcal disease cases in Maryland in the 1990s. The genomes of isolates from the early and late populations were highly similar, with 1231 of 1776 shared genes exhibiting 100% amino acid identity and an average π_N  =  0.0033 and average π_S  =  0.0216. However, differences were found in predicted proteins that affect pilin structure and antigen profile and in predicted proteins involved in iron acquisition and uptake. The observed changes are consistent with acquisition of new alleles through horizontal gene transfer. Changes in antigen profile due to the genetic differences found in this study likely allowed the late population to emerge due to escape from population immunity. These findings may predict which antigenic factors are important in the cyclic epidemiology of meningococcal disease.

Ion torrent personal genome machine sequencing for genomic typing of Neisseria meningitidis for rapid determination of multiple layers of typing information

Neisseria meningitidis causes invasive meningococcal disease in infants, toddlers, and adolescents worldwide. DNA sequence-based typing, including multilocus sequence typing, analysis of genetic determinants of antibiotic resistance, and sequence typing of vaccine antigens, has become the standard for molecular epidemiology of the organism. However, PCR of multiple targets and consecutive Sanger sequencing provide logistic constraints to reference laboratories. Taking advantage of the recent development of benchtop next-generation sequencers (NGSs) and of BIGSdb, a database accommodating and analyzing genome sequence data, we therefore explored the feasibility and accuracy of Ion Torrent Personal Genome Machine (PGM) sequencing for genomic typing of meningococci. Three strains from a previous meningococcus serogroup B community outbreak were selected to compare conventional typing results with data generated by semiconductor chip-based sequencing. In addition, sequencing of the meningococcal type strain MC58 provided information about the general performance of the technology. The PGM technology generated sequence information for all target genes addressed. The results were 100% concordant with conventional typing results, with no further editing being necessary. In addition, the amount of typing information, i.e., nucleotides and target genes analyzed, could be substantially increased by the combined use of genome sequencing and BIGSdb compared to conventional methods. In the near future, affordable and fast benchtop NGS machines like the PGM might enable reference laboratories to switch to genomic typing on a routine basis. This will reduce workloads and rapidly provide information for laboratory surveillance, outbreak investigation, assessment of vaccine preventability, and antibiotic resistance gene monitoring.

Characterization of Neisseria meningitidis isolates from Egypt using multilocus sequence typing

To characterize Neisseria meningitidis isolates collected from cerebrospinal fluid of meningitis cases in Egypt (1998-2003) as part of surveillance studies, 67 isolates were serogrouped, tested for antibiotic sensitivity and analyzed using multilocus sequence typing (MLST). Results show that isolates expressing serogroup B (50.7%) and serogroup A (34.3%) antigens were predominant in Egypt during the surveillance period, possibly due to suppression of other serogroups by meningococcal vaccines in current use. Intermediate resistance to penicillin was observed in 71% of the isolates, suggesting a need for physicians to shift to third-generation cephalosporins during the empirical treatment of infection. Recurrent lineages of N. meningitidis in Egypt appear to originate from Europe and other Middle Eastern countries. Of 19 sequence types detected, five were unique to Africa and 10 were not observed previously in the MLST database. The information obtained illustrates the changing dynamics of meningitis after vaccine introduction in Egypt.

The effect of iron availability on transcription of the Neisseria meningitidis fHbp gene varies among clonal complexes

Factor H binding protein (fHbp) is a major antigenic component of novel vaccines designed to protect against meningococcal disease. Prediction of the potential coverage of these vaccines is difficult, as fHbp is antigenically variable and levels of expression differ among isolates. Transcriptional regulation of the fHbp gene is poorly understood, although evidence suggests that oxygen availability is involved. In this study iron accessibility was found to affect fHbp transcription. However, regulation differed among meningococcal clonal complexes (ccs). For the majority of isolates, increased iron concentrations upregulated transcription. This effect was enhanced by the presence of a 181 bp insertion element upstream of fHbp, associated with isolates belonging to cc4 and cc5. Conversely, meningococci belonging to cc32 showed iron-repressed control of fHbp, as regulation was dominated by cotranscription with the iron-repressed upstream gene cbbA. These results highlight the complexity of fHbp regulation and demonstrate that control of transcription can vary among genetic lineages.

Genotypic and phenotypic characterization of carriage and invasive disease isolates of Neisseria meningitidis in Finland

The relationship between carriage and the development of invasive meningococcal disease is not fully understood. We investigated the changes in meningococcal carriage in 892 military recruits in Finland during a nonepidemic period (July 2004 to January 2006) and characterized all of the oropharyngeal meningococcal isolates obtained (n = 215) by using phenotypic (serogrouping and serotyping) and genotypic (porA typing and multilocus sequence typing) methods. For comparison, 84 invasive meningococcal disease strains isolated in Finland between January 2004 and February 2006 were also analyzed. The rate of meningococcal carriage was significantly higher at the end of military service than on arrival (18% versus 2.2%; P < 0.001). Seventy-four percent of serogroupable carriage isolates belonged to serogroup B, and 24% belonged to serogroup Y. Most carriage isolates belonged to the carriage-associated ST-60 clonal complex. However, 21.5% belonged to the hyperinvasive ST-41/44 clonal complex. Isolates belonging to the ST-23 clonal complex were cultured more often from oropharyngeal samples taken during the acute phase of respiratory infection than from samples taken at health examinations at the beginning and end of military service (odds ratio [OR], 6.7; 95% confidence interval [95% CI], 2.7 to 16.4). The ST-32 clonal complex was associated with meningococcal disease (OR, 17.8; 95% CI, 3.8 to 81.2), while the ST-60 clonal complex was associated with carriage (OR, 10.7; 95% CI, 3.3 to 35.2). These findings point to the importance of meningococcal vaccination for military recruits and also to the need for an efficacious vaccine against serogroup B isolates.

Extensive genomic variation within clonal complexes of Neisseria meningitidis

Meningococcal disease is a widely distributed complex disease affecting all age categories. It can cause severe meningitis and septicemia, especially in unvaccinated infants and young children. The causative agent, Neisseria meningitidis (Nm), can be phenotypically and genetically differentiated into serogroups and sequence types (STs) and has a highly dynamic population structure. To obtain a deeper understanding of the epidemiology of Nm, we sequenced seven Nm genomes. Large-scale genomic analysis was conducted with these 7 Nm genomes, 27 additional Nm genomes from GenBank, and 4 other Neisseria genomes. We observed extensive homologous recombination in all gene functional categories among different Nm genomes. Homologous recombination is so frequent that it has resulted in numerous chimeric open reading frames, including genes in the capsule biosynthesis cluster and loci targeted by commercial vaccines. Our results reveal that, despite widespread use, evolutionary relationships inferred from the standard seven-gene multilocus sequence typing (MLST) method could not predict virulence gene content or strain phenotype. In fact, up to 28% of the virulence-associated genes could differ between strains of identical STs. Consistent with previous studies, we found that allelic recombination is also associated with alterations in antibiotic susceptibility. Overall, these findings emphasize the extensive genomic plasticity of Nm and the limitations of standard molecular methods to quantify this genotypic and phenotypic diversity.

Recombination and population structure in Salmonella enterica

Salmonella enterica is a bacterial pathogen that causes enteric fever and gastroenteritis in humans and animals. Although its population structure was long described as clonal, based on high linkage disequilibrium between loci typed by enzyme electrophoresis, recent examination of gene sequences has revealed that recombination plays an important evolutionary role. We sequenced around 10% of the core genome of 114 isolates of enterica using a resequencing microarray. Application of two different analysis methods (Structure and ClonalFrame) to our genomic data allowed us to define five clear lineages within S. enterica subspecies enterica, one of which is five times older than the other four and two thirds of the age of the whole subspecies. We show that some of these lineages display more evidence of recombination than others. We also demonstrate that some level of sexual isolation exists between the lineages, so that recombination has occurred predominantly between members of the same lineage. This pattern of recombination is compatible with expectations from the previously described ecological structuring of the enterica population as well as mechanistic barriers to recombination observed in laboratory experiments. In spite of their relatively low level of genetic differentiation, these lineages might therefore represent incipient species.

Neisseria meningitidis is structured in clades associated with restriction modification systems that modulate homologous recombination

Molecular data on a limited number of chromosomal loci have shown that the population of Neisseria meningitidis (Nm), a deadly human pathogen, is structured in distinct lineages. Given that the Nm population undergoes substantial recombination, the mechanisms resulting in the evolution of these lineages, their persistence in time, and the implications for the pathogenicity of the bacterium are not yet completely understood. Based on whole-genome sequencing, we show that Nm is structured in phylogenetic clades. Through acquisition of specific genes and through insertions and rearrangements, each clade has acquired and remodeled specific genomic tracts, with the potential to impact on the commensal and virulence behavior of Nm. Despite this clear evidence of a structured population, we confirm high rates of detectable recombination throughout the whole Nm chromosome. However, gene conversion events were found to be longer within clades than between clades, suggesting a DNA cleavage mechanism associated with the phylogeny of the species. We identify 22 restriction modification systems, probably acquired by horizontal gene transfer from outside of the species/genus, whose distribution in the different strains coincides with the phylogenetic clade structure. We provide evidence that these clade-associated restriction modification systems generate a differential barrier to DNA exchange consistent with the observed population structure. These findings have general implications for the emergence of lineage structure and virulence in recombining bacterial populations, and they could provide an evolutionary framework for the population biology of a number of other bacterial species that show contradictory population structure and dynamics.

Inference of homologous recombination in bacteria using whole-genome sequences

Bacteria and archaea reproduce clonally, but sporadically import DNA into their chromosomes from other organisms. In many of these events, the imported DNA replaces an homologous segment in the recipient genome. Here we present a new method to reconstruct the history of recombination events that affected a given sample of bacterial genomes. We introduce a mathematical model that represents both the donor and the recipient of each DNA import as an ancestor of the genomes in the sample. The model represents a simplification of the previously described coalescent with gene conversion. We implement a Monte Carlo Markov chain algorithm to perform inference under this model from sequence data alignments and show that inference is feasible for whole-genome alignments through parallelization. Using simulated data, we demonstrate accurate and reliable identification of individual recombination events and global recombination rate parameters. We applied our approach to an alignment of 13 whole genomes from the Bacillus cereus group. We find, as expected from laboratory experiments, that the recombination rate is higher between closely related organisms and also that the genome contains several broad regions of elevated levels of recombination. Application of the method to the genomic data sets that are becoming available should reveal the evolutionary history and private lives of populations of bacteria and archaea. The methods described in this article have been implemented in a computer software package, ClonalOrigin, which is freely available from http://code.google.com/p/clonalorigin/.

Independent evolution of the core and accessory gene sets in the genus Neisseria: insights gained from the genome of Neisseria lactamica isolate 020-06

BACKGROUND

The genus Neisseria contains two important yet very different pathogens, N. meningitidis and N. gonorrhoeae, in addition to non-pathogenic species, of which N. lactamica is the best characterized. Genomic comparisons of these three bacteria will provide insights into the mechanisms and evolution of pathogenesis in this group of organisms, which are applicable to understanding these processes more generally.

RESULTS

Non-pathogenic N. lactamica exhibits very similar population structure and levels of diversity to the meningococcus, whilst gonococci are essentially recent descendents of a single clone. All three species share a common core gene set estimated to comprise around 1190 CDSs, corresponding to about 60% of the genome. However, some of the nucleotide sequence diversity within this core genome is particular to each group, indicating that cross-species recombination is rare in this shared core gene set. Other than the meningococcal cps region, which encodes the polysaccharide capsule, relatively few members of the large accessory gene pool are exclusive to one species group, and cross-species recombination within this accessory genome is frequent.

CONCLUSION

The three Neisseria species groups represent coherent biological and genetic groupings which appear to be maintained by low rates of inter-species horizontal genetic exchange within the core genome. There is extensive evidence for exchange among positively selected genes and the accessory genome and some evidence of hitch-hiking of housekeeping genes with other loci. It is not possible to define a 'pathogenome' for this group of organisms and the disease causing phenotypes are therefore likely to be complex, polygenic, and different among the various disease-associated phenotypes observed.

Impact of recombination on bacterial evolution

Genetic exchange plays a defining role in the evolution of many bacteria. The recent accumulation of nucleotide sequence data from multiple members of diverse bacterial genera has facilitated comparative studies that have revealed many features of this process. Here we focus on genetic exchange that has involved homologous recombination and illustrate how nucleotide sequence data have furthered our understanding of: (i) the frequency of recombination; (ii) the impact of recombination in different parts of the genome; and (iii) patterns of gene flow within bacterial populations. Summarizing the results obtained for a range of bacteria, we survey evidence indicating that the extent and nature of recombination vary widely among microbiological species and often among lineages assigned to the same microbiological species. These results have important implications in studies ranging from epidemiological investigations to examination of the bacterial species problem.

Host association of Campylobacter genotypes transcends geographic variation

Genetic attribution of bacterial genotypes has become a major tool in the investigation of the epidemiology of campylobacteriosis and has implicated retail chicken meat as the major source of human infection in several countries. To investigate the robustness of this approach to the provenance of the reference data sets used, a collection of 742 Campylobacter jejuni and 261 Campylobacter coli isolates obtained from United Kingdom-sourced chicken meat was established and typed by multilocus sequence typing. Comparative analyses of the data with those from other isolates sourced from a variety of host animals and countries were undertaken by genetic attribution, genealogical, and population genetic approaches. The genotypes from the United Kingdom data set were highly diverse, yet structured into sequence types, clonal complexes, and genealogical groups very similar to those seen in chicken isolates from the Netherlands, the United States, and Senegal, but more distinct from isolates obtained from ruminant, swine, and wild bird sources. Assignment analyses consistently grouped isolates from different host animal sources regardless of geographical source; these associations were more robust than geographic associations across isolates from three continents. We conclude that, notwithstanding the high diversity of these pathogens, there is a strong signal of association of multilocus genotypes with particular hosts, which is greater than the geographic signal. These findings are consistent with local and international transmission of host-associated lineages among food animal species and provide a foundation for further improvements in genetic attribution.