<title>Introduction</title> <p>The opportunistic pathogen <italic>Listeria monocytogenes</italic> causes life-threatening infections in animal and in human populations at risk. This facultative intracellular bacterium is widespread in the environment and infections occur through ingestion of contaminated food ##REF##1943998##[1]##,##REF##6401354##[2]##. Although the species <italic>L. monocytogenes</italic> has long been known to be genetically diverse ##REF##2498876##[3]##, with strains showing differences in their virulence potential ##REF##3087749##[4]##–##UREF##0##[7]##, detailed knowledge of strain diversity and evolution is still lacking.</p> <p>Several methods have been used to differentiate <italic>L. monocytogenes</italic> strains ##REF##16687581##[8]##. The <italic>Listeria</italic> serotyping scheme ##UREF##1##[9]## based on somatic (O) and flagellar (H) antigens currently represents a common language for <italic>L. monocytogenes</italic> isolate typing and investigations into the ecological distribution, epidemiology and virulence of strains. Unfortunately, serotyping discriminates only 13 serotypes, many of which are known to represent genetically diverse groups of strains, and only four serotypes (1/2a, 1/2b, 1/2c, and 4b) cause almost all cases of listeriosis in humans ##REF##1943998##[1]##. Given its higher discriminatory power, pulsed-field gel electrophoresis (PFGE) is considered accurate for epidemiological investigations and of help for surveillance and control of listeriosis ##REF##8074531##[10]##,##REF##9761339##[11]##, but fingerprint-based methods such as PFGE or ribotyping ##REF##7883880##[12]## are difficult to standardize. Hence, inter-laboratory comparisons necessitate considerable harmonization ##REF##11322701##[13]##, which limits knowledge at the global scale. In addition, these widely used methods provide only limited information on the phylogenetic relationships among strains, which is a serious limitation to understand the evolution of important phenotypic traits such as virulence. Sequence-based or SNP-based approaches appear as promising tools for strain typing and phylogeny in <italic>L. monocytogenes</italic>
##REF##12574278##[14]##–##REF##18256228##[17]##. Multilocus sequence typing (MLST) ##REF##9501229##[18]##–##REF##11544367##[20]## can accurately define the clonal framework of bacterial species. MLST has been shown to discriminate among <italic>L. monocytogenes</italic> isolates ##REF##12574278##[14]##,##REF##16878601##[21]##,##REF##16077098##[22]##, but has not yet been applied on a large scale, and an overview of the clonal structure of <italic>L. monocytogenes</italic> is currently not available. The molecular factors that determine ecological differences among strains are also poorly understood.</p> <p>One salient feature of the population structure of <italic>L. monocytogenes</italic> is the distinction of three phylogenetic lineages. Initially, two major lineages were distinguished, mainly based on multilocus enzyme electrophoresis and PFGE ##REF##2498876##[3]##,##REF##8074531##[10]##,##REF##7883880##[12]##,##REF##1937753##[23]##,##REF##10879987##[24]##, with a third lineage being subsequently recognized based on virulence gene variation, ribotyping and DNA arrays ##REF##7496516##[25]##–##REF##14742555##[28]##. Lineage I includes isolates of serotypes 4b, 1/2b, 3b, 4d and 4e, whereas lineage II includes serotypes 1/2a, 1/2c, 3a and 3c. Lineage III contains serotypes 4a and 4c, as well as serotype 4b as was recently discovered ##REF##15262937##[27]##. The relative virulence and contribution of the three lineages and their serotypes to food contamination and clinical burden is subject of debate ##REF##2498876##[3]##, ##REF##9199440##[26]##, ##REF##15262937##[27]##, ##REF##11320113##[29]##–##REF##17188773##[32]##. As each lineage is genetically heterogeneous, a precise delineation of <italic>L. monocytogenes</italic> clones is needed to determine which ones mostly contribute to human or animal infection ##REF##17085705##[16]##,##UREF##2##[33]##,##REF##17215339##[34]##, and this knowledge would set a landmark for further studies on the biological characteristics of the clones and the evolution of molecular mechanisms by which they cause disease ##REF##17204401##[35]##.</p> <p>Several virulence genes play an important role in the virulence of <italic>L. monocytogenes</italic> strains ##REF##11432815##[36]##,##REF##16710323##[37]##. Internalin (InlA) is a surface protein that mediates the entry of <italic>L. monocytogenes</italic> into various non-phagocytic human eukaryotic cells expressing its receptor E-cadherin ##REF##1905979##[38]##,##REF##8601315##[39]## and plays a key role in the crossing of the intestinal barrier, enabling the bacterium to reach the host bloodstream ##REF##11387478##[40]##. Almost all isolates causing listeriosis in humans express a full-length functional InlA, whereas isolates expressing a truncated form are frequently found in food items and the environment and are associated with a lower virulence potential ##REF##15143478##[5]##. Currently, the ecological factors that drive the evolution of these apparently attenuated strains are unknown. Evolution of virulence would be best understood by mapping the variation of virulence genes such as <italic>inlA</italic>, onto the phylogenetic framework of the genomes in which they are presently distributed.</p> <p>The aims of this study were to provide a robust phylogenetic framework based on MLST analysis of a highly diverse isolate collection and determine (i) the population structure of <italic>L. monocytogenes</italic>; (ii) the evolutionary origin and stability of serotypes; and (iii) the patterns of variation of the virulence gene <italic>inlA</italic> with respect to the evolution of the core genome.</p>
<title>Materials and Methods</title> <title>Bacterial isolates</title> <p>A total of 360 <italic>Listeria monocytogenes</italic> and four <italic>L. innocua</italic> isolates were selected from the collections of the French National Reference Centre for <italic>Listeria</italic> and the WHO Collaborative Centre for foodborne listeriosis (##SUPPL##1##
Table S1
##). These 360 <italic>L. monocytogenes</italic> isolates were subdivided in three subsets, each being included in order to address specific questions: (i) a diversity subset of 171 isolates, which included representative isolates of the distinct <italic>L. monocytogenes</italic> serotypes, atypical strains from lineage III, isolates that caused major epidemics throughout the world, strains for which the complete genome sequence is available, 75 historical strains collected from 1924 to 1966 and belonging to H.P.R. Seeliger <italic>Listeria</italic> Culture Collection (Würzburg, Germany), isolates from the environment, food or animals, and research strains from several countries used in previous studies involving the Institut Pasteur <italic>Listeria</italic> laboratory (##SUPPL##1##
Table S1
##); (ii) 126 isolates selected from maternal-fetal cases, collected prospectively and exhaustively from 1987 to 2005 (i.e., 5 to 10 epidemiologically non-related isolates randomly selected per year), and which were included to probe the temporal dynamics of clone prevalence (‘MF chronological’ subset in ##SUPPL##1##
Table S1
##); and (iii) 63 isolates from year 2000, including 25 from bacteremia, 20 from central nervous system (CNS) infection, and 18 from maternal-fetal infection, which were included to investigate the possible association of specific clones with clinical forms (subset ‘Human clinical, 2000’ in ##SUPPL##1##
Table S1
##).</p> <p>Isolates were identified as <italic>L. monocytogenes</italic> using API <italic>Listeria</italic> strips (BioMerieux, La Balme Les Grottes, France). Identification was confirmed and subdivided into serotypes by classical serotyping ##UREF##1##[9]##, which distinguishes 13 serotypes, and multiplex PCR ##REF##15297538##[41]##, which groups <italic>L. monocytogenes</italic> isolates into four major groups (IIA, IIB, IIC et IVB) corresponding to groups of serotypes (##SUPPL##1##
Table S1
##).</p> <title>Multilocus Sequence Typing</title> <p>The MLST scheme used to characterize <italic>Listeria</italic> strains is based on the sequence analysis of the following seven housekeeping genes: <italic>acbZ</italic> (ABC transporter), <italic>bglA</italic> (beta-glucosidase), <italic>cat</italic> (catalase), <italic>dapE</italic> (Succinyl diaminopimelate desuccinylase), <italic>dat</italic> (D-amino acid aminotransferase), <italic>ldh</italic> (lactate deshydrogenase), and <italic>lhkA</italic> (histidine kinase). This MLST scheme was adapted from the MLST system proposed by Salcedo and colleagues ##REF##12574278##[14]##, with the following modifications. First, the template for gene <italic>ldh</italic> was extended from 354 to 453 nucleotides, thus improving strain discrimination. Second, gene templates were shortened because the extremities of the previous templates correspond to parts of the PCR primer sequences, thus possibly not corresponding totally to the genomic sequence of the isolates analyzed. Third, we incorporated universal sequencing tails to the PCR primers (##TAB##0##
Table 1
##), which allows to sequence PCR fragments of all genes using only two primers. DNA extraction was performed by the boiling method ##REF##15297538##[41]##. The PCR amplification conditions were as follows: an initial cycle of 94°C for 4 min; 25 amplification cycles, each consisting of 94°C for 30 s, 52°C for 30 s (except for <italic>bglA</italic> which has an annealing temperature of 45°C), and 72°C for 2 min; and a final incubation at 72°C for 10 min. The PCR products were purified by ultrafiltration (Millipore, France) and were sequenced on both strands with Big Dye v.1.1 chemistry on an ABI3730XL sequencer (Applied BioSystems).</p> <title>
<italic>inlA</italic> gene sequencing</title> <p>The 2,400 bp long <italic>inlA</italic> gene was sequenced from 157 isolates (##SUPPL##1##
Table S1
##) representing the clonal diversity of <italic>L. monocytogenes</italic> (see below). DNA extraction was performed with the Wizard® kit (Promega Corporation, USA). The PCR amplification conditions were as follows: an initial cycle at 94°C for 5 min; 35 amplification cycles, each consisting of 94°C for 30 s, 55.2°C for 30 s, and 72°C for 1 min 30; and a final incubation at 72°C for 10 min. We used external primers for amplification and internal primers for sequencing (##TAB##0##
Table 1
##), which was performed as described above.</p> <title>Data analysis</title> <p>For each MLST locus, an allele number was given to each distinct sequence variant, and a distinct sequence type (ST) number was attributed to each distinct combination of alleles at the seven genes. Numbers were initially based on highest frequency for the frequent alleles and STs, and were subsequently incremented arbitrarily. In order to define the relationships among strains at the microevolutionary level, we performed allelic profile-based comparisons using a minimum spanning tree (MST) analysis with the BioNumerics v5.10 software (Applied-Maths, Sint Maartens-Latem, Belgium). MST analysis links profiles so that the sum of the distances (number of distinct alleles between two STs) is minimized ##REF##15292152##[42]##. Strains were grouped into clonal complexes (clonal families), defined as groups of profiles differing by no more than one gene from at least one other profile of the group ##REF##15152204##[19]##. Accordingly, singletons were defined as STs having at least two allelic mismatches with all other STs.</p> <p>Neighbor-joining tree analysis was performed using MEGA v4 ##REF##11751241##[43]## or SplitsTree v4b06 ##REF##16221896##[44]##. Calculations of recombination tests were performed using RDP3 ##REF##15377507##[45]##. Nucleotide diversity indices were calculated using DNAsp v4 ##REF##10089204##[46]##. ClonalFrame analysis ##REF##17151252##[47]## was performed with 50,000 burn-in iterations and 100,000 subsequent iterations.</p> <p>To test for phylogenetic congruence among genes, one strain of all 39 STs with allelic mismatch distance &gt;0.65 was used in order to exclude the expected congruence among genes at small evolutionary scale due to common clonal descent, as proposed previously ##REF##12682154##[48]##. Neighbor-joining trees were generated using PAUP* v4 ##UREF##3##[49]## for each gene individually and for the concatenated sequence of the seven genes. For each gene, the differences in log likelihood (Δ−ln L) were computed using PAUP* between the tree for that gene and the trees constructed using the other genes, with branch lengths optimized ##REF##11136255##[50]##. These differences were compared to those obtained for 200 randomly generated trees.</p> <p>The relative contribution of recombination and mutation on the short term was calculated using software MultiLocus Analyzer (Brisse, unpublished) and the simplest implementation of the clonal diversification method ##REF##10555280##[51]##,##REF##7973728##[52]##. For each pair of allelic profiles that are closely related, the number of nucleotide changes between the alleles that differ is counted. A single nucleotide difference is considered to be likely caused by mutation, whereas more than one mutation in the same gene portion is considered to derive from recombination, as it is considered unlikely that two mutations would occur on the same gene while the other genes remain identical. No correction was made for single nucleotide differences possibly introduced by recombination.</p> <p>We used the linkage model in S<sc>tructure</sc>
##REF##12930761##[53]## to identify groups with distinct allele frequencies ##REF##12930761##[53]##. This procedure assigns a probability of ancestry for each polymorphic nucleotide for a given number of groups, <italic>K</italic>, and also estimates <italic>q</italic>, the combined probability of ancestry from each of the <italic>K</italic> groups for each individual isolate. We chose three groups for this report because repeated analyses (200,000 iterations, following a burn-in period of 100,000 iterations) with <italic>K</italic> between 1 and 10 showed that the model probability increased dramatically between <italic>K</italic> = 2 and <italic>K</italic> = 3 and only slowly thereafter.</p> <p>The population recombination rate was estimated by a composite-likelihood method with LD<sc>hat</sc>
##REF##11901136##[54]##. LD<sc>hat</sc> employs a parametric approach, based on the neutral coalescent, to estimate the scaled parameter 2<italic>N</italic>
<sub>e</sub>
<italic>r</italic> where <italic>N<sub>e</sub></italic> is the effective population size, and <italic>r</italic> is the rate at which recombination events separate adjacent nucleotides. The crossing-over model L was used for the analysis of biallelic sites.</p> <p>We also tested for presence of positively selected sites using the software <sc>omega</sc>M<sc>ap</sc>
##REF##16387887##[55]##. This program applies a coalescence-based Bayesian strategy that co-estimates the rate of synonymous vs. non-synonymous substitions ω and the population recombination rate ρ, thus circumventing the high rate of false positives arising from incongruent phylogenies ##REF##12871927##[56]##. The following prior distributions were used for the analyses: μ, κ and Φ<sub>indel</sub>: improper inverse, ω: inverse with range 0.0001–10, ρ: inverse with range 0.01–10. The variable block model was chosen for both ω and ρ, with block sizes of 10 and 30, respectively. We created 10 subsets of 50 randomly drawn <italic>inlA</italic> sequences each, and analyzed each subset with 50,000 iterations and 10 reorderings. The first 20,000 sequences were discarded as burn-in period.</p> <title>Nucleotide sequences</title> <p>Sequences generated in this study are available at <ext-link ext-link-type="uri" xlink:href="http://www.pasteur.fr/mlst">www.pasteur.fr/mlst</ext-link> for the seven MLST genes. <italic>inlA</italic> sequences have been deposited in GenBank/EMBL/DDBJ databases under the accession numbers FM178779 to FM178796 and FM179771 to FM179785. Alleles of the seven MLST genes were deposited under the accession numbers FM180227 to FM180445.</p>
<title>Results</title> <title>The majority of clinical isolates of <italic>L. monocytogenes</italic> belong to seven distinct clones</title> <p>The seven gene portions, sequenced in the 360 <italic>L. monocytogenes</italic> isolates, harbored a total of 438 polymorphisms (13.3%; range 7.01%–17.7% per gene) consisting in bi-allelic (404 sites), tri-allelic (32 sites) or four-allelic (2 sites) single nucleotide polymorphisms (SNPs). The average nucleotide diversity π was 2.91%, ranging from 1.18% to 5.98% per gene (##TAB##1##
Table 2
##). The GC% observed in all alleles ranged from 36.5% to 43.3%, consistent with the 39% value observed across the entire <italic>L. monocytogenes</italic> EGDe genome ##REF##11679669##[57]##. The 126 resulting allelic profiles (or sequence types, STs) were distributed into twenty-three clonal complexes (CC) and 22 singletons (##FIG##0##
Figure 1
##). Five CCs (CC2 to CC4, CC7 and CC9) consisted of a central prevalent genotype associated with several much less-frequent single locus variants (SLVs). CC1 was slightly more diverse, as its central genotype had two SLVs that themselves were associated to other variants. ST5 stood out among all singletons by its high frequency. Each of these CCs and singletons is likely to have descended from a single ancestral bacterium, i.e. corresponds to a clone. Remarkably, the seven above-mentioned CCs were well demarcated, as they differed by at least four genes out of seven among themselves (with the exception of CC2 and CC3, with three mismatches between one pair of STs) and by at least three mismatches from all other STs (##FIG##0##
Figure 1, inset A
##). Together, these seven clones comprised 58 (47%) STs and 245 (69%) isolates, and included 73% of the 252 recent (after 1987) clinical isolates. Five of these clones belonged to lineage I (see below) and comprised 177 of 203 (87%) isolates of this lineage. Other frequent clones were CC6, CC8 and CC101, together representing 32 (9%) additional isolates. Reference strains of large outbreaks and genome sequencing project strains were mapped on the disclosed MLST diversity (##FIG##0##
Figure 1
##
<bold>; </bold>
##SUPPL##1##
Table S1
##); for example, ST1, ST6 and ST11 include reference strains of epidemic clones I, II and III ##UREF##2##[33]##,##REF##17215339##[34]##, respectively.</p> <p>Remarkably, most isolates within a given clone had the same serotype, or a restricted set of serotypes. CC1 and CC2 were dominated by isolates of serotype 4b, and included all isolates of serotypes 4d and 4e. CC3 comprised a large proportion (19/42, 45%) of isolates of serotype 1/2b, and included all isolates of serotypes 3b, and all but one (ST75) isolates of serotype 7. These results suggest that serotypes 4d and 4e each derived at least twice from 4b ancestors, consistent with previous data ##REF##17085705##[16]##,##REF##14742555##[28]##, whereas isolates of serotypes 3b and 7 (excepted ST75) may be regarded as serotypic variants of serotype 1/2b CC3 isolates. CC4 (serotype 4b), ST5 (1/2b), CC6 (4b), CC7 (1/2a), CC101 (1/2a) and CC102 (4b) were each homogeneous with respect to serotype. Finally, CC9 included all isolates of serotype 1/2c, indicating that this serotype is genetically homogeneous. Notably, the virulent strain EGDe of serotype 1/2a also fell into CC9. EGDe only differs from ST9 (1/2c) by <italic>dapE</italic> (allele dapE-20 instead of dapE-4 in ST9), while it differs from all other 1/2a strains by several genes. CC9 also comprised the only included isolate of serotype 3c.</p> <p>Among isolates from human cases of listeriosis, we sought to determine the possible association between clones and clinical sources of the isolates. To eliminate the possible effect of the temporal variation (see below), we compared <italic>L. monocytogenes</italic> isolates from a single year (year 2000) and the three major clinical presentations in humans: bacteremia (n = 25), CNS infections (n = 20), and maternal-fetal infections (n = 18). These isolates corresponded to 28 STs, distributed into 7 CCs and 13 singletons (##SUPPL##1##
Table S1
##). There was no association of particular CCs or ST with clinical presentation: the 11 STs with more than one isolate were encountered in at least two clinical sources, and isolates from prevalent CCs or STs were equally isolated from the three clinical forms (##SUPPL##1##
Table S1
##).</p> <p>Possible trends in the relative prevalence of CCs over time were investigated based on 126 isolates from maternal-fetal cases of listeriosis, collected from 1987 to 2005 (##SUPPL##1##
Table S1
##). These isolates (##SUPPL##1##
Table S1
##) fell into 43 STs and were grouped into 7 CCs and 14 singletons. Four CCs (CC1 to CC4) and two singletons (ST5 and ST9) comprised more than 10 isolates. Numbers of isolates of each of these clones over the 19 year period showed distinct patterns of temporal dynamics: while CC1 (4b) and ST9 (1/2c) were sampled equally over the entire period, CC3 (1/2b-3b-7) shows a clear decrease in prevalence (16 isolates before 1995, 2 isolates after; Chi2 p&lt;0.001). In contrast, ST5 (1/2b) was isolated only once before 1997 but 12 times in the second period (p = 0.02). Similarly, CC2 (4b) showed an apparent increase in prevalence (2 vs. 9, p = 0.034).</p> <title>Homologous recombination is rare in <italic>Listeria monocytogenes</italic>
</title> <p>Divergence among genotypes appeared to be mainly driven by the progressive accumulation of mutations over time, as strains diverge from their common ancestor (##FIG##0##
Figure 1, inset B
##). Congruence among the seven individual gene phylogenies obtained for the distantly related STs ##REF##12682154##[48]## was statistically significant (p&lt;0.005), as assessed by the likelihood method ##REF##11136255##[50]##. Similarly, the short-term contribution of recombination to genotypic diversity was modest, as <italic>L. monocytogenes</italic> alleles are five times more likely to change by mutation than by recombination (<italic>r/m</italic> = 0.197). In addition, the <italic>r/m</italic> rate for nucleotides was 0.59, indicating that nucleotides are approximately twice more likely to change by mutation than by recombination. As an independent approach, the composite likelihood of <italic>r/m</italic>
##REF##11901136##[54]## on the concatenated sequence of the seven genes was 0.62 for lineage I, 0.47 for lineage II and 0 for lineage III. <italic>r/m</italic> values of some of the observed housekeeping genes exceeded 1, but lacked statistical significance (##TAB##2##
Table 3
##). Consistently, <italic>r/m</italic> was 0.81 as estimated using ClonalFrame ##REF##17151252##[47]##.</p> <p>In order to determine which lineages underwent recombination events that left a detectable footprint in extant strains, the nucleotide polymorphisms within the seven gene fragments were analyzed with <sc>structure</sc>
##REF##12930761##[53]##,##REF##10835412##[58]##, a Bayesian method that attempts to identify the ancestral sources of nucleotides. The ancestry of each isolate can be estimated as the summed probabilities of derivation from each ancestral group over all polymorphic nucleotides. <sc>structure</sc> recognized three clusters of strains within <italic>L. monocytogenes</italic>, which were largely homogeneous in terms of their ancestral sources of polymorphism (##FIG##1##
Figure 2A
##). However, a number of isolates are likely to have a mixed origin (##FIG##1##
Figure 2A
##), and this was confirmed statistically using RDP3 on the concatenated sequences (##SUPPL##2##
Table S2
##).</p> <title>Phylogenetic structure of <italic>L. monocytogenes</italic> and evolutionary origin of serotypes</title> <p>Because recombination events, even if they are rare, can strongly distort phylogenetic reconstruction, we took into account potential recombination events using ClonalFrame (##FIG##2##
Figure 3
##). The majority-rule consensus tree revealed three major branches, which could be equated to the three currently recognized <italic>L. monocytogenes</italic> lineages I, II and III, as deduced from serotyping data and inclusion of reference strains. In particular, strains with serotypes 4b and 1/2b fell into lineage I, serotypes 1/2a and 1/2c were associated with lineage II, whereas serotypes 4a and 4c belonged to lineage III. The neighbor-joining (NJ) method (##FIG##1##
Figure 2
##) retained the three major lineages, which were also consistent with the three major groups revealed by S<sc>tructure</sc>. However, the obtained branching pattern was conspicuously distinct for those isolates that underwent recombination events. The most conspicuous example was isolate CLIP85, which was clearly associated with lineage III (##FIG##2##
Figure 3
##), but not in the NJ tree (##FIG##1##
Figure 2B
##). This difference could be attributed to horizontal transfer of <italic>lhkA</italic> from lineage II into CLIP85, as detected with high probability by ClonalFrame (##FIG##2##
Figure 3E
##). Likewise, strains that were inferred to have mixed ancestries (##FIG##1##
Figure 2A
##) were placed at the tip of relatively longer branches on the NJ tree (##FIG##1##
Figure 2B
##) than on the tree derived from ClonalFrame tree (##FIG##2##
Figure 3
##).</p> <p>One exceptional isolate, CLIP98 (serotype 1/2a) isolated from a human blood infection in Canada, was placed at the tip of a long branch, thus representing an apparent fourth lineage. Individual gene genealogies based on the neighbor-joining method also placed CLIP98 outside the three lineages, except for genes <italic>dat</italic> and <italic>lhkA,</italic> which clearly associated CLIP98 with lineage II (not shown). Close inspection of the sequence alignment showed that a large proportion of nucleotide changes that distinguished CLIP98 from lineage II strains were clustered in a small number of short segments and corresponded to nucleotide bases also observed in <italic>L. innocua</italic> strains.</p> <p>The phylogenetic relationships within lineage I (##FIG##2##
Figure 3B
##) suggest that serotype 4b is monophyletic, since all strains of this serotype formed a unique branch. Differently, serotype 1/2b is paraphyletic, as it was encountered in two distinct branches, one of which is branching off early in the history of lineage I. Within lineage II, serotype 1/2a was paraphyletic, whereas 1/2c was monophyletic (##FIG##2##
Figure 3C
##). Notably, the sequenced strain EGDe (1/2a) appears to branch off just before the evolutionary change from 1/2a to 1/2c.</p> <title>Evolution of <italic>inlA</italic> coding sequences: local recombination and convergence of truncated forms in defined clones</title> <p>The 2,400-nt coding sequence of virulence factor InlA showed 162 (6.7%) polymorphic sites and 33 alleles undergoing a distinctive pattern of evolution (##FIG##3##
Figure 4
##). First, in contrast to housekeeping genes, phylogenetic analysis of <italic>inlA</italic> sequences revealed a conspicuous pattern of intragenic homologous recombination. Visual inspection of the distribution of polymorphic sites in <italic>inlA</italic> revealed a mosaic structure (##FIG##3##
Figure 4
##), with several regions having estimated recombination rates one order of magnitude higher (ρ≥0.03) than for baseline regions (ρ≈0.003). Thus, each <italic>L. monocytogenes inlA</italic> sequence represents a composite assembly of short sequences with distinct evolutionary history, likely the result of multiple horizontal gene transfer events. Notably, in no case did we find fully- or nearly identical <italic>inlA</italic> sequences in unrelated STs, showing that horizontal transfer of entire <italic>inlA</italic> alleles is infrequent or non-existing, and that the entire <italic>inlA</italic> coding sequence is clone-specific. The short-term and long-term impacts of localized recombination in the <italic>inlA</italic> sequence were contrasted. Over the short term, <italic>inlA</italic> sequences clearly evolved more rapidly than housekeeping genes. For example, <italic>inlA</italic> sequences within lineage II evolved more than twice as fast (c = 2.13, r = 18%) as MLST genes. In contrast, over the long term, <italic>inlA</italic> sequence divergence was restricted, as housekeeping genes were on average more divergent between the three major lineages than are <italic>inlA</italic> sequences (e.g., 4.8% and 1.3%, respectively, between lineages I and III strains). This constraint on <italic>inlA</italic> sequence divergence resulted in the lack of phylogenetic demarcation of lineages I and III (##SUPPL##0##
Figure S1
##), contrasting sharply with housekeeping genes-based phylogeny (##FIG##1##
Figures 2
##
<bold> and </bold>
##FIG##2##
3
##). Thus, import of sequence stretches from other clones accelerates diversification of clones, while homogenizing <italic>inlA</italic> sequences among distantly related strains.</p> <p>Second, when the entire length of <italic>inlA</italic> was considered, purifying selection against amino-acid changes appeared more relaxed than for housekeeping genes, with a Ka/Ks ratio for <italic>inlA</italic> (0.094) higher than for the seven concatenated MLST genes (0.039 for the same 157 isolates). However, the distribution of ω (the Bayesian estimate of the rate of synonymous vs. non-synonymous substitutions) along the sequence was heterogeneous, with a highly constrained LRR-region (ω≈0.04) and moderately constrained Ig-like and B-repeats regions, with peak values of 0.5. It is worth mentioning that no single stretch of the molecule displayed a significant signature of positive selection, a result that contradicts previous analyses in which recombination was not incorporated ##REF##17660431##[59]##. All but one amino-acid changes found in the LRR domain were located in repeats 1 to 6 (##FIG##3##
Figure 4
##), suggesting a stronger constraint on repeats 7 to 15, which are more extensively involved in interactions with E-cadherin ##REF##12526809##[60]##.</p> <p>Truncated forms of InlA have been described and associated with reduced virulence ##REF##15143478##[5]##, ##REF##17085705##[16]##, ##REF##17660431##[59]##, ##REF##9632615##[61]##–##REF##15066811##[64]##. We found four distinct <italic>inlA</italic> alleles that had premature stop codons (PMSC) at positions 492, 539, 577 and 685 (##FIG##0##
Figures 1
##
<bold> and </bold>
##FIG##3##
4
##
<bold>; </bold>
##TAB##3##
Table 4
##). Remarkably, three of these four PMSCs occurred in isolates that belonged to CC9 (##FIG##0##
Figure 1
##). In addition, out of eight previously reported <italic>inlA</italic> sequences leading to truncated forms, five were identical or nearly identical (&lt;2 SNPs) to the <italic>inlA</italic> sequences that are specific of clone CC9, strongly suggesting that they first occurred in isolates of clone CC9 as well. The remaining PMSC (at codon 492) was observed in ST121 (1/2a), and three other PMSCs from previous reports were also observed in <italic>inlA</italic> alleles unrelated to those in CC9 (##TAB##3##
Table 4
##).</p>
<title>Discussion</title> <p>The phylogenetic structure of <italic>L. monocytogenes</italic> was investigated to provide a framework for the evolutionary history, epidemiology and virulence of this model pathogen. We based our analysis on an update of the MLST scheme proposed by Salcedo et al. ##REF##12574278##[14]##. Alternative sets of genes were used previously ##REF##14766571##[15]##, ##REF##15528725##[65]##–##REF##14715765##[67]##, but were either not extensively validated ##REF##15528725##[65]##, biased towards high levels of nucleotide diversity between two particular lineages ##REF##15066813##[66]## or based at least in part on virulence genes ##REF##14766571##[15]##,##REF##16077098##[22]##,##REF##14715765##[67]##. Virulence-associated genes generally provide improved discrimination among strains, but may reflect ecological adaptation and selection. In contrast, housekeeping genes are considered more appropriate to obtain an unbiased view of the population structure, as their polymorphisms can be considered nearly neutral and are less subject to horizontal transfer.</p> <title>Clonal structure of <italic>L. monocytogenes</italic>
</title> <p>A majority of <italic>L. monocytogenes</italic> isolates belonged to limited number of major clones. Although these were defined based on allelic profiles, the same groupings were obtained by MST analysis based on nucleotide sequences (not shown), with the exception of the three STs (ST35, ST15 and ST74) inferred to derive by a recombination event that changed more than two SNPs. Hence, no significant loss of information was incurred by collapsing nucleotide sequence information into allelic profiles, as expected given limited amounts of recombination. Interestingly, the major clones were almost always separated from other strains by at least three allelic mismatches, indicating ancient divergence. Hence, a relaxed criterion (e.g., two allelic mismatches) would have little impact on the assignment of <italic>L. monocytogenes</italic> isolates to particular clones. Given the neat delineation of the major clones and the fact that they account for a large proportion of clinical <italic>L. monocytogenes</italic> isolates, we propose that these genetic entities represent reference units for future studies on strain virulence, ecology or epidemiology.</p> <p>For global population biology and international surveillance purposes, a definitive strain typing scheme is greatly needed ##REF##9501229##[18]##. The present MLST data represent a unifying language on clone characterization in <italic>L. monocytogenes</italic> and are freely available for comparison at <ext-link ext-link-type="uri" xlink:href="http://www.pasteur.fr/mlst">www.pasteur.fr/mlst</ext-link>. Other sequence-based strain characterization methods have been developed ##REF##14766571##[15]##–##REF##18256228##[17]##,##REF##17215339##[34]##. Future determination of the relative power of MLST and these methods for discrimination among <italic>L. monocytogenes</italic> strains, as well as establishment of the correspondence among the sequence types they define, is required for optimal communication and will provide <italic>Listeria</italic> specialists with a choice of methods that may be suited for distinct purposes (e.g., fine-scale epidemiology versus long-term population biology).</p> <p>Although our strain collection was not designed to address the important question of ecological or epidemiological differences among <italic>L. monocytogenes</italic> clones, we found different distributions of serotypes and clones among animal, environmental and clinical isolates that are consistent with previous reports ##REF##9199440##[26]##,##REF##15262937##[27]##; for example, serotypes 4b clones mostly included isolates from clinical or food sources, whereas they are rare among isolates from animals (4 out of 21) or the environment (none out of 9). However, further MLST studies should be performed with ecologically representative collections of isolates. Likewise, our initial temporal analysis of maternal-fetal isolates suggests the existence of temporal shifts in prevalence of clones over relatively short periods of time (19 years), but a larger longitudinal survey is needed to provide a clearer picture of the temporal dynamics of these clones. Finally, the disclosed diversity is based largely on isolates from France, and a worldwide collection may therefore reveal additional diversity. However, it is noteworthy that many reference strains from outbreaks in other countries and continents belonged to clones defined by French isolates (##FIG##0##
Figure 1
##). It is also important to remember that <italic>L. monocytogenes</italic> is an environmental saprophyte, which does not need to infect animals for survival and propagation. The diversity of clinical isolates may thus only represent a particular subset of the entire diversity of this species.</p> <title>Evolution of serotypes</title> <p>The evolutionary relationships of serotypes within the major lineages have not been previously defined, although this knowledge is particularly important for correct interpretation of serotyping data. We found that all serotype 4b strains belonged to a unique branch, consistent with early MLEE data ##REF##2498876##[3]## and recent sequence data ##REF##17085705##[16]##. This result indicates that serotype 4b appeared only once during the evolution of lineage I. Hence, the apparent increased potential of serotype 4b strains to cause outbreaks may be explained by genetic characteristics that evolved ancestrally, before the diversification of serotype 4b into several clones that appear highly successful among clinical isolates. The fact that the 4b branch is nested within a larger diversity made of 1/2b strains suggests that serotype 1/2b is more ancestral than 4b, and possibly represents the ancestral serotype of lineage I. Likewise, within lineage II, paraphyletic serotype 1/2a clearly stands as the most likely ancestor of lineage II, with serotype 1/2c having evolved more recently (##FIG##2##
Figure 3C
##). As a consequence, our results contradict the DNA array-derived hypothesis that serotype 1/2c represents an ancestral state ##REF##14742555##[28]##. Notably, strain EGDe (1/2a) appears to branch off from the ancestor of CC9 1/2c strains just before the evolutionary shift from 1/2a to 1/2c. Hence, strain EGDe may represent a genomic state close to the evolutionary link between 1/2a and 1/2c strains. While exhibiting an antigenic structure that retained ancestral characteristics, its full genomic content ##REF##11679669##[57]## is likely to be more related to that of 1/2c strains than to most other 1/2a strains, consistent with DNA array data ##REF##14742555##[28]##,##REF##15297538##[41]##.</p> <p>Groups of distinct serotypes observed in isolates with the same ST or CC reflect relatively recent evolution of antigenic structures. Under our evolutionary scenario (##FIG##0##
Figures 1
##
<bold> and </bold>
##FIG##1##
2
##), serotypes 7 and 3b derived from serotype 1/2b in the branch corresponding to clone 3, whereas serotypes 4d and 4e each derived at least twice from distinct 4b ancestors, as also proposed recently ##REF##17085705##[16]##. Likewise, serotype 3c may be derived from 1/2c strains (although the reverse, from 3c to 1/2c, cannot be excluded). All these evolutionary changes in serotype involve modification of the somatic antigens ##UREF##4##[68]##. In contrast, H (flagellar) antigens appear stable over evolutionary time, as we inferred a single within-lineage change of the H antigen, namely from antigen A to antigen D, corresponding to the evolution of 1/2c [antigenic formula I,II,(III):B,D] from 1/2a [I,II,(III):A,B] ##UREF##4##[68]##. Knowledge and subsequent characterization of the genetic determinants of somatic and flagellar antigenic structures ##UREF##5##[69]##,##REF##3151210##[70]## will provide molecular details on serotype evolution.</p> <title>Phylogenetic structure of <italic>L. monocytogenes</italic> and recombination</title> <p>The distinct phylogenetic methods used herein consistently identified three major lineages of <italic>L. monocytogenes</italic> strains, in agreement with a wide set of previous studies based on alternative markers ##REF##7496516##[25]##–##REF##15262937##[27]##. The only exception was isolate CLIP98, placed at the tip of a long branch with weak relatedness to any lineage. However, we suggest that CLIP98 does not represent the first disclosed member of a fourth <italic>L. monocytogenes</italic> lineage. Instead, the fact that two genes of CLIP98 were typical of lineage II, whereas numerous polymorphisms in the other five genes were shared with <italic>L. innocua</italic>, suggests that CLIP98 could have a ‘hybrid’ genome derived from lineage II, but which has received from <italic>L. innocua</italic> donors enough recombined fragments to now appear poorly affiliated with its ancestral lineage. Because putative recombined segments are very small and often consist of only two or three SNPs, these imports were overlooked by ClonalFrame. As three nucleotides of CLIP98 (between positions <italic>cat</italic> 364–368) were also uniquely shared with lineage III strains, CLIP98 may correspond to a recombinant strain with multiple ancestries, possibly due to an increased capacity for incorporation of foreign DNA. Further genomic characterization of this isolate is needed to clarify these aspects.</p> <p>
<italic>L. monocytogenes</italic> lineages differed from <italic>L. innocua</italic>, their closest relative, by 11% of nucleotide positions on average. This is consistent with the monophyly of <italic>L. monocytogenes</italic> and does not support the hypothesis of a descent of <italic>L. innocua</italic> as a whole from <italic>L. monocytogenes</italic>
##REF##14742555##[28]##. Amounts of diversity within the three lineages were 0.33%, 0.61% and 1.25%, respectively, consistent with data obtained by others ##REF##16077098##[22]##,##REF##15262937##[27]##. In contrast, the nucleotide divergence was 4.99% between lineages I and II, 5.3% between lineages I and III, and 7.57% between lineages II and III. Notably, there was not a single common allele among the three lineages, as well as between them and <italic>L. innocua</italic>, whereas strains within a particular lineage generally shared at least one allele. Hence, consistent with DNA hybridization data ##UREF##6##[71]##, the three major lineages correspond to clearly demarcated sequence clusters that fulfill the separateness criteria and divergence levels used in other bacterial groups to distinguish species ##REF##9336922##[72]##–##REF##17062419##[75]##. However, as noted earlier ##REF##15262937##[27]##, the issue of taxonomic revision of <italic>L. monocytogenes</italic> needs careful evaluation. In particular, improved sampling (especially from diversity-rich sources such as the environment) would be necessary to challenge the neat demarcation among lineages and characterize their ecology.</p> <p>The rate of homologous recombination within bacterial species can differ widely from one species to another ##REF##8506277##[76]##,##REF##16503511##[77]## and has a profound impact on the validity of phylogenetic analyses ##REF##1938920##[78]##, on the evolutionary stability of genotypes ##REF##9770535##[79]##, on biological features such as virulence ##REF##17204401##[35]##,##REF##16689791##[80]## and on interpretation of typing data ##REF##11544367##[20]##. Previous reports have indicated that <italic>L. monocytogenes</italic> undergoes low levels of recombination among housekeeping genes ##REF##2498876##[3]##,##REF##12574278##[14]##. Here, we quantified the relative impact of recombination and mutation at various levels of phylogenetic depth and found similar estimates with independent approaches. To our knowledge, <italic>L. monocytogenes</italic> is one of the bacterial species with the lowest rate of recombination ##REF##11136255##[50]##,##REF##16503511##[77]##,##REF##12754228##[81]##. Highly restricted levels of recombination were disclosed in all three major lineages, contrasting with a previous proposal that rates of recombination vary among lineages ##REF##15066813##[66]##. Full genome sequencing ##REF##11679669##[57]## suggested that competence genes are present in <italic>L. monocytogenes</italic> EDGe (CC9), but the regulatory genes for their expression have not been identified, and a noncoding RNA (RliE) may regulate negatively the <italic>L. monocytogenes</italic> orthologs required for competence in <italic>Bacillus subtilis</italic>
##REF##17259222##[82]##. Yet, evidence that <italic>L. monocytogenes</italic> retained the ability for localized recombination is clearly provided by the <italic>inlA</italic> gene encoding InlA, in agreement with previous reports ##REF##16077098##[22]##,##REF##15262937##[27]##,##REF##17660431##[59]##. For this gene, recombination events may contribute to the acquisition by the recipient strain of a selective advantage, for example by escaping the immune response while retaining the ability to interact with its receptor. Such exchange could possibly occur in the intestinal lumen, in which multiple <italic>L. monocytogenes</italic> strains may coexist. Hence, rapid diversification of <italic>inlA</italic> sequences contrasts with the inferred high stability of clonal backgrounds defined by housekeeping genes. This illustrates how the phylogenetic structure based on MLST genes provides a scaffold, which sheds light onto the evolution of individual genes exposed to selective pressures, such as virulence genes.</p> <title>InlA evolution</title> <p>Currently, the ecological significance of loss of a full-length InlA is not understood, and the clonal background in which these forms evolve have not been defined precisely. Little is known about the ecology of <italic>L. monocytogenes</italic> clones, but a realistic scenario is that different clonal families might be adapted to different niches, and their occurrence as mammalian pathogens may be of limited significance for their evolutionary success in the long term. Among the four alleles of gene <italic>inlA</italic> identified among ST9 isolates, the one corresponding to the non-truncated form (<italic>inlA</italic>-8) can be inferred to be ancestral, as the three <italic>inlA</italic> alleles corresponding to truncated forms (inlA-9, inlA-10 and inlA-16) differed by only one mutation from inlA-8, whereas they differed by two mutations from each other (##TAB##3##
Table 4
##). In addition, <italic>inlA</italic>-8 was also found in strain EGDe (ST35), which was inferred to branch off before diversification of other CC9 members (##FIG##2##
Figure 3
##). It is intriguing that InlA, an important bacterial factor for host colonization, was repeatedly lost by convergent evolution in the genetically homogeneous 1/2c ST9 genotype. Such a pattern can be explained either by a relaxed selective constraint on maintaining InlA function, or by a selective advantage provided by the loss of a functional InlA protein, in the ecological niche occupied by members of ST9. Determining the natural habitat of ST9 may provide clues as to why the expression of a virulence trait may in fact turn out to be disadvantageous in particular environments.</p>
<p>Conceived and designed the experiments: ALM SB. Performed the experiments: MR RL. Analyzed the data: MR TW FH SB. Wrote the paper: MR TW FH ML ALM SB.</p> <p>
<italic>Listeria monocytogenes</italic> is a model organism for cellular microbiology and host–pathogen interaction studies and an important food-borne pathogen widespread in the environment, thus representing an attractive model to study the evolution of virulence. The phylogenetic structure of <italic>L. monocytogenes</italic> was determined by sequencing internal portions of seven housekeeping genes (3,288 nucleotides) in 360 representative isolates. Fifty-eight of the 126 disclosed sequence types were grouped into seven well-demarcated clonal complexes (clones) that comprised almost 75% of clinical isolates. Each clone had a unique or dominant serotype (4b for clones 1, 2 and 4, 1/2b for clones 3 and 5, 1/2a for clone 7, and 1/2c for clone 9), with no association of clones with clinical forms of human listeriosis. Homologous recombination was extremely limited (r/m&lt;1 for nucleotides), implying long-term genetic stability of multilocus genotypes over time. Bayesian analysis based on 438 SNPs recovered the three previously defined lineages, plus one unclassified isolate of mixed ancestry. The phylogenetic distribution of serotypes indicated that serotype 4b evolved once from 1/2b, the likely ancestral serotype of lineage I. Serotype 1/2c derived once from 1/2a, with reference strain EGDe (1/2a) likely representing an intermediate evolutionary state. In contrast to housekeeping genes, the virulence factor internalin (InlA) evolved by localized recombination resulting in a mosaic pattern, with convergent evolution indicative of natural selection towards a truncation of InlA protein. This work provides a reference evolutionary framework for future studies on <italic>L. monocytogenes</italic> epidemiology, ecology, and virulence.</p> <title>Author Summary</title> <p>
<italic>Listeria monocytogenes</italic> is a pathogen transmitted through contaminated food and is responsible for severe infections, including meningitis and abortion in animals and humans. It is known that many distinct strains of this pathogen exist, and that they differ in their virulence and epidemic potential. Unfortunately, there is currently no standard definition of strains and no comprehensive overview of their evolution. To tackle these serious limitations to the control of listeriosis and to improve knowledge of how virulence evolves, we characterized a large collection of isolates with sequence-based genotyping methods. We were thus able to identify precisely the most prevalent clones of <italic>L. monocytogenes</italic>, i.e., groups of isolates that descend from a single ancestral bacterium, which can now be characterized further for diagnostic purposes and determination of their precise ecology and virulence potential. We also determined how these clones evolved from their common ancestor and the evolutionary history by which they acquired their phenotypic characteristics, such as antigenic structures. Finally, we show that some particular strains tend to lose a virulence factor that plays a crucial role in infection in humans. This is a rare example of evolution towards reduced virulence of pathogens, and the discovery of the selective forces behind this phenomenon may have important epidemiological and biological implications.</p>
<title>Supporting Information</title>
<p>We thank L. Diancourt and C. Tran for expert technical support and A. Leclercq and C. Buchrieser for their expertise on <italic>Listeria</italic> strains.</p>
<fig id="ppat-1000146-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.g001</object-id><label>Figure 1</label><caption><title>Minimum spanning tree analysis of 360 <italic>L. monocytogenes</italic> and four <italic>L. innocua</italic> strains based on MLST data.</title><p>Each circle corresponds to a sequence type (ST). Grey zones surround STs that belong to the same clonal complex (CC; 24 CCs are visible in total). ST numbers are given inside the circles and are enlarged for the central genotypes that define the major CCs (e.g., ST9 defines the central genotype of CC9). The three major lineages are highlighted by polygons. Four <italic>L. innocua</italic> sequence types are also represented (black circles). The lines between STs indicate inferred phylogenetic relationships and are represented as bold, plain, discontinuous and light discontinuous depending on the number of allelic mismatches between profiles (1, 2, 3 and 4 or more, respectively); note that discontinuous links are only indicative, as alternative links with equal weight may exist. There were no common alleles between the three major lineages, <italic>L. innocua</italic>, ST161 (CLIP98) and ST157 (CLIP85); they are arbitrarily linked through ST7 by default. Circles and sectors were colored based on serotyping data according to the provided legend; in addition, rare serotypes (3a, 3c, 4d, 4e) are indicated directly on the Figure. Note that for simplicity, the serotype of strains that were serotyped by the PCR method (##SUPPL##1##Table S1##) was equated to the most frequent serotype of each PCR group (e.g., 1/2a for PCR group IIA). STs in which truncated forms of InlA were found are indicated by a black triangle, with the position of the premature stop codons given after letter Δ. The ST of reference genome strains is indicated. The positioning of H7858 (ECII) is based on 6 genes only, as gene <italic>dat</italic> is incomplete. Inset A. Crosslinks corresponding to one or two allelic mismatches are indicated. Note the absence of links among major clonal complexes, indicative of their neat demarcation. Circles were colored by grey levels according to the number of isolates. Inset B. Correlation between the number of allelic mismatches (number of distinct alleles between MLST profiles) and the average number of nucleotide differences at distinct alleles. Note the regular positive trend, which indicates that <italic>L. monocytogenes</italic> genotypes diverge predominantly by a mutational process ##REF##12754228##[81]##. Allelic mismatch values of 7 correspond mostly to inter-lineages comparisons.</p></caption></fig> <fig id="ppat-1000146-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.g002</object-id><label>Figure 2</label><caption><title>Homologous recombination is rare, but distorts phylogenetic reconstruction.</title><p>A) Proportions of ancestry in seven housekeeping genes of <italic>L. monocytogenes</italic> strains from three ancestral populations as inferred by the linkage model of <sc>structure</sc>. This plot shows one vertical line for each isolate in which the proportions of ancestry from the three sources are color-coded. For example, strain CLIP85 was inferred as having mixed ancestry, with approximately 75% of nucleotides originated from lineage III (blue), whereas 25% of them were inferred to have an origin in lineage II (green). A number of strains from lineage I have a small proportion of nucleotides with ancestry in lineage III, while strains of lineage II had some nucleotides from lineages I (red) or III (blue). The case of CLIP98 is particular, as it was inferred as deriving from lineage II by imports mainly from <italic>L. innocua</italic> (see text). B) Neighbor-joining phylogenetic analysis of concatenated housekeeping gene sequences, using the Tamura-Nei+G+I model. The three major <italic>L. monocytogenes</italic> lineages are recognized. Together, they included all strains except CLIP85 and CLIP98. Bootstrap support of lineages is given at corresponding branches. An asterisk (*) marks strains that were recognized by <sc>structure</sc> to contain a fraction of nucleotides imported from another ancestral population, with corresponding branches colored according to the source of the recombined nucleotides. Recombination events detected independently using with RDP3 are numbered from 1 to 7 (referring to ##SUPPL##2##Table S2##), and the involved genes are indicated.</p></caption></fig> <fig id="ppat-1000146-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.g003</object-id><label>Figure 3</label><caption><title>Phylogeny obtained with ClonalFrame and serotype relationships.</title><p>A) 50% consensus tree obtained after 100,000 iterations (after 50,000 burn-in) for the 130 distinct <italic>Listeria</italic> STs. <italic>L. monocytogenes</italic> appears monophyletic, with three distinct lineages and one strain (CLIP98) considered as a fourth lineage. Note the close association of CLIP85 with other lineage III strains. B) Detailed view of the inferred relationships within lineage I. Note the monophyly of serotype 4b. C) Detailed view of the inferred relationships within lineage II. Note that all strains of serotype 1/2c (CC9) are nested inside the diversity of 1/2a, and that serotype 1/2c seems to have evolved from 1/2a just after the split from strain EDGe. D) Detailed view of the inferred relationships within lineage III. E) Events on the branch that separates CLIP85 (ST157) from the rest of lineage III. Note the high number of nucleotide changes in <italic>lhkA</italic> (seventh gene), inferred by clonalFrame to correspond to a single recombination event.</p></caption></fig> <fig id="ppat-1000146-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.g004</object-id><label>Figure 4</label><caption><title>InlA polymorphisms.</title><p>A) Distribution of the polymorphisms along the 2,400 nt of gene <italic>inlA</italic> in the 33 distinct <italic>inlA</italic> alleles encountered. The scale above the graph is in amino-acids (AA). InlA functional domains are represented as distinct blocks: signal peptide (SP), alpha-domain linker, leucine-rich repeats (LRR), Ig-like, B-repeat, Pre-anchor (PA) and cell wall anchor (CWA). Vertical bars above these blocks, correspond to synonymous nucleotide polymorphisms; below these blocks, non-synonymous polymorphisms resulting in amino-acid changes. Note that the LRR domain, especially repeats 7 to 15, is highly conserved. Triangles indicate the position of premature stop codons (PMSC) observed in this study and previous reports; black triangles: PMSCs observed in clonal complex CC9; red: PMSCs in other clones; see ##TAB##3##Table 4## for details. B) Deduced amino-acid polymorphisms in InlA. Lineages in which the <italic>inlA</italic> alleles were found are indicated on the left (1, 2 or 3). Blocks of amino-acids that are identical to the sequence in reference strain EGDe (allele 8) are color-shaded. Note the mosaic pattern, with blocks of polymorphisms shared between distinct groups of alleles when scrolling along the sequence.</p></caption></fig>
<table-wrap id="ppat-1000146-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.t001</object-id><label>Table 1</label><caption><title>PCR and sequencing primers used.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Locus</td><td align="left" rowspan="1" colspan="1">Putative function of gene</td><td align="left" rowspan="1" colspan="1">Forward primer</td><td align="left" rowspan="1" colspan="1">Reverse primer</td><td align="left" rowspan="1" colspan="1">Location<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">Annealing temperature (°C)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>abcZ</italic>
</td><td align="left" rowspan="1" colspan="1">ABC transporter</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTATCGCTGCTGCCACTTTTATCCA</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCTCAAGGTCGCCGTTTAGAG</named-content>
</td><td align="left" rowspan="1" colspan="1">2,828,236 to 2,830,008</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>bglA</italic>
</td><td align="left" rowspan="1" colspan="1">Beta glucosidase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTAGCCGACTTTTTATGGGGTGGAG</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCCGATTAAATACGGTGCGGACATA</named-content>
</td><td align="left" rowspan="1" colspan="1">343,221 to 344,636</td><td align="left" rowspan="1" colspan="1">45</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>cat</italic>
</td><td align="left" rowspan="1" colspan="1">Catalase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTAATTGGCGCATTTTGATAGAGA</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCAGATTGACGATTCCTGCTTTTG</named-content>
</td><td align="left" rowspan="1" colspan="1">2,871,318 to 2,872,784</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>dapE</italic>
</td><td align="left" rowspan="1" colspan="1">Succinyl diaminopimelate desuccinylase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTACGACTAATGGGCATGAAGAACAAG</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCATCGAACTATGGGCATTTTTACC</named-content>
</td><td align="left" rowspan="1" colspan="1">287,853 to 288,992</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>dat</italic>
</td><td align="left" rowspan="1" colspan="1">D-amino acid aminotransferase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTAGAAAGAGAAGATGCCACAGTTGA</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCTGCGTCCATAATACACCATCTTT</named-content>
</td><td align="left" rowspan="1" colspan="1">1,661,588 to 1,662,457</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ldh</italic>
</td><td align="left" rowspan="1" colspan="1">L-lactate dehydrogenase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTAGTATGATTGACATAGATAAAGA</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCTATAAATGTCGTTCATACCAT</named-content>
</td><td align="left" rowspan="1" colspan="1">214,486 to 215,427</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>lhkA</italic>
</td><td align="left" rowspan="1" colspan="1">Histidine kinase</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">GTTTTCCCAGTCACGACGTTGTAAGAATGCCAACGACGAAACC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTGTGAGCGGATAACAATTTCTGGGAAACATCAGCAATAAAC</named-content>
</td><td align="left" rowspan="1" colspan="1">1,538,498 to 1,539,937</td><td align="left" rowspan="1" colspan="1">52</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Sequencing primers for above genes</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Forward: <named-content content-type="gene">GTT TTC CCA GTC ACG ACG TTG T</named-content>
</td><td align="left" rowspan="1" colspan="1">Reverse: <named-content content-type="gene">TTG TGA GCG GAT AAC AAT TTC</named-content>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>inlA</italic>
</td><td align="left" rowspan="1" colspan="1">Internalin</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CGGATGCAGGAGAAAATCC</named-content>
</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CTTTCACACTATCCTCTCC</named-content>
</td><td align="left" rowspan="1" colspan="1">454,534 to 456,936</td><td align="left" rowspan="1" colspan="1">55</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>inlA internal sequencing primers</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">F1: <named-content content-type="gene">GATATAACTCCACTTGGG</named-content>
</td><td align="left" rowspan="1" colspan="1">R1: <named-content content-type="gene">GCTCTAAGTTAGTGAGTGCG</named-content>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>inlA internal sequencing primers</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">F2: <named-content content-type="gene">GTGGACGGCAAAGAAAC</named-content>
</td><td align="left" rowspan="1" colspan="1">R2: <named-content content-type="gene">GAGATGTTGTTACACCGTC</named-content>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="ppat-1000146-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.t002</object-id><label>Table 2</label><caption><title>Polymorphism of seven housekeeping protein-coding genes among <italic>L. monocytogenes</italic> isolates.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene</td><td align="left" rowspan="1" colspan="1">Template size</td><td align="left" rowspan="1" colspan="1">No. alleles</td><td align="left" rowspan="1" colspan="1">No. (%) polymorphic sites</td><td align="left" rowspan="1" colspan="1">Ks</td><td align="left" rowspan="1" colspan="1">Ka</td><td align="left" rowspan="1" colspan="1">Ka/Ks</td><td align="left" rowspan="1" colspan="1">π</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>abcZ</italic>
</td><td align="left" rowspan="1" colspan="1">537</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">51 (9.49)</td><td align="left" rowspan="1" colspan="1">0.09624</td><td align="left" rowspan="1" colspan="1">0.00139</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">0.0208</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>bglA</italic>
</td><td align="left" rowspan="1" colspan="1">399</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">28 (7.01)</td><td align="left" rowspan="1" colspan="1">0.05412</td><td align="left" rowspan="1" colspan="1">0.00058</td><td align="left" rowspan="1" colspan="1">0.0107</td><td align="left" rowspan="1" colspan="1">0.0118</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>cat</italic>
</td><td align="left" rowspan="1" colspan="1">486</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">47 (9.67)</td><td align="left" rowspan="1" colspan="1">0.09882</td><td align="left" rowspan="1" colspan="1">0.00287</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">0.0221</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>dapE</italic>
</td><td align="left" rowspan="1" colspan="1">462</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">82 (17.7)</td><td align="left" rowspan="1" colspan="1">0.15317</td><td align="left" rowspan="1" colspan="1">0.00894</td><td align="left" rowspan="1" colspan="1">0.058</td><td align="left" rowspan="1" colspan="1">0.0358</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>dat</italic>
</td><td align="left" rowspan="1" colspan="1">471</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">70 (14.9)</td><td align="left" rowspan="1" colspan="1">0.31024</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">0.045</td><td align="left" rowspan="1" colspan="1">0.0598</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ldh</italic>
</td><td align="left" rowspan="1" colspan="1">453</td><td align="left" rowspan="1" colspan="1">71</td><td align="left" rowspan="1" colspan="1">79 (17.4)</td><td align="left" rowspan="1" colspan="1">0.10767</td><td align="left" rowspan="1" colspan="1">0.00235</td><td align="left" rowspan="1" colspan="1">0.0218</td><td align="left" rowspan="1" colspan="1">0.0232</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>lhkA</italic>
</td><td align="left" rowspan="1" colspan="1">480</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">81 (16.9)</td><td align="left" rowspan="1" colspan="1">0.16468</td><td align="left" rowspan="1" colspan="1">0.00437</td><td align="left" rowspan="1" colspan="1">0.0265</td><td align="left" rowspan="1" colspan="1">0.0289</td></tr><tr><td align="left" rowspan="1" colspan="1">Concatenate, 353 strains</td><td align="left" rowspan="1" colspan="1">3,288</td><td align="left" rowspan="1" colspan="1">121</td><td align="left" rowspan="1" colspan="1">438 (13.3)</td><td align="left" rowspan="1" colspan="1">0.12885</td><td align="left" rowspan="1" colspan="1">0.0049</td><td align="left" rowspan="1" colspan="1">0.038</td><td align="left" rowspan="1" colspan="1">0.0291</td></tr><tr><td align="left" rowspan="1" colspan="1">Concatenate, Lineage I (199 strains)</td><td align="left" rowspan="1" colspan="1">3,288</td><td align="left" rowspan="1" colspan="1">48</td><td align="left" rowspan="1" colspan="1">53 (1.61)</td><td align="left" rowspan="1" colspan="1">0.0124</td><td align="left" rowspan="1" colspan="1">0.00067</td><td align="left" rowspan="1" colspan="1">0.054</td><td align="left" rowspan="1" colspan="1">0.0033</td></tr><tr><td align="left" rowspan="1" colspan="1">Concatenate, Lineage II (133 strains)</td><td align="left" rowspan="1" colspan="1">3,288</td><td align="left" rowspan="1" colspan="1">61</td><td align="left" rowspan="1" colspan="1">143 (4.35)</td><td align="left" rowspan="1" colspan="1">0.02481</td><td align="left" rowspan="1" colspan="1">0.00076</td><td align="left" rowspan="1" colspan="1">0.0306</td><td align="left" rowspan="1" colspan="1">0.0061</td></tr><tr><td align="left" rowspan="1" colspan="1">Concatenate, Lineage III (19 strains)</td><td align="left" rowspan="1" colspan="1">3,288</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">156 (4.74)</td><td align="left" rowspan="1" colspan="1">0.04896</td><td align="left" rowspan="1" colspan="1">0.00258</td><td align="left" rowspan="1" colspan="1">0.0527</td><td align="left" rowspan="1" colspan="1">0.0125</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="ppat-1000146-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.t003</object-id><label>Table 3</label><caption><title>Comparison of mutation rates (<italic>μ</italic>) and recombination rates (r) per base.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Group</td><td colspan="3" align="left" rowspan="1">Θ</td><td colspan="3" align="left" rowspan="1">
<italic>ρ</italic>
</td><td colspan="3" align="left" rowspan="1">r/<italic>μ</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">II</td><td align="left" rowspan="1" colspan="1">III</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">II</td><td align="left" rowspan="1" colspan="1">III</td><td align="left" rowspan="1" colspan="1">I</td><td align="left" rowspan="1" colspan="1">II</td><td align="left" rowspan="1" colspan="1">III</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">ldh</td><td align="left" rowspan="1" colspan="1">0.00519</td><td align="left" rowspan="1" colspan="1">0.01553</td><td align="left" rowspan="1" colspan="1">0.00629</td><td align="left" rowspan="1" colspan="1">0.15011</td><td align="left" rowspan="1" colspan="1">0.01325 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.00625</td><td align="left" rowspan="1" colspan="1">28.923</td><td align="left" rowspan="1" colspan="1">0.8529 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.99364</td></tr><tr><td align="left" rowspan="1" colspan="1">dapE</td><td align="left" rowspan="1" colspan="1">0.00295</td><td align="left" rowspan="1" colspan="1">0.1342</td><td align="left" rowspan="1" colspan="1">0.1762</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">bglA</td><td align="left" rowspan="1" colspan="1">0.00342</td><td align="left" rowspan="1" colspan="1">0.00790</td><td align="left" rowspan="1" colspan="1">0.00499</td><td align="left" rowspan="1" colspan="1">0.01013 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.01000 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.0296 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.7900 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0 <xref ref-type="table-fn" rid="nt105">**</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">lhkA</td><td align="left" rowspan="1" colspan="1">0.00121</td><td align="left" rowspan="1" colspan="1">0.00121</td><td align="left" rowspan="1" colspan="1">0.04512</td><td align="left" rowspan="1" colspan="1">0.01042</td><td align="left" rowspan="1" colspan="1">0.05208</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">8.609</td><td align="left" rowspan="1" colspan="1">43.044</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">dat</td><td align="left" rowspan="1" colspan="1">0.0165</td><td align="left" rowspan="1" colspan="1">0.00248</td><td align="left" rowspan="1" colspan="1">0.00422</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">abcZ</td><td align="left" rowspan="1" colspan="1">0.00290</td><td align="left" rowspan="1" colspan="1">0.00437</td><td align="left" rowspan="1" colspan="1">0.01179</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.00372 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.00372</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.8523 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.6165</td></tr><tr><td align="left" rowspan="1" colspan="1">cat</td><td align="left" rowspan="1" colspan="1">0.00200</td><td align="left" rowspan="1" colspan="1">0.01273</td><td align="left" rowspan="1" colspan="1">0.00527</td><td align="left" rowspan="1" colspan="1">0.00823</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">4.1152</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">Concat</td><td align="left" rowspan="1" colspan="1">0.00296</td><td align="left" rowspan="1" colspan="1">0.00834</td><td align="left" rowspan="1" colspan="1">0.01352</td><td align="left" rowspan="1" colspan="1">0.00182 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.00395 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.6165 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0.4741 <xref ref-type="table-fn" rid="nt105">**</xref>
</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="ppat-1000146-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1000146.t004</object-id><label>Table 4</label><caption><title>Premature stop codons identified in internalin gene <italic>inlA</italic> of <italic>Listeria monocytogenes</italic> isolates.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">AA position of stop codon</td><td align="left" rowspan="1" colspan="1">AA position of mutation</td><td align="left" rowspan="1" colspan="1">Nt position of mutation</td><td align="left" rowspan="1" colspan="1">
<italic>inlA</italic> allele</td><td align="left" rowspan="1" colspan="1">ST (CC)</td><td align="left" rowspan="1" colspan="1">Isolates or mutation name</td><td align="left" rowspan="1" colspan="1">Reference</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">6 (1 nt del.→frameshift)</td><td align="left" rowspan="1" colspan="1">Id. to inlA-8 except for pos. 6</td><td align="left" rowspan="1" colspan="1">ST9 (CC9)</td><td align="left" rowspan="1" colspan="1">3 strains</td><td align="left" rowspan="1" colspan="1">Orsi et al., 2006</td></tr><tr><td align="left" rowspan="1" colspan="1">460</td><td align="left" rowspan="1" colspan="1">513</td><td align="left" rowspan="1" colspan="1">1,539 (1 nt del.→frameshift)</td><td align="left" rowspan="1" colspan="1">undefined</td><td align="left" rowspan="1" colspan="1">lineage II</td><td align="left" rowspan="1" colspan="1">NV8</td><td align="left" rowspan="1" colspan="1">Rousseaux et al., 2004</td></tr><tr><td align="left" rowspan="1" colspan="1">492</td><td align="left" rowspan="1" colspan="1">492</td><td align="left" rowspan="1" colspan="1">1,474 (C→T)</td><td align="left" rowspan="1" colspan="1">inlA-15</td><td align="left" rowspan="1" colspan="1">ST121</td><td align="left" rowspan="1" colspan="1">LM71322, LM74705</td><td align="left" rowspan="1" colspan="1">This study</td></tr><tr><td align="left" rowspan="1" colspan="1">492</td><td align="left" rowspan="1" colspan="1">492</td><td align="left" rowspan="1" colspan="1">1,474 (C→T)</td><td align="left" rowspan="1" colspan="1">Id. to inlA-7 (over pos. 859–1,591) except pos. 1,474</td><td align="left" rowspan="1" colspan="1">ST20 ?</td><td align="left" rowspan="1" colspan="1">H1</td><td align="left" rowspan="1" colspan="1">Olier et al. 2002; Olier 200et al., 2003</td></tr><tr><td align="left" rowspan="1" colspan="1">519</td><td align="left" rowspan="1" colspan="1">513</td><td align="left" rowspan="1" colspan="1">1,539 (1 nt del.→frameshift)</td><td align="left" rowspan="1" colspan="1">1 SNP and 1 del. to inlA-8 (over pos. 859–1,591)</td><td align="left" rowspan="1" colspan="1">ST9 (CC9)</td><td align="left" rowspan="1" colspan="1">NV7</td><td align="left" rowspan="1" colspan="1">Rousseaux et al., 2004</td></tr><tr><td align="left" rowspan="1" colspan="1">539</td><td align="left" rowspan="1" colspan="1">539</td><td align="left" rowspan="1" colspan="1">1,615 (C→T)</td><td align="left" rowspan="1" colspan="1">inlA-10</td><td align="left" rowspan="1" colspan="1">ST9</td><td align="left" rowspan="1" colspan="1">LM57179</td><td align="left" rowspan="1" colspan="1">This study</td></tr><tr><td align="left" rowspan="1" colspan="1">577</td><td align="left" rowspan="1" colspan="1">545</td><td align="left" rowspan="1" colspan="1">1,636 (1 nt del.→frameshift)</td><td align="left" rowspan="1" colspan="1">inlA-16</td><td align="left" rowspan="1" colspan="1">ST9</td><td align="left" rowspan="1" colspan="1">LM6186, LM70179, LM73771, LM95474</td><td align="left" rowspan="1" colspan="1">This study</td></tr><tr><td align="left" rowspan="1" colspan="1">577</td><td align="left" rowspan="1" colspan="1">546</td><td align="left" rowspan="1" colspan="1">1,637 (1 nt del.→frameshift)</td><td align="left" rowspan="1" colspan="1">undefined</td><td align="left" rowspan="1" colspan="1">ST9 (CC9)</td><td align="left" rowspan="1" colspan="1">LO28</td><td align="left" rowspan="1" colspan="1">Jonquieres et al., 1998</td></tr><tr><td align="left" rowspan="1" colspan="1">606</td><td align="left" rowspan="1" colspan="1">606</td><td align="left" rowspan="1" colspan="1">1,818 (T→A)</td><td align="left" rowspan="1" colspan="1">Id. To inlA-29 (over pos. 1,642 to 2,347) except pos. 1,818)</td><td align="left" rowspan="1" colspan="1">ST5 or close to ST5</td><td align="left" rowspan="1" colspan="1">Mutation type 1</td><td align="left" rowspan="1" colspan="1">Nightingale et al., 2005</td></tr><tr><td align="left" rowspan="1" colspan="1">656</td><td align="left" rowspan="1" colspan="1">656</td><td align="left" rowspan="1" colspan="1">1,966 (C→T)</td><td align="left" rowspan="1" colspan="1">3 SNPs to inlA-29 (over pos. 1,642 to 2,347)</td><td align="left" rowspan="1" colspan="1">ST5 or close to ST5</td><td align="left" rowspan="1" colspan="1">Mutation type 2</td><td align="left" rowspan="1" colspan="1">Nightingale et al., 2005</td></tr><tr><td align="left" rowspan="1" colspan="1">685</td><td align="left" rowspan="1" colspan="1">685</td><td align="left" rowspan="1" colspan="1">2,054 (G→A)</td><td align="left" rowspan="1" colspan="1">inlA-9</td><td align="left" rowspan="1" colspan="1">ST9</td><td align="left" rowspan="1" colspan="1">LM15839, LM77097, LM77571, LM86391, LM93096</td><td align="left" rowspan="1" colspan="1">This study</td></tr><tr><td align="left" rowspan="1" colspan="1">685</td><td align="left" rowspan="1" colspan="1">685</td><td align="left" rowspan="1" colspan="1">2,054 (G→A)</td><td align="left" rowspan="1" colspan="1">Id. to inlA-9 (over pos. 859–2,189)</td><td align="left" rowspan="1" colspan="1">ST9 (CC9)</td><td align="left" rowspan="1" colspan="1">NV5</td><td align="left" rowspan="1" colspan="1">Rousseaux et al., 2004</td></tr><tr><td align="left" rowspan="1" colspan="1">700</td><td align="left" rowspan="1" colspan="1">700</td><td align="left" rowspan="1" colspan="1">2,100 (C→G)</td><td align="left" rowspan="1" colspan="1">Id. to inlA-8 (over pos. 1,642 to 2,347) except pos. 2,100</td><td align="left" rowspan="1" colspan="1">ST9 (CC9)</td><td align="left" rowspan="1" colspan="1">Mutation type 3</td><td align="left" rowspan="1" colspan="1">Nightingale et al., 2005</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="ppat.1000146.s001"><label>Figure S1</label><caption><p>inlA phylogeny. Neighbor-joining tree of the 33 inlA alleles (midpoint rooting). Alleles resulting in truncated forms of the protein InlA are indicated by a delta letter followed by the position of the stop codon. Isolates in which the alleles were observed are indicated in ##SUPPL##1##Table S1##. The lineages in which the alleles were observed are indicated by colored vertical bars on the right. Note the lack of separation of lineages I and III, as opposed to the results obtained with housekeeping genes (##FIG##1##Figures 2## and ##FIG##2##3##).</p><p>(0.25 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="ppat.1000146.s002"><label>Table S1</label><caption><p>Strains</p><p>(0.09 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="ppat.1000146.s003"><label>Table S2</label><caption><p>RDP analysis</p><p>(0.03 MB XLS)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>a</label><p>Positions correspond to complete genome sequence of strain EGDe (NC003210).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><p>Ks: No. of synonymous changes per synonymous site. Ka: No. of non-synonymous changes per non-synonymous site.</p></fn><fn id="nt103"><p>π: nucleotide diversity.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><label>a</label><p>Concat., concatenated data set. Values for rho (ρ) were obtained by dividing the per-locus recombination rate estimate from LDhat by the sequence length.</p></fn><fn id="nt105"><p><bold>**:</bold> Estimates that are significant at the 5% level. Θ and <italic>ρ</italic> correspond to the population estimates of mutation and recombination rates, respectively.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><p>ST: sequence type; CC: clonal complex; pos.: position; Id.: identical; nt: nucleotide; del. Deleted; SNP: single nucleotide polymorphism.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This study received financial support from Institut Pasteur and the Institut de Veille Sanitaire (St-Maurice, France), and the GPH Program “Towards new therapeutics against low GC% Gram-positive bacteria”.</p></fn></fn-group>
<graphic id="ppat-1000146-t001-1" xlink:href="ppat.1000146.t001"/> <graphic xlink:href="ppat.1000146.g001"/> <graphic id="ppat-1000146-t002-2" xlink:href="ppat.1000146.t002"/> <graphic id="ppat-1000146-t003-3" xlink:href="ppat.1000146.t003"/> <graphic xlink:href="ppat.1000146.g002"/> <graphic xlink:href="ppat.1000146.g003"/> <graphic xlink:href="ppat.1000146.g004"/> <graphic id="ppat-1000146-t004-4" xlink:href="ppat.1000146.t004"/>
<media xlink:href="ppat.1000146.s001.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="ppat.1000146.s002.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="ppat.1000146.s003.xls"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["7"], "element-citation": ["\n"], "surname": ["Brosch", "Catimel", "Milon", "Buchrieser", "Vindel"], "given-names": ["R", "B", "G", "C", "E"], "year": ["1993"], "article-title": ["Virulence heterogeneity of "], "italic": ["Listeria monocytogenes"], "source": ["J Food Prot"], "volume": ["56"], "fpage": ["293"], "lpage": ["301"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["Seeliger", "Hohne"], "given-names": ["HP", "K"], "year": ["1979"], "article-title": ["Serotyping of "], "italic": ["Listeria monocytogenes"], "source": ["Methods Microbiol"], "volume": ["13"], "fpage": ["31"], "lpage": ["49"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Kathariou", "Torrence", "Isaacson"], "given-names": ["S", "ME", "RE"], "year": ["2003"], "article-title": ["Foodborne outbreaks of listeriosis and epidemic-associated lineages of "], "italic": ["Listeria monocytogenes"], "source": ["Microbial food safety in animal agriculture"], "publisher-loc": ["Ames, IA"], "publisher-name": ["Iowa State University Press"]}, {"label": ["49"], "element-citation": ["\n"], "surname": ["Swofford"], "given-names": ["DL"], "year": ["1998"], "source": ["PAUP*: phylogenetic analysis using parsimony (and other methods). 4.0b10 ed"], "publisher-loc": ["Sunderland, Massachusetts"], "publisher-name": ["Sinauer"]}, {"label": ["68"], "element-citation": ["\n"], "surname": ["Seeliger", "Jones"], "given-names": ["HPR", "D"], "year": ["1986"], "source": ["Genus "], "italic": ["Listeria"], "publisher-loc": ["Baltimore"], "publisher-name": ["Williams & Wilkins"], "fpage": ["1235"], "lpage": ["1245"]}, {"label": ["69"], "element-citation": ["\n"], "surname": ["Fiedler", "Seger", "Schrettenbrunner", "Seeliger"], "given-names": ["F", "J", "A", "HPR"], "year": ["1984"], "article-title": ["The biochemistry of murein and cell wall teichoic acids in the genus "], "italic": ["Listeria"], "source": ["Syst Appl Microbiol"], "volume": ["5"], "fpage": ["360"], "lpage": ["376"]}, {"label": ["71"], "element-citation": ["\n"], "surname": ["Rocourt", "Grimont", "Grimont", "Seeliger"], "given-names": ["J", "F", "PAD", "HPR"], "year": ["1982"], "article-title": ["DNA relatedness among serovars of "], "italic": ["Listeria monocytogenes"], "source": ["Curr Microbiol"], "volume": ["7"], "fpage": ["383"], "lpage": ["388"]}]
PLoS Pathog. 2008 Sep 5; 4(9):e1000146