<title>Background</title> <p>Positive selection and recombination are two evolutionary forces that are clearly important in the evolution of many microorganisms [##REF##12270895##1##, ####REF##17475002##2##, ##REF##16621913##3##, ##REF##16473049##4##, ##REF##17690297##5##, ##REF##16585510##6##, ##REF##16537397##7##, ##REF##17675366##8##, ##REF##17660431##9####17660431##9##]. A number of studies of natural bacterial populations have found evidence for positive selection in specific genes, including in <italic>Escherichia coli </italic>[##REF##11244580##10##], <italic>Neisseria meningitides </italic>[##REF##12270895##1##,##REF##15020403##11##], and <italic>Listeria monocytogenes </italic>[##REF##16473049##4##,##REF##17660431##9##,##REF##16077098##12##]. Recent whole-genome analyses of <italic>E. coli </italic>[##REF##16585510##6##,##REF##17675366##8##] and <italic>Streptococcus </italic>[##REF##17475002##2##] have also confirmed the importance of positive selection during evolution of these pathogens. One study specifically suggested that, in bacteria, up to 2 × 10<sup>-5 </sup>mutations per genome, per generation, are beneficial [##REF##17690297##5##] and another study reported that more than half of the amino acid substitutions between <italic>E. coli </italic>and <italic>Salmonella enterica </italic>appear to have been fixed by positive selection [##REF##16621913##3##]. Furthermore, gains in fitness associated with nonsynonymous changes have also been confirmed in <italic>in vitro </italic>experiments [##REF##12538876##13##,##REF##15489515##14##]. Lateral gene transfer (LGT), followed by incorporation of homologous DNA into the genome, appears to be common in many bacteria and occurrence of homologous recombination has been described in many microorganisms [##REF##12270895##1##,##REF##17475002##2##,##REF##16473049##4##,##REF##17660431##9##,##REF##15020403##11##,##REF##10331248##15##,##REF##11514442##16##]. Bacterial populations can differ considerably in frequency of recombination though; while some populations appear to be panmictic (e.g., <italic>Helicobacter pylori </italic>[##REF##9770535##17##]), others seem to show much more limited recombination (e.g., <italic>Borrelia burgdorferi </italic>[##REF##8234271##18##]).</p> <p>In absence of recombination, positive selection can be inefficient due to clonal interference and/or genetic load. In the case of "clonal interference", advantageous mutations that arise in different lineages of the same population compete against each other for fixation, which can slow down the fixation of advantageous mutations, and can result in loss of advantageous mutations. "Genetic load" refers to the increase in frequency or fixation in the population of disadvantageous mutations that are linked to advantageous mutations. Recombination not only allows advantageous mutations present in different lineages to be combined and fixed in the same lineage, thus preventing clonal interference [##REF##5980116##19##, ####REF##9720276##20##, ##REF##15342539##21##, ##REF##9888858##22##, ##REF##17713986##23####17713986##23##], but also can break the linkage between the advantageous and disadvantageous mutations, thus counteracting "genetic load" [##REF##15715837##24##, ####REF##11641490##25##, ##REF##16957730##26##, ##REF##12242249##27####12242249##27##]. Positive selection may also play an important role in facilitating maintenance of fragments introduced by recombination in a given population if these fragments confer a selective advantage to the recipient organism.</p> <p>The genus <italic>Listeria </italic>includes both mammalian pathogenic species (i.e., <italic>L. monocytogenes</italic>, a human and animal pathogen and <italic>L. ivanovii</italic>, an animal pathogen) as well as non-pathogenic species (e.g., <italic>L. innocua, L. welshimeri</italic>) [##REF##15709360##28##]. <italic>L. monocytogenes </italic>is a facultative intracellular foodborne pathogen, which can cause severe invasive human disease with case mortality rates of 20% [##REF##10511517##29##]. Adaptive immunity against <italic>L. monocytogenes </italic>is believed to be mainly cellular-mediated [##REF##10899025##30##], although natural antibodies also seem to play a role in protection [##REF##10591647##31##,##REF##10510340##32##]. <italic>L. monocytogenes </italic>also has the ability to grow under a wide range of environmental stress conditions, including temperatures ranging from 0°C to 45°C [##REF##2108109##33##,##REF##4956900##34##], pH ranging from 4 to 9.6 [##REF##12501427##35##,##UREF##0##36##] and salt concentration of up to 10% [##REF##11016615##37##], facilitating its foodborne transmission. <italic>L. monocytogenes </italic>isolates form a structured population with at least four phylogenetic lineages, including lineages I and II, which are common and lineages IIIA/C and IIIB, which are rare [##REF##16077098##12##,##REF##16514149##38##]. Although isolates from all four lineages have been associated with human listeriosis, most human listeriosis cases and outbreaks have been associated with lineage I isolates, in particular those of serotype 4b [##REF##15033265##39##,##REF##11320113##40##]. Lineage II isolates, on the other hand, seem to be overrepresented among isolates from foods and environmental sources, and underrepresented among human clinical cases [##REF##16416906##41##,##REF##15466521##42##]. These findings suggest that lineage I isolates are more virulent than lineage II isolates, which has been supported by a risk assessment [##REF##16496574##43##] as well as by observations that lineage I isolates, on average, show higher measures of tissue culture pathogenicity as compared to lineage II isolates [##REF##15466521##42##,##REF##11157227##44##,##REF##9215449##45##]. In addition, a considerable proportion of lineage II isolates, but only few lineage I isolates, are virulence-attenuated due to nonsense and frameshift mutations in virulence genes resulting in truncated proteins [##REF##17660431##9##,##REF##16332872##46##,##REF##16204519##47##]. Combined, these observations have led to the conclusions that <italic>L. monocytogenes </italic>lineage I may be host-adapted, while lineage II may represent an environmentally-adapted group [##REF##16416906##41##].</p> <p>Interestingly, the pathogenic <italic>L. monocytogenes </italic>is most closely related to the non-pathogenic <italic>L. innocua</italic>. Consequently, the <italic>L. innocua/L. monocytogenes </italic>lineage within the genus <italic>Listeria </italic>has been used as a model system to study the evolution of pathogenicity characteristics, including through comparative genome analyses [##REF##12648839##48##,##REF##11679669##49##]. While gene presence/absence patterns in these two sister species have been probed through both genome sequencing [##REF##11679669##49##] and macroarray [##REF##14742555##50##] studies, facilitating identification of confirmed and putative virulence genes, evolutionary patterns of the core genome of the <italic>L. innocua/L. monocytogenes </italic>lineage have not yet been comprehensively studied. We used genome sequences available for <italic>L. innocua </italic>[##REF##11679669##49##] as well as for two <italic>L. monocytogenes </italic>lineage I and two lineage II strains [##REF##11679669##49##,##REF##15115801##51##] to investigate the contributions of recombination and positive selection to the evolution of the core genome in these <italic>Listeria </italic>lineages and to gain a better understanding of mechanism that may be important in the evolution of core genome genes during diversification of bacterial pathogens.</p>
<title>Methods</title> <title>Genome data</title> <p>Full genome sequence data for four <italic>L. monocytogenes </italic>isolates and one <italic>L. innocua </italic>isolate were used for this study (Table ##TAB##0##1##). Protein and gene sequence data for these five isolates were retrieved from the Comprehensive Microbial Resource [##UREF##1##52##]. The <italic>L. monocytogenes </italic>isolates represented two serotype 4b lineage I strains (F2365, H7858) as well as two serotype 1/2a lineage II strains (EGD-e, F6854). While F2365 shows at least 20 authentic mutations resulting in premature stop codons and demonstrates some atypical invasion characteristics [##REF##17340887##53##], the genome for this strain was included in our analyses to provide us with increased power for our analyses and an appropriate number of sequences to perform the lineage specific analyses for positive selection. Genes with premature stop codons were excluded from the analyses for positive selection and recombination; presence of these genes in F2365 thus did not affect our analyses. To identify orthologous genes found in all five genomes (i.e., genes representing members of the <italic>L. monocytogenes</italic>/<italic>L. innocua </italic>lineage core genome), the predicted protein sequences of each gene from each genome were clustered using BLAST and TribeMCL [##REF##11917018##54##]. Gene clusters were initially identified using TribeMCL (run with the inflation value set at 2) using BLAST cutoff values of 1e-150, followed by identification of clusters containing less conserved genes (using BLAST cutoff values of 1e-100, 1e-50, and 1e-30). This stepwise approach was used to minimize inclusion of multiple genes from the same genome in a given cluster; the majority of the clusters identified had no more than one gene from each genome. Four clusters contained four sequences from each EGD-e, CLIP 11262, F2365, and H7858 as well as two identical or nearly sequences from F6854 (these four sequences were considered paralogs in F6854); only the F6854 sequence that matched the length of the other genes in this cluster was retained. Only clusters containing five sequences, one from each genome, were further analyzed.</p> <p>Orthologs grouped in the same clusters were aligned using the Clustal W method [##REF##7984417##55##]. Alignments were scanned for frameshift mutations, presence of stop codons, and gene sequences with unequal length; in addition, number of informative sites, average nucleotide diversity (π), overall identity, and identity in the first and last 15 nucleotides were obtained for each cluster. Alignments identified as having low identities or containing sequences with different lengths were manually evaluated using the program BioEdit [##UREF##2##56##] and trimmed or otherwise edited if necessary. If alignments contained frameshift mutations generated by indels (insertion/deletion) followed by another indel that restored the original frame, the alignment was edited by removing the region between the frameshift mutations. Final alignments were used for positive selection and recombination analyses as detailed below.</p> <p>Each gene cluster was also assigned to one of 19 COGs (Clusters of Orthologous Groups of proteins) or the category "not in COG" based on the EGD-e genome annotation available in the NCBI genome database. The effective number of codons used in a gene (Nc), a measure of the codon bias, was assessed using the program "chips" implemented in the EMBOSS package [##REF##10827456##57##]. Nucleotide diversity and number of informative sites were obtained from PhiPack outputs (see below under "Recombination analyses"). Genes also were classified as encoding (i) cell wall proteins, (ii) secreted proteins or (iii) membrane proteins based on the classifications listed in the LEGER database [##REF##16381897##58##].</p> <title>Positive selection analysis</title> <p>Genes under positive selection were identified using codeml as implemented in PAML version 3.15 [##REF##9367129##59##]. The models implemented in PAML allow for identification of genes under positive selection (as well as specific sites that are under positive selection in a gene) even if the overall d<sub>N</sub>/d<sub>S </sub>ratio (ω) for a gene is &lt; 1. We employed two types of tests implemented in PAML to identify genes under positive selection. An overall test for positive selection (Test Overall; TO) was carried out using the null model M1a (Nearly-neutral) and the alternative model M2a (Positive selection) [##REF##15514074##60##]; this test identifies genes under positive selection in any or all of the branches of a given phylogeny. To identify genes that are under positive selection in specific branches of the <italic>L. monocytogenes/L. innocua </italic>phylogeny, the branch-site test2 described by Zhang et al. [##REF##16107592##61##] was used. This test was used to identify genes under positive selection in three branches (Fig. ##FIG##0##1##), including (i) the ancestral branch of <italic>L. monocytogenes </italic>lineage I (<underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>I</underline><underline>A</underline>ncestral; TLM1A), (ii) the ancestral branch of <italic>L. monocytogenes </italic>lineage II (<underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>II</underline><underline>A</underline>ncestral; TLM2A), and (iii) the branch separating <italic>L. monocytogenes </italic>and <italic>L. innocua </italic>(<underline>T</underline>est <italic><underline>L</underline>. innocua</italic>/<italic><underline>L</underline>. <underline>m</underline>onocytogenes</italic>; TLI/LM). Because the sequence of the <italic>L. innocua </italic>and <italic>L. monocytogenes </italic>ancestor is unknown, TLI/LM cannot differentiate between positive selection in <italic>L. innocua </italic>and positive selection in the ancestor of <italic>L. monocytogenes</italic>. No test was performed to test for evidence of positive selection among genes within a given lineage. Initially, one universal phylogenetic tree (Fig. ##FIG##0##1##), representing the consensus tree of the 2267 genes analyzed, was used for all PAML analyses. For all genes that were identified as being under positive selection, gene-specific trees were constructed and TO and branch-specific PAML analyses were re-run if the gene-specific tree differed from the consensus tree. For eight genes, gene-specific trees differed from the universal tree. PAML analysis with gene-specific trees confirmed positive selection for five genes, while for three genes analyses with the gene-specific trees did not find any evidence for positive selection (these genes were thus not considered to be under positive selection).</p> <p>For each test, nested models (one null model that does not allow for positive selection and one alternative model that allows for positive selection) were compared using a Likelihood Ratio Test (LRT) as described by Yang et al. [##REF##10790415##62##]. For each model, three replicates were generated and the maximum likelihood values for each model were used in the LRT. Genes with negative LRT values were re-run 10 times and the maximum values for each model were used for LRT. Persistent negative LRT values were rounded to zero (<italic>P </italic>= 1). For all branch-specific tests, one degree of freedom was used to calculate <italic>p</italic>-values, while for the overall test, two degrees of freedom were used to calculate <italic>p</italic>-values.</p> <title>Recombination analyses</title> <p>Four tests were used initially to assess gene clusters for evidence of intragenic recombination, including Sawyer's test implemented in GENECONV version 1.81 [##UREF##3##63##] as well as Neighbor Similarity Score (NSS), Maximum χ<sup>2</sup>, and the Pairwise Homoplasy Index (PHI), the last three implemented in PhiPack [##REF##16489234##64##]. For GENECONV analyses, the parameter g-scale was set to 1; this setting allows polymorphisms within the recombinant fragment, increasing the likelihood of the test to identify ancient recombination events or events where the donated recombinant fragment is similar, but not identical to the recombinant sequence in the alignment. In the GENECONV analyses, only inner fragments were considered. For Maximum χ<sup>2</sup>, a fixed window size of 2/3 the number of polymorphic sites was used. For PHI, a window size of 50 nucleotides was used. <italic>p</italic>-values were estimated using 10,000 permutations of the alignment for GENECONV and 1,000 permutations for NSS, Maximum χ<sup>2 </sup>and PHI. Therefore, for all recombination tests, the <italic>p</italic>-values represent the proportion of test statistics of the permuted alignments that were at least as extreme as the observed test statistic.</p> <p>ClonalFrame version 1.1 [##REF##17151252##65##] was used on selected genes to estimate recombination breakpoints and to help identify the most likely recipients in a given recombination event. ClonalFrame assumes that recombination events generate new polymorphisms in the population and is most useful for data sets where the donor of the recombinant fragment is not present in the data set [##REF##17151252##65##]. Nevertheless, the program can also be used to identify recombination between sequences in the data set, although it might underestimate the amount of recombination in alignments where the donor and recipient are closely related [##REF##17151252##65##]. Analyses were based on two independent runs of the program both using the same settings (100,000 burn-in iteration and data collection for an additional 100,000 iterations; default settings were used for all other parameters). A 95% consensus tree was obtained from these two runs, and only those branches present in the 95% consensus tree were analyzed for recombination.</p> <title>Statistical analyses</title> <p>Correction for multiple testing was performed using the procedure reported by Benjamini &amp; Hochberg [##UREF##4##66##] as implemented in the program Q-Value [##REF##12883005##67##] with the proportion of expected true null hypotheses set to 1 (π<sub>0 </sub>= 1). For each <italic>p</italic>-value, the <italic>q-</italic>value (the expected proportion of false positives among the significant tests) was calculated. Corrections were performed separately for each test (e.g., GENECONV, NSS, etc.; TO, TLI/LM, etc.) to account for testing of multiple genes (i.e., 2237 genes). As the tests used for positive selection are already conservative [##REF##16107592##61##], a false discovery rate (FDR) of 20% was used for the positives selection analyses. For recombination analyses, an FDR of 10% was used to compensate the fact that no correction for multiple tests (Sawyer's test, NSS, Maximum χ<sup>2 </sup>and PHI) was carried out due to the high correlation among the tests.</p> <p>Correlation between COGs and positive selection, recombination, and gene parameters (e.g., gene length, codon bias, nt diversity) were carried out using chi-square tests, Fisher's exact tests, and U-tests implemented in SAS. For association between genes in a given COG category, COGs that were numerically overrepresented among genes under recombination or positive selection were tested for the significance of associations using one-sided tests; Bonferroni corrections were performed based on the number of one-sided tests performed. Significance was set at 5%.</p> <title>Confirmation of positive selection and recombination patterns in selected genes in a larger isolate population</title> <p>To probe whether the positive selection and recombination patterns determined using genome wide analyses on five isolates were representative for larger populations, a diverse set of 40 additional <italic>L. monocytogenes </italic>isolates (Additional file ##SUPPL##0##1##) was assembled and used to determine the sequences for five genes (i.e., <italic>cheA</italic>, <italic>phoP</italic>, <italic>lmo0693</italic>, <italic>flaR </italic>and <italic>lmo2537</italic>) for positive selection and recombination analyses. Isolates were selected to represent the genetic diversity of <italic>L</italic>. <italic>monocytogenes</italic>, including lineage I (19 isolates), II (13 isolates), IIIA/C (5 isolates), and IIIB (3 isolates), as well as diverse sources (e.g., foods, human clinical cases; see Additional file ##SUPPL##0##1##). Experiments with biohazardous materials were approved by the Cornell Institutional Biosafety Committee (MUA #15520). Nucleotide sequences for these five genes have been deposited in GenBank as alignments (PopSet accession numbers 164520363, 164520263, 164520173, 164520083, 164519993).</p> <p>PCR amplification of the five selected genes was carried out using primers and conditions described in Additional file ##SUPPL##1##2##. PCR fragments were purified using Exonuclease I (0.5 U/μl) and Shrimp alkaline phosphatase (0.05 U/μl) (USB, NEB) and sequenced (at the Biotechnology Resource Center, Cornell University) using Big Dye Terminator chemistry and AmpliTaq-FS DNA Polymerase and an automated 3730 DNA Analyzer. Sequences were proofread and aligned using Clustal W using Seqman and Megalign, as implemented in Lasergene 7.2.1. Alignments were used for positive selection and recombination analyses as described above.</p> <title>Swarming assays</title> <p>Swarming assays were performed with six isolates representing each of the six different mutations leading to premature stop codons (i.e., FSL C1-057, FSL F2-649, FSL E1-123, FSL F2-663, FSL S4-766 and FSL F2-086), three isolates bearing full length <italic>flaR </italic>(i.e. 10403S, FSL F2-661 and FSL J1-208), and a non-motile isogenic 10403S Δ<italic>flaA </italic>strain [##REF##16982842##68##], which harbors an in-frame deletion of the gene that encodes the flagellin subunit in <italic>L. monocytogenes</italic>. We furthermore constructed an isogenic in-frame <italic>flaR </italic>null mutant in <italic>L. monocytogenes </italic>10403S background using Splicing-by-Overlap (SOEing) PCR and allelic exchange, as previously described [##REF##2357375##69##] for use in swarming assays.</p> <p>Swarming abilities of <italic>L. monocytogenes </italic>strains were evaluated on semi-soft agar. Strains were initially grown for 24 h at 37°C on BHI agar and colonies were used to stab-inoculate BHI semi-soft agar (0.4%). Swarming ability was assessed by measuring colony area using SigmaScan Pro 5.0 (SPSS Inc., Chicago, IL) for each strain after incubation at room temperature for 48 h. Swarming area for each mutant strain was normalized to the swarming area for strain 10403S, which was set at 100%.</p>
<title>Results</title> <title>Initial identification and characterization of the <italic>L. monocytogenes/L. innocua </italic>core genome</title> <p>Using BLAST and TribeMCL (as detailed in the "Methods"), we identified 2267 orthologous genes present in all five genomes, representing an initial definition of the core genome for the <italic>L. monocytogenes/L. innocua </italic>lineage. The orthologs were highly syntenic in the two <italic>L. monocytogenes </italic>lineages and <italic>L. innocua </italic>(Fig. ##FIG##1##2A##). The 2267 orthologous genes identified represent 76% of the coding genes identified in <italic>L. innocua </italic>CLIP11262 and 80% of the coding genes identified in <italic>L. monocytogenes </italic>EGD-e and F2365 (i.e., the two closed <italic>L. monocytogenes </italic>genomes). Thirty of the 2267 genes in the core genome had ≤ 1 informative site and were thus not used in subsequent analyses; a final set of 2237 genes was thus used in all genome-wide analyses (shown in green in Fig. ##FIG##1##2B##).</p> <p>Genes in COGs "Energy production and conversion", "Amino acid transport and metabolism", "Carbohydrate transport and metabolism", "Replication, recombination and repair", "Defense mechanisms", and "Cell wall/membrane biogenesis" showed a significant tendency to be longer than genes in other COGs (Table ##TAB##1##2##). Although genes categorized into the COGs "Amino acid transport and metabolism", "General functional prediction", "Defense mechanisms", "Coenzyme transport and metabolism", "Replication, recombination and repair", and "Cell wall/membrane biogenesis" showed a significant tendency for a greater number of informative sites than the genes in other COGs (Table ##TAB##1##2##), only genes in the last three COGs showed a significant tendency for a higher average genetic diversity (π)(Table ##TAB##1##2##).</p> <p>Genes categorized into the COGs "Energy production and conversion", "Translation", and "Posttranslational modification, protein turnover, chaperone" showed a significant tendency for a higher codon bias (Table ##TAB##1##2##), possibly reflecting their housekeeping roles and higher expression rates [##REF##6760125##70##, ####REF##15537809##71##, ##REF##17295928##72####17295928##72##]. Conversely, genes in the COGs "Coenzyme transport and metabolism", "Transcription", and "Signal transduction mechanisms" showed a tendency for lower codon bias (Table ##TAB##1##2##), possibly because these genes are not constitutively expressed and are not highly expressed in the cell, lowering the constraint for preferential codons usage [##REF##6760125##70##, ####REF##15537809##71##, ##REF##17295928##72####17295928##72##].</p> <title>A considerable number of <italic>L. monocytogenes </italic>and <italic>L. innocua </italic>genes show evidence for recombination</title> <p>Among the 2237 orthologs tested, 1097 genes (representing approx. 49% of the genes in the <italic>L. monocytogenes/L. innocua </italic>core genome) showed significant evidence (FDR &lt; 10%) for recombination in at least one of the four recombination tests used (these genes are shown in red in Fig. ##FIG##1##2b##). GENECONV, NSS, Maximum χ<sup>2 </sup>and PHI identified 508, 460, 900, and 252 orthologs, respectively, with significant evidence for recombination. A total of 516, 282, 156, and 143 orthologs showed significant evidence for recombination in one, two, three and all four tests, respectively. Genes with evidence for recombination showed a tendency to have longer alignments (<italic>P </italic>&lt; 0.001; One-sided U-Test), lower codon bias (<italic>P </italic>= 0.013), higher nucleotide diversity (<italic>P </italic>&lt; 0.001), and more informative sites (<italic>P </italic>&lt; 0.001) then genes with no evidence for recombination. These findings are consistent with the expectation that, by chance, shorter genes are less likely to be involved in intragenic recombination, but also the observation that shorter sequences provides less power in the analyses for evidence of recombination [##REF##16489234##64##,##REF##11557798##73##].</p> <p>When genes that encode for (i) cell wall proteins, (ii) secreted proteins or (iii) membrane proteins (based on the listings in LEGER) were tested for their prevalence among genes with evidence for recombination, these three genes classes (i.e., cell wall proteins, secreted proteins, and membrane proteins) were under-represented (<italic>P </italic>= 0.002, <italic>P </italic>= 0.001, and <italic>P </italic>= 0.013, respectively; one-sided Fisher's exact test) among the 1097 genes with significant evidence for recombination, suggesting that these genes are less likely to have experienced recombination.</p> <p>Genes in three COGs are overrepresented among the genes that show evidence for a history of recombination (Table ##TAB##2##3##). For example, the "Carbohydrate transport and metabolism" COG was significantly overrepresented among the 1097 genes with evidence for recombination in at least one of the four tests (<italic>P </italic>= 0.012). Genes in this COG were also significantly more likely to have low <italic>p</italic>-values (indicative of significant evidence for recombination) as compared to genes in the other COGs for each of the four tests (as determined by U-tests; see Table ##TAB##2##3##). For three tests (NSS, Maximum χ<sup>2</sup>, and PHI), genes in the "Amino acid transport and metabolism" COG were significantly more likely to have low <italic>p</italic>-values (indicative of significant evidence for recombination), as compared to genes in the other COGs, respectively (Table ##TAB##2##3##). These data may indicate that recombination allows <italic>L. monocytogenes </italic>and <italic>L. innocua </italic>to rapidly generate and acquire diversity in genes involved in carbohydrate and amino acid transport and metabolism, which may facilitate adaptation to environments that differ in nutrient availability (e.g., host and non-host associated environments). As these two COGs also showed a tendency to have longer genes, the association between these COGs and recombination could also be due to an increased power to detect recombination.</p> <title>Thirty-six <italic>L. monocytogenes </italic>and <italic>L. innocua </italic>genes show evidence for positive selection</title> <p>PAML identified 36 genes under positive selection (FDR &lt; 20%) with either the overall test (TO) or the branch specific tests, including one gene (<italic>lmo2178</italic>) identified with two branch specific tests and one gene (<italic>lmo0782</italic>) identified with TO and two branch specific tests (Table ##TAB##3##4##). Three genes were identified as being under positive selection with the overall test (TO). Seven were identified with the <italic>L. innocua </italic>test (TLI/LM), 9 with the <italic>L. monocytogenes </italic>lineage I ancestor test (TLM1A), and 20 with the <italic>L. monocytogenes </italic>lineage II ancestor test (TLM2A). Genes under positive selection showed a tendency to have more informative sites and to be longer than genes not under positive selection (<italic>P </italic>= 0.001 and <italic>P </italic>= 0.031; one-sided U-Test), probably as these two factors may increase the power of the test for positive selection. Positive selection was not associated with nucleotide diversity or codon bias.</p> <p>While eight of the genes under positive selection (i.e., <italic>lmo0098</italic>, <italic>lmo0653</italic>, <italic>lmo0782</italic>. <italic>lmo0785</italic>, <italic>lmo1424</italic>, <italic>lmo1529</italic>, <italic>lmo2215</italic>, and <italic>lmo2596</italic>) encode membrane proteins [##REF##15966022##74##], neither genes that encode for cell wall proteins nor genes that encode for secreted proteins or membrane proteins (based on the listings in LEGER) were significantly overrepresented among the 36 genes under positive selection (one-sided Fisher's exact test). One COG, "Signal transduction mechanisms", had a significant association with positive selection (nominal <italic>P </italic>= 0.008; one-sided Fisher's exact test) though, suggesting an enrichment for genes under positive selection in this category. However, after correction for multiple comparisons, the association is not significant (<italic>P </italic>= 0.098; Bonferroni correction). Because of the low number of genes under positive selection, it was not possible to assess the association between positive selection and most COGs. We thus assessed whether the distribution of the <italic>p</italic>-values for each test deviates from the random distribution for any of the COGs using a U-test. After Bonferroni correction, none of the COGs showed evidence for association with lower <italic>p</italic>-values (indicating evidence for positive selection) for either the lineage I (TLM1A) or the lineage II branch test (TLM2A). While three COGs (i.e., "Cell wall/membrane biogenesis", "Coenzyme transport and metabolism", and "Amino acid transport and metabolism") were associated with lower <italic>p</italic>-values for the TLI/LM test (see Fig. ##FIG##2##3##), only the association for the "Cell wall/membrane biogenesis" COG was significant after Bonferroni correction (nominal <italic>P </italic>= 0.004; Bonferroni corrected <italic>P </italic>= 0.036). Importantly, genes in this COG were not significantly associated with evidence for recombination, suggesting that a significant tendency for these genes to be under positive selection was not driven by an enrichment of genes with a history of recombination.</p> <p>Among the 36 genes that showed evidence for positive selection, 29 genes also showed evidence for recombination, including five genes for which only one of the four recombination tests was significant. Statistical analyses showed that genes with evidence for recombination were overrepresented among the 36 genes found to be under positive selection (chi-square, <italic>P </italic>&lt; 0.001). Among the seven genes with evidence for positive selection and no evidence for recombination, five and two genes showed evidence for positive selection in the lineage II ancestral branch and the <italic>L. monocytogenes/L. innocua </italic>branch, respectively; none of these genes showed evidence for positive selection in the lineage I ancestral branch.</p> <title>Core genome genes encoding MHC antigen do not show evidence for positive selection</title> <p>The six core genome genes encoding antigens known to induce adaptive cellular immunity against <italic>L. monocytogenes </italic>in mice [##REF##7506732##75##, ####REF##8543821##76##, ##REF##8758895##77##, ##REF##9584149##78##, ##REF##8758896##79##, ##REF##1353418##80##, ##REF##7500019##81##, ##REF##9725229##82####9725229##82##] were evaluated for evidence for positive selection. Epitopes in these antigens are presented to CD8<sup>+ </sup>or CD4<sup>+ </sup>T cells through MHC class Ia, Ib or II molecules. One MHC antigen (a putative 23 aa leader peptide encoded by a transcription attenuator upstream of <italic>lmo2165</italic>) showed no nonsynonymous changes and was thus not formally tested for positive selection. The other five MHC antigens (p60, encoded by <italic>iap</italic>; LemA, encoded by <italic>lemA</italic>; a lipoprotein encoded by <italic>lmo1388</italic>; an extracellular solute-binding protein encoded by <italic>lmo0135</italic>; and a protein with unknown function encoded by <italic>lmo1602</italic>) showed no evidence for positive selection. Moreover, the protein alignment of these five antigens showed no amino acid changes in the major epitope regions.</p> <title>Lineage II strains harbor more recombinant fragments than lineage I strains</title> <p>While the four recombination tests detailed above showed that a considerable number of genes in the <italic>L. monocytogenes/L. innocua </italic>core genome had evidence for recombination, these analyses did not allow for easy determination of the recipient strain and thus did not permit us to test the hypothesis that <italic>L. monocytogenes </italic>lineages differ in their frequency of gene fragment acquisition by recombination. ClonalFrame allows for identification of the recipient strains in recombination events and was thus used to analyze a set of 40 randomly selected genes (<italic>clpX</italic>, <italic>lmo0343</italic>, <italic>lmo0405</italic>, <italic>pflC</italic>, <italic>phoP</italic>, <italic>lmo1436</italic>, <italic>lmo1460</italic>, <italic>lmo1537</italic>, <italic>hemC</italic>, <italic>ccpA</italic>, <italic>lmo1623</italic>, <italic>lmo1790</italic>, <italic>lmo2262</italic>, <italic>pepC</italic>, <italic>lmo2391</italic>, <italic>trxB</italic>, <italic>lmo0190</italic>, <italic>lmo0860</italic>, <italic>lmo0877</italic>, <italic>lmo1087</italic>, <italic>proA</italic>, <italic>lmo0992</italic>, <italic>smbA</italic>, <italic>lmo1401</italic>, <italic>lmo1420</italic>, <italic>opuCC</italic>, <italic>trpD</italic>, <italic>lmo1693</italic>, <italic>purK</italic>, <italic>lmo1825</italic>, <italic>panB</italic>, <italic>lmo0028</italic>, <italic>lmo2175</italic>, <italic>lmo2348</italic>, <italic>lmo2566</italic>, <italic>lmo0487</italic>, <italic>lmo0878</italic>, <italic>lmo1004</italic>, <italic>lmo1011</italic>, and <italic>cbiH</italic>) for evidence of recombination and to determine the recipient strains in the recombination events that were identified. Due to computational constraints, testing larger number of genes was not easily feasible. The 40 genes were randomly chosen from a set of 1227 genes that showed no evidence for positive selection, had at least 5 informative sites and had alignment lengths between 600 and 1400 nucleotides. Among the 40 genes selected, 20 showed evidence for recombination in the genome-wide analysis. As recipient lineages cannot be reliably determined for recombination events in the lineage I and II ancestral branches, we only analyzed recombination events in the external branches (Fig. ##FIG##3##4##). Eleven recombination events were identified in the two lineage II strains (five in F6854 and six in EGD-e), while no recombination events were identified in the two lineage I strains (H7858 or F2365); in this data set, lineage I is thus significantly less likely to have recombination events as compared to lineage II (<italic>P </italic>&lt; 0.001; Fisher's exact test).</p> <title>Analyses of five gene sequences obtained for 40 <italic>L. monocytogenes </italic>isolates confirm positive selection and recombination patterns observed in the genome-wide analyses</title> <p>Five genes (Table ##TAB##4##5##) were selected for sequencing in a set of 40 <italic>L. monocytogenes </italic>isolates (representing lineage and source diversity; see Additional file ##SUPPL##0##1##) to confirm the positive selection and recombination patterns observed in the genome-wide analyses. The five genes chosen for these analyses included (i) <italic>cheA </italic>(<italic>lmo0692</italic>), which showed significant evidence for recombination (all tests significant) and positive selection with TLM1A; (ii) <italic>lmo0693</italic>, which showed evidence for positive selection with TLM2A and had no evidence for recombination; (iii) <italic>flaR </italic>(<italic>lmo1412</italic>), which showed evidence for positive selection with TLM2A and no evidence for recombination; (iv) <italic>lmo2537</italic>, which showed significant evidence for recombination with Maximum χ<sup>2 </sup>but no evidence for positive selection; and (v) <italic>phoP </italic>(<italic>lmo2501</italic>), which showed evidence for recombination with GENECONV but showed no evidence for positive selection. Positive selection and recombination analyses were performed with gene alignments containing 45 sequences (40 gene sequences determined here as well as the respective gene sequences from the four <italic>L. monocytogenes </italic>and the one <italic>L. innocua </italic>genome [Table ##TAB##0##1##]).</p> <p>Three genes (<italic>lmo2537</italic>, <italic>phoP</italic>, and <italic>cheA</italic>) showed significant evidence for recombination with all four recombination analyses (i.e., GENECONV, NSS, Maximum χ<sup>2</sup>, and PHI), consistent with the genome-wide analyses, which also found evidence for recombination in these genes. While <italic>flaR </italic>showed no evidence for recombination in the genome-wide analyses, analyses of the 45 <italic>flaR </italic>sequences showed significant evidence for recombination with Maximum χ<sup>2</sup>, PHI, and GENECONV (Table ##TAB##4##5##); these findings suggest that at least some of the additional <italic>flaR </italic>sequences represent recombinant alleles. Based on the 45 sequences, <italic>lmo0693 </italic>showed marginally significant evidence for recombination with NSS (P = 0.033) and no evidence for recombination with the other three methods, largely consistent with the genome-wide recombination analyses, which found no evidence for recombination in this gene.</p> <p>ClonalFrame analysis on a concatenated alignment of the five genes (for all 45 isolates) identified 7, 2, 2, and 3 recombination events in the external branches for <italic>cheA, phoP, flaR</italic>, and <italic>lmo2537</italic>. No recombination events were identified in external branches for <italic>lmo0693</italic>. The recombination events identified in branches other than the ancestral branches leading to a given lineage involved 1, 7, 3, and 3 lineage I, II, IIIA/C, and IIIB branches as recipients in recombination events, respectively. Lineage II strains were significantly more likely to be involved as recipients in recombination events as compared to lineage I strains (<italic>P </italic>= 0.013, Fisher's exact test). ClonalFrame also allowed us to estimate the relative rate of recombination over mutation for these five genes. While, on average, mutations are 4.5 times more common than recombination events (95% IC = {3.0; 7.3}), a recombination event is 1.9 times more likely to change a single nucleotide than a point mutation (95% IC = {1.27; 2.6}).</p> <p>Analyses of positive selection on the 45 sequences for the five genes yielded results similar to those obtained in the genome-wide analyses. As in the genome-wide analyses, we found no evidence for positive selection in the 45 sequences for <italic>lmo2537 </italic>and <italic>phoP</italic>. For <italic>lmo0693 </italic>and <italic>flaR</italic>, the lineage II ancestral branch was identified as evolving by positive selection (<italic>P </italic>&lt; 0.001 and <italic>P </italic>= 0.002, respectively; the same branches were identified as being under positive selection in the genome-wide analyses). In <italic>lmo0693</italic>, two adjacent aa sites (17 and 18) were identified as being under positive selection, including one site identified with a posterior probability &gt; 95% (Table ##TAB##4##5##). In <italic>flaR</italic>, three amino acid sites were identified as evolving by positive selection in the lineage II ancestral branch, including one site identified with a posterior probability &gt; 95% (Table ##TAB##4##5##). While the genome-wide analysis for <italic>cheA </italic>found evidence for positive selection in the lineage I ancestral branch, analysis of the 45 sequences did not find evidence for positive selection in this branch (<italic>P </italic>= 0.207), but found significant evidence (<italic>P </italic>&lt; 0.001) for positive selection in the ancestral branch that leads to lineages I and IIIA/C isolates (branch B, Fig. ##FIG##4##5##), a branch that was not present in the tree used in the genome-wide analyses (due to the absence of lineage IIIA/C strains in that data set). A single amino acid site was identified as evolving by positive selection in <italic>cheA </italic>(posterior probability &gt; 95%; Table ##TAB##4##5##); isolates in lineage I and IIIA/C have glutamine at this site, while all other isolates, including <italic>L. innocua </italic>CLIP11262, bear an alanine. This amino acid site lies within the P2 domain of the CheA protein, which is the binding site for the response regulator, CheY; <italic>cheY </italic>itself seems to be extremely conserved (with no nonsynonymous change among the five strains included in the genome-wide analyses).</p> <p>Interestingly, all three aa sites identified in <italic>flaR</italic>, <italic>lmo0693</italic>, and <italic>cheA </italic>as being under positive selection with posterior probabilities &gt; 95%, involved nucleotide substitutions at all three codon positions (e.g., <italic>cheA </italic>aa site 140 GCC → CAG). These findings suggest that these codons evolved rapidly in the branches under positive selection. The nucleotide diversity in synonymous sites ranged from 0.141 to 0.151 among the three genes with evidence for positive selection while the nucleotide diversity in the synonymous sites of <italic>lmo2537 </italic>and <italic>phoP</italic>, which showed no significant evidence for positive selection, were 0.338 and 0.449, respectively. This suggests that significant tests for positive selection were not due to higher rates of synonymous substitutions, which could result in an underestimation of the number of synonymous changes due to recurrent mutations in synonymous sites.</p> <title>Unrelated <italic>L. monocytogenes </italic>isolates carry independently acquired premature stop codons in flaR</title> <p>Among the 45 <italic>L. monocytogenes flaR </italic>sequences, we identified six sequences with distinct premature stop codons in <italic>flaR </italic>(Fig. ##FIG##4##5##), including five caused by frameshift mutations and one due to a nonsense mutation. The mutations leading to the premature stop codon were found in different regions of the sequence (Additional file ##SUPPL##2##3##). While one of these frameshift mutations was found in two isolates (FSL S4-766 and FSL N4-015, two closely related lineage II strains [Fig. ##FIG##4##5##]), all other mutations were only identified in one isolate each. <italic>flaR </italic>premature stop codons were found in different <italic>L. monocytogenes </italic>lineages, including lineage I (one frameshift and the nonsense mutation), lineage II (three different frameshift mutations), and lineage IIIB (one frameshift mutation) isolates. The isolates carrying premature stop codons in <italic>flaR </italic>were from human clinical cases (n = 4), animals (n = 2) and natural environment (n = 1). No statistical association could be identified between the presence of premature stop codon and source of the isolate (Fisher's exact test).</p> <p>Since <italic>flaR </italic>has been previously shown to be involved in <italic>L. monocytogenes </italic>motility [##REF##8825084##83##], swarming experiments were conducted to assess the swarming ability of six isolates presenting each of the six different mutations leading to premature stop codons, three isolates bearing full length <italic>flaR </italic>and an isogenic <italic>flaR </italic>null mutant. Isolates harboring naturally occurring premature stop codons in <italic>flaR </italic>showed, on average, significantly reduced swarming areas (<italic>P </italic>&lt; 0.001; One-sided U-Test) and FSL F2-649, which harbor a premature stop codon in <italic>flaR </italic>due to a nonsense mutation, had the smallest swarming area among all natural isolates (<italic>P </italic>&lt; 0.001; One-sided U-Test). Even though all isolates with <italic>flaR </italic>premature stop codons still showed some swarming and the 10403S Δ<italic>flaR </italic>showed no significant reduction in swarming ability, these findings suggest that the frameshift mutations in <italic>flaR </italic>may affect motility, at least in some strains.</p>
<title>Discussion</title> <p>We have chosen two closely related species within the genus <italic>Listeria</italic>, i.e., the pathogen <italic>L. monocytogenes </italic>and the non-pathogenic <italic>L. innocua</italic>, as a model system to further probe the evolution of bacterial pathogens using a comparative genomics approach. While comparative genomics approaches have provided important data on the importance of gene acquisitions in the evolution of bacterial pathogens in general [##REF##17475002##2##,##REF##16968231##84##,##REF##12620865##85##] and in the evolution of <italic>L. monocytogenes </italic>in particular (e.g., by identifying a number of virulence genes associated with pathogenic <italic>Listeria </italic>spp.), our understanding of the contributions of recombination and positive selection to the evolution of the <italic>L. monocytogenes/L. innocua </italic>core genome has been limited so far.</p> <title>Recombination and positive selection both contribute to evolution of <italic>L. monocytogenes</italic>, but the relative contribution of these evolutionary forces differs among <italic>L. monocytogenes </italic>lineages</title> <p>Overall, almost half of the genes present in <italic>L. monocytogenes</italic>/<italic>L. innocua </italic>core genome showed evidence for intragenic recombination (location of these genes is shown in Fig. ##FIG##1##2##). Only a much smaller proportion of the genes in the core genome (1.7%) showed evidence for positive selection. By comparison, a recent genomic study on the genus <italic>Streptococcus </italic>reported evidence for recombination in 18 to 37% of the genes in the core genome and evidence for positive selection in 11 to 34% of the genes in the core genome. This study was able to use a larger number of genome sequences (i.e., 26 genomes) and used different criteria for selection of genes with significant evidence for positive selection and recombination (i.e. only genes significant by all three tests for recombination were considered as recombinants) though [##REF##17475002##2##]. While a comparison of the frequency of genes with evidence of recombination between <italic>Listeria </italic>and <italic>Helicobacter pylori</italic>, which appears to be panmictic [##REF##9770535##17##], will be of interest, no genome wide studies of recombination in <italic>Helicobacter </italic>have been reported to date. In both gene-specific and multilocus sequence typing (MLST) studies, evidence has previously been found for contributions of recombination to the evolution of different <italic>L. monocytogenes </italic>genes, even though many of the genes identified as having an apparent history of recombination appear to be lineage and species specific, including a number of genes in the internalin family [##REF##16473049##4##,##REF##17660431##9##] and the <italic>prfA </italic>virulence gene cluster [##REF##16077098##12##]. While some studies also have identified housekeeping genes with significant evidence for recombination [##REF##16077098##12##], the magnitude of recombination in <italic>Listeria </italic>on a genome-wide level has not previously been known. Interestingly, it seems that although homologous recombination between closely related strains of <italic>L. monocytogenes </italic>is common, non-homologous recombination seems to be rare given the high synteny of the different <italic>Listeria </italic>genomes [[##REF##15115801##51##], this study] and the relatively small number of strain and lineage specific genes (e.g., [##REF##8825084##83##] and [##REF##15115801##51##] genes specific to serotype 1/2a and 4b strains, respectively [##REF##15115801##51##]). While others have also previously observed positive selection in <italic>L. monocytogenes</italic>, all previously identified genes under positive selection were virulence genes (i.e., <italic>actA</italic>, <italic>inlA</italic>, <italic>inlB</italic>, <italic>inlC</italic>, <italic>inlC2 </italic>and <italic>inlF</italic>) [##REF##16473049##4##,##REF##17660431##9##,##REF##16077098##12##], which were specific to <italic>L. monocytogenes </italic>or selected <italic>L. monocytogenes </italic>lineages. We are thus the first to identify genes in the core genome which are under positive selection, indicating that positive selection does not just act on accessory and virulence genes.</p> <p>Overall, most genes in the core genome found to be under positive selection also showed evidence for recombination; an association of recombination and positive selection was also previously observed in a genome-wide study on the evolution of genes in the genus <italic>Streptococcus </italic>[##REF##17475002##2##]. Occurrence of both positive selection and recombination in specific genes has also previously been reported for selected <italic>L. monocytogenes </italic>virulence genes [##REF##16473049##4##,##REF##17660431##9##] as well as in other microorganisms [##REF##12270895##1##,##REF##17475002##2##,##REF##11244580##10##,##REF##15020403##11##,##REF##11514442##16##].</p> <p>While it has been shown that high recombination rates can lead to false positives in overall analysis using PAML, such as the TO test [##REF##12871927##86##], no studies have evaluated the effect of recombination on the branch-site tests and it is unclear how recombination affects this analysis. The fact that, in our analyses, more than 1000 genes showed evidence for recombination but no evidence for positive selection suggests that the positively selected genes identified harbor distinct features that allowed their identification. As previously pointed out [##REF##17475002##2##,##REF##17660431##9##], some horizontally transferred fragments may also be more likely to be under positive selection, and positive selection may be important for fixation of recombinant gene and/or for adaptation of the newly acquired genes or alleles to a different function.</p> <p>In analyses of sequence data for 40 genes, lineage II strains were significantly more likely to be identified as recipients of DNA fragments by horizontal gene transfer, indicating either more frequent lateral gene transfer or more frequent fixation of recombinant fragments. Higher frequency of recombinant fragments among lineage II strains has also previously been described for <italic>L. monocytogenes </italic>specific virulence genes (e.g., internalin genes [##REF##16473049##4##,##REF##17660431##9##]) and for some housekeeping genes [##REF##15066813##87##]. The genetic basis of the promiscuity of lineage II isolates is currently unknown but it might be related to the observation that lineage II isolates appear to be overrepresented in foods, farms and natural environments [##REF##16416906##41##], where <italic>Listeria </italic>bacteriophages (listeriophages) may be common due to the frequent presence of <italic>L. monocytogenes </italic>in some of these environments [##REF##15294773##88##,##REF##15553634##89##]. While listeriophages can perform generalized transduction [##REF##10652092##90##], transduction between isolates of serotypes 1/2a and 4b has not been shown. Most listeriophages appear to be serotype-specific and phage host-specificity could account for the differences in recombination frequency between lineages. Alternatively, lineage I strains might have a more effective restriction systems for degradation of foreign DNA, a more efficient mismatch repair system that avoids incorporation of foreign DNA fragments into the chromosome, or may show reduced competency.</p> <p>Our genome-wide analyses also indicated that positive selection occurred in more genes in the ancestral branch of lineage II as compared to the ancestral branch of lineage I. Interestingly, one gene with evidence for positive selection in the ancestral branch of lineage II encodes a putative transcriptional antiterminator of the BglG family (<italic>lmo0297</italic>). As another antiterminator of the BglG family (BvrA), which is only present in lineage II isolates [##REF##14742555##50##], has been shown to be involved in cellobiose-dependent repression of PrfA-dependent virulence genes in <italic>L. monocytogenes </italic>[##REF##10438775##91##], <italic>lmo0297 </italic>may have evolved to facilitate specific transcriptional repression functions in lineage II. Positive selection in the lineage II ancestral branch was further confirmed for two genes (<italic>flaR</italic>, <italic>lmo0693</italic>) in a larger isolate set. As lineage II strains are found in many different environments, including natural environment, farms, foods, animals with clinical disease, as well as human clinical cases (although less frequently than lineage I isolates) [##REF##16416906##41##], one could hypothesize that lineage II isolates are exposed to a more diverse repertoire of distinct selective pressures than lineage I strains, which appear to be less common in natural environments and foods and seem to be more likely to be adapted to human or mammalian hosts [##REF##16416906##41##,##REF##15466521##42##].</p> <title>Diversification, by multiple mechanisms, of cell wall/membrane biogenesis and motility-related genes may play a particularly important role in the evolution of <italic>L. monocytogenes</italic></title> <p>In our study, genes involved in cell wall/membrane biogenesis showed a significant tendency to be identified as being under positive selection in the <italic>L. monocytogenes</italic>/<italic>L. innocua </italic>branch, suggesting that positive selection in these genes contributed to the divergence of these two species. <italic>L. monocytogenes </italic>genes encoding proteins involved in cell wall metabolism and encoding cell wall-anchored proteins have also previously been identified as harboring more nonsynonymous changes than genes involved in other functions [##REF##15115801##51##], further supporting an important role for diversification of these gene categories in the evolution of <italic>L. monocytogenes</italic>. Specific surface associated proteins identified as being under positive selection included genes involved in transport of carbohydrates such as <italic>mptD</italic>, which appears to be involved in resistance to class IIa bacteriocins [##REF##11739758##92##,##REF##12177330##93##]. Class IIa bacteriocins are antimicrobial peptides frequently produced by lactic acid bacteria in foods [##REF##12177330##93##], and resistance to these compounds are likely to confer an advantage to <italic>L. monocytogenes </italic>isolates in foods and environments. In addition, a putative glycosyl transferase (<italic>lmo2121</italic>) was identified as having evolved by positive selection in lineage I; glycosyl transferases could be associated with the differences in the somatic antigens in different serotypes and might be associated with strain differences in virulence and immunogenicity [##REF##15115801##51##,##REF##11015420##94##, ####REF##11029438##95##, ##REF##9882654##96####9882654##96##]. Functionally, active and rapid evolution of genes involved in cell wall/membrane biogenesis is likely to be important to allow bacteria to adapt to different and possibly rapidly changing environments, including, but not limited to, competing microorganisms as well as innate and adaptive immune system effectors in different host species. Our findings are consistent with a recent genome-wide study, which showed that most <italic>E. coli </italic>proteins that undergo positive selection are exposed on the cell surface, including a number of proteins known to interact with bacterial, host, or phage surface molecules [##REF##17675366##8##].</p> <p>A number of genes involved in flagellar synthesis and motility (e. g., <italic>flaR</italic>, <italic>lmo0693 </italic>and <italic>cheA</italic>) also showed evidence for positive selection, including when these genes were analyzed in a larger set of more diverse isolates. <italic>lmo0693</italic>, <italic>cheA </italic>and several other genes involved in motility and chemotaxis have also been found to be expressed more highly in four lineage I strains as compared to two lineage II strains [##UREF##5##97##], further suggesting a difference in regulation of motility and chemotaxis functions between the two lineages. Motility and flagellar expression appear to contribute to both biofilm formation [##REF##17416647##98##] and host cells invasion [##REF##16982842##68##,##REF##16113269##99##,##REF##15155625##100##] in <italic>L. monocytogenes</italic>, and diversification in these genes is likely to be important for adaptation to different host or non-host environments. Interestingly, seven isolates showed premature stop codons in <italic>flaR</italic>, including six due to frameshift mutations and one due to a nonsense mutation. While FlaR has been reported to be a histone-like protein that regulates transcription of the flagellin gene <italic>flaA </italic>in <italic>L. monocytogenes </italic>serotype 1/2c lineage II strain LO28 (as determined in a transposon mutant; [##REF##8825084##83##]), a non-polar <italic>flaR </italic>null mutant in a serotype 1/2a lineage II strain (generated in our study reported here) did not show reduced motility at room temperature, even though the LO28 <italic>flaR </italic>mutant was reported to be non-motile at this temperature [##REF##8825084##83##]. As shown here, natural isolates with premature stop codons in <italic>flaR </italic>showed, though, on average, significant reduced ability to swarm as compared to isolates harboring a full-length gene. <italic>flaR </italic>thus seems to have strain or perhaps serotype or lineage specific functions in <italic>L. monocytogenes </italic>and <italic>flaR </italic>inactivation in some strains is likely to be recent since most frameshift/nonsense mutations are isolate and not clade-specific. Alternatively, some or all frameshift and nonsense mutations could be reversible and <italic>flaR </italic>might be phase-variable, a tempting hypothesis as phase-variable flagella-related genes have been described in a number of other bacteria [##REF##17505524##101##, ####REF##10899861##102##, ##REF##15066026##103##, ##REF##2166334##104####2166334##104##]. Interestingly, previous studies have shown that <italic>inlA</italic>, which encodes another <italic>L. monocytogenes </italic>surface protein that promotes mammalian host cell invasion, also carries several different premature stop codons in both lineage I and II strains [##REF##17660431##9##,##REF##16332872##46##,##REF##9632615##105##, ####REF##12055305##106##, ##REF##12595435##107##, ##REF##15066811##108####15066811##108##]. The isolates with truncated InlA seem to be more common in foods [##REF##16332872##46##] and show significantly reduced invasiveness for human intestinal epithelial cells [##REF##16332872##46##], suggesting reduced virulence of these strains. Diversification by both positive selection and gene inactivation of genes encoding surface molecules with a role in virulence thus appears to be broadly important in the adaptation of <italic>L. monocytogenes </italic>to host and non-host associated environments.</p>
<title>Conclusion</title> <p>Our analyses reported here indicate that both recombination and positive selection contribute to the evolution of the <italic>L. monocytogenes/L. innocua </italic>core genome. While considerably more genes appear to be affected by recombination, positive selection still appears to play an important role in the evolution of both genes in the core genome (this study) as well as <italic>L. monocytogenes </italic>virulence genes that are not part of the core genome [##REF##16473049##4##,##REF##17660431##9##,##REF##16077098##12##]. The list of genes identified as being under positive selection hopefully can be used by the scientific community to advance the discovery of genetic factors that allow this organism to adapt to diverse environments and hosts. In particular, our data suggest important roles for positive selection and diversification of genes encoding proteins associated with the cell wall and membrane biosynthesis on the evolution of <italic>L. monocytogenes</italic>.</p> <p>Overall, genes in lineage I isolates were less likely to be affected by either recombination or positive selection, possibly reflecting that this lineage has experienced a recent bottleneck, as previously proposed [##REF##15066813##87##]. Frequent recombination in combination with positive selection of some genes in lineage II strains, on the other hand, may be important for the evolution of this generalist lineage, which is present in many different environments and host, including human clinical cases (although less common than lineage I). In combination with previous studies that have shown considerable differences in frequency of recombination and positive selection among different <italic>Streptococcus </italic>lineages and species [##REF##17475002##2##], our findings further show that even closely related bacterial lineages may differ in mechanisms contributing to their evolution.</p>
<title>Background</title> <p>The genus <italic>Listeria </italic>includes two closely related pathogenic and non-pathogenic species, <italic>L. monocytogenes </italic>and <italic>L. innocua</italic>. <italic>L. monocytogenes </italic>is an opportunistic human foodborne and animal pathogen that includes two common lineages. While lineage I is more commonly found among human listeriosis cases, lineage II appears to be overrepresented among isolates from foods and environmental sources. This study used the genome sequences for one <italic>L. innocua </italic>strain and four <italic>L. monocytogenes </italic>strains representing lineages I and II, to characterize the contributions of positive selection and recombination to the evolution of the <italic>L. innocua</italic>/<italic>L. monocytogenes </italic>core genome.</p> <title>Results</title> <p>Among the 2267 genes in the <italic>L. monocytogenes/L. innocua </italic>core genome, 1097 genes showed evidence for recombination and 36 genes showed evidence for positive selection. Positive selection was strongly associated with recombination. Specifically, 29 of the 36 genes under positive selection also showed evidence for recombination. Recombination was more common among isolates in lineage II than lineage I; this trend was confirmed by sequencing five genes in a larger isolate set. Positive selection was more abundant in the ancestral branch of lineage II (20 genes) as compared to the ancestral branch of lineage I (9 genes). Additional genes under positive selection were identified in the branch separating the two species; for this branch, genes in the role category "Cell wall and membrane biogenesis" were significantly more likely to have evidence for positive selection. Positive selection of three genes was confirmed in a larger isolate set, which also revealed occurrence of multiple premature stop codons in one positively selected gene involved in flagellar motility (<italic>flaR</italic>).</p> <title>Conclusion</title> <p>While recombination and positive selection both contribute to evolution of <italic>L. monocytogenes</italic>, the relative contributions of these evolutionary forces seem to differ by <italic>L. monocytogenes </italic>lineages and appear to be more important in the evolution of lineage II, which seems to be found in a broader range of environments, as compared to the apparently more host adapted lineage I. Diversification of cell wall and membrane biogenesis and motility-related genes may play a particularly important role in the evolution of <italic>L. monocytogenes</italic>.</p>
<title>Abbreviations</title> <p>COG: Clusters of Orthologous Groups of proteins; LRT: Likelihood Ratio Test; TO: Test Overall; an overall test for positive selection was carried out using the null model M1a (Nearly-neutral) and the alternative model M2a in PAML; TLM1A: <underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>I</underline><underline> A</underline>ncestral; this describes the branch-site test2 used to test for evidence of positive selection in the ancestral branch of <italic>L. monocytogenes </italic>lineage I.; TLM2A: <underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>II </underline><underline>A</underline>ncestral; this describes the branch-site test2 used to test for evidence of positive selection in the ancestral branch of <italic>L. monocytogenes </italic>lineage II; TLI/LM: <underline>T</underline>est <italic><underline>L</underline>. <underline>i</underline>nnocua</italic>/<italic><underline>L</underline>. <underline>m</underline>onocytogenes</italic>; this describes the branch-site test2 used to test for evidence of positive selection in the branch separating <italic>L. monocytogenes </italic>and <italic>L. innocua</italic>.</p> <title>Authors' contributions</title> <p>RHO outlined, performed, and interpreted the phylogenetic and statistical analyses, and drafted the manuscript. QS performed orthologous gene clustering and alignment, and implemented the analysis on the parallel computer cluster. MW supervised the project, participated in the design of the study and data interpretation, and finalized the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was partially supported by USDA Special Research Grant 2005-34459-15625 (to MW); the computer cluster used in the data analysis is partially funded by Microsoft. We thank Alphina Ho and Jasmine Badamo for help with PCR amplification and sequencing of some of the genes.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Unrooted phylogenetic tree of the five strains used in the genome-wide analyses.</bold> This tree represents the consensus tree of 2267 gene cluster alignments used for analyses. Branches tested for positive selection are indicated as TLM1A (<underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>I</underline><underline>A</underline>ncestral), TLM2A (<underline>T</underline>est <italic><underline>L</underline>. <underline>m</underline>onocytogenes </italic>lineage <underline>II</underline><underline>A</underline>ncestral), and TLI/LM (<underline>T</underline>est <italic><underline>L</underline>. <underline>i</underline>nnocua</italic>/<italic><underline>L</underline>. <underline>m</underline>onocytogenes</italic>).</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Schematic representations of the genes analyzed including (A) relative position of the orthologs in EGD-e, CLIP 11262 and F2365; and (B) circular chromosome of EGD-e.</bold> In "B" all protein coding, tRNA, and rRNA genes are shown in blue, brown, and purple; all genes analyzed are shown in green and genes with evidence for recombination (at least one test significant) are shown in red. There was no evidence for spatial clustering of genes with evidence for recombination (<italic>P </italic>= 0.957; U-test).</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Cumulative distribution of the <italic>p</italic>-values obtained from TLI/LM for selected COGs.</bold> "Overall" represents all <italic>p</italic>-values regardless of the COG classification. Genes involved in "Cell-wall/membrane biosynthesis", "Coenzyme transport and metabolism", and "Amino acid transport and metabolism" showed a tendency to have lower <italic>p</italic>-values in comparison with all genes analyzed, while genes involved in "Transcription" and "Inorganic ion transport and metabolism" showed a tendency to have higher <italic>p</italic>-values than all genes analyzed.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Recombination events identified by Clonal Frame, using the concatenated alignment of 40 randomly selected genes, in the external branches of the <italic>L. monocytogenes </italic>strains (A) H7858, (B) F2365, (C) F6854, and (D) EGD-e. </bold>Each of the 40 genes is represented between gray vertical lines. The order of the genes (left to right) is as follow: <italic>clpX </italic>(<italic>lmo1268</italic>), <italic>lmo0343</italic>, <italic>lmo0405</italic>, <italic>pflC </italic>(<italic>lmo1407</italic>), <italic>phoP </italic>(<italic>lmo2501</italic>), <italic>lmo1436</italic>, <italic>lmo1460</italic>, <italic>lmo1537</italic>, <italic>hemC </italic>(<italic>lmo1556</italic>), <italic>ccpA </italic>(<italic>lmo1599</italic>), <italic>lmo1623</italic>, <italic>lmo1790</italic>, <italic>lmo2262</italic>, <italic>pepC </italic>(<italic>lmo2338</italic>), <italic>lmo2391</italic>, <italic>trxB </italic>(<italic>lmo2478</italic>), <italic>lmo0190</italic>, <italic>lmo0860</italic>, <italic>lmo0877</italic>, <italic>lmo1087</italic>, <italic>proA </italic>(<italic>lmo1259</italic>), <italic>lmo0992</italic>, <italic>smbA </italic>(<italic>lmo1313</italic>), <italic>lmo1401</italic>, <italic>lmo1420</italic>, <italic>opuCC </italic>(<italic>lmo1426</italic>), <italic>trpD </italic>(<italic>lmo1631</italic>), <italic>lmo1693</italic>, <italic>purK </italic>(<italic>lmo1774</italic>), <italic>lmo1825</italic>, <italic>panB </italic>(<italic>lmo1902</italic>), <italic>lmo0028</italic>, <italic>lmo2175</italic>, <italic>lmo2348</italic>, <italic>lmo2566</italic>, <italic>lmo0487</italic>, <italic>lmo0878</italic>, <italic>lmo1004</italic>, <italic>lmo1011</italic>, <italic>cbiH </italic>(<italic>lmo1199</italic>). "x" indicate substitutions inferred to have occurred in the respective branches. Red lines represent the probability for each nucleotide to have been imported by means of recombination. Values at the bottom represent the position in the alignment in kilobases.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Phylogenetic consensus tree generated by ClonalFrame from the concatenated alignment of <italic>cheA</italic>, <italic>flaR</italic>, <italic>lmo0693</italic>, <italic>phoP </italic>and <italic>lmo2537 </italic>for 45 isolates.</bold> The 95% consensus phylogeny was obtained from two independent runs of ClonalFrame. This phylogeny clearly shows that the <italic>L. monocytogenes </italic>isolates form four distinct clusters, with lineage IIIA and IIIC (IIIA/C) isolates forming a sister group to lineage I isolates, while lineage IIIB isolates form an independent cluster that diverged earlier from the other isolates. Internal branches that showed evidence for recombination are labeled from A to J. Isolates with premature stop codons in <italic>flaR </italic>are marked with*.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Strains and genomes analyzed.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Strain</th><th align="center">Serotype</th><th align="center">Lineage</th><th align="center">Species</th><th align="center">Genome size (nt)</th><th align="center">No. of CDS<sup>(1)</sup></th><th align="center">Ref.</th></tr></thead><tbody><tr><td align="left">EGD-e</td><td align="center">1/2a</td><td align="center">II</td><td align="center"><italic>L. monocytogenes</italic></td><td align="center">2,944,528</td><td align="center">2846</td><td align="center">[##REF##11679669##49##]</td></tr><tr><td align="left">F6854</td><td align="center">1/2a</td><td align="center">II</td><td align="center"><italic>L. monocytogenes</italic></td><td align="center">2,953,211</td><td align="center">2945</td><td align="center">[##REF##15115801##51##]</td></tr><tr><td align="left">F2365</td><td align="center">4b</td><td align="center">I</td><td align="center"><italic>L. monocytogenes</italic></td><td align="center">2,905,310</td><td align="center">2821</td><td align="center">[##REF##15115801##51##]</td></tr><tr><td align="left">H7858</td><td align="center">4b</td><td align="center">I</td><td align="center"><italic>L. monocytogenes</italic></td><td align="center">2,893,921</td><td align="center">3007</td><td align="center">[##REF##15115801##51##]</td></tr><tr><td align="left">CLIP 11262</td><td align="center">6a</td><td align="center">-</td><td align="center"><italic>L. innocua</italic></td><td align="center">3,011,209</td><td align="center">2968</td><td align="center">[##REF##11679669##49##]</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Associations between COGs and descriptive variables.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">COG</th><th align="center">Number of gene analyzed</th><th align="center" colspan="5">Bonferroni-corrected P value for one-sided U-test testing for associations between genes in a given COG and<sup>(1)</sup></th></tr><tr><th/><th/><th colspan="5"><hr/></th></tr><tr><th/><th/><th align="center">&gt; length</th><th align="center">&gt; nt diversity</th><th align="center">&gt; Number of Informative sites</th><th align="center">&gt; Codon bias<sup>(2)</sup></th><th align="center">&lt; Codon bias<sup>(2)</sup></th></tr></thead><tbody><tr><td align="left">Energy production and conversion</td><td align="center">100</td><td align="center">&lt; 0.002</td><td align="center">ND</td><td align="center">NS</td><td align="center">0.007</td><td align="center">ND</td></tr><tr><td align="left">Cell cycle control, mitosis and meiosis</td><td align="center">20</td><td align="center">NS</td><td align="center">ND</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Amino acid transport and metabolism</td><td align="center">181</td><td align="center">&lt; 0.002</td><td align="center">NS</td><td align="center">&lt; 0.002</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Nucleotide transport and metabolism</td><td align="center">63</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Carbohydrate transport and metabolism</td><td align="center">186</td><td align="center">&lt; 0.002</td><td align="center">ND</td><td align="center">NS</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Coenzyme transport and metabolism</td><td align="center">87</td><td align="center">NS</td><td align="center">&lt; 0.001</td><td align="center">0.004</td><td align="center">ND</td><td align="center">0.047</td></tr><tr><td align="left">Lipid transport and metabolism</td><td align="center">51</td><td align="center">NS</td><td align="center">NS</td><td align="center">NS</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Translation</td><td align="center">91</td><td align="center">NS</td><td align="center">ND</td><td align="center">ND</td><td align="center">&lt; 0.001</td><td align="center">ND</td></tr><tr><td align="left">Transcription</td><td align="center">187</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">&lt; 0.002</td></tr><tr><td align="left">Replication, recombination and repair</td><td align="center">92</td><td align="center">&lt; 0.002</td><td align="center">&lt; 0.001</td><td align="center">&lt; 0.002</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Cell wall/membrane biogenesis</td><td align="center">76</td><td align="center">&lt; 0.002</td><td align="center">&lt; 0.001</td><td align="center">&lt; 0.002</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Cell motility</td><td align="center">40</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Posttranslational modification, protein turnover, chaperones</td><td align="center">52</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">0.021</td><td align="center">ND</td></tr><tr><td align="left">Inorganic ion transport and metabolism</td><td align="center">104</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Secondary metabolites biosynthesis, transport and catabolism</td><td align="center">31</td><td align="center">ND</td><td align="center">NS</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">General function prediction only</td><td align="center">254</td><td align="center">NS</td><td align="center">NS</td><td align="center">0.008</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Function unknown</td><td align="center">161</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Signal transduction mechanisms</td><td align="center">104</td><td align="center">NS</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td><td align="center">&lt; 0.002</td></tr><tr><td align="left">Intracellular trafficking and secretion</td><td align="center">37</td><td align="center">NS</td><td align="center">ND</td><td align="center">ND</td><td align="center">NS</td><td align="center">ND</td></tr><tr><td align="left">Defense mechanisms</td><td align="center">51</td><td align="center">&lt; 0.002</td><td align="center">NS</td><td align="center">0.020</td><td align="center">ND</td><td align="center">NS</td></tr><tr><td align="left">Not in COGs</td><td align="center">511</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">ND</td><td align="center">NS</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Association between COGs and recombination</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">COG</th><th align="center"><italic>P</italic>-values for association with recombination<sup>(1)</sup></th><th align="center" colspan="4"><italic>P</italic>-values for association with recombination test<sup>(2)</sup></th></tr><tr><th/><th/><th colspan="4"><hr/></th></tr><tr><th/><th/><th align="center">GENECONV</th><th align="center">NSS</th><th align="center"><bold>Maximum </bold>χ<sup>2</sup></th><th align="center">PHI</th></tr></thead><tbody><tr><td align="left">Carbohydrate transport and metabolism</td><td align="center">0.012</td><td align="center">0.001</td><td align="center">&lt; 0.001</td><td align="center">0.014</td><td align="center">0.032</td></tr><tr><td align="left">Amino acid transport and metabolism</td><td align="center">NS</td><td align="center">NS</td><td align="center">0.001</td><td align="center">0.026</td><td align="center">0.002</td></tr><tr><td align="left">Defense mechanisms</td><td align="center">NS</td><td align="center">NS</td><td align="center">0.022</td><td align="center">NS</td><td align="center">NS</td></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>Genes under positive selection.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Gene locus</th><th align="center">Gene description (gene symbol)</th><th align="center">COG<sup>(1)</sup></th><th align="center">Recombination<sup>(2)</sup></th><th align="center">Branch under pos. selection</th><th align="center">Q-value</th><th align="center">ω<sup>(3)</sup></th><th align="center">p<sup>(4)</sup></th><th align="center">BEB (<italic>P </italic>&gt; 95%)<sup>(5)</sup></th></tr></thead><tbody><tr><td align="left">Lmo0098</td><td align="center">PTS system, mannose/fructose/sorbose family, IID component (<italic>mptD</italic>)</td><td align="center">NCOG</td><td align="center">GCV; MAX</td><td align="center">LI/LM</td><td align="center">0.170</td><td align="center">472.58</td><td align="center">0.004</td><td align="center">-</td></tr><tr><td align="left">Lmo0099</td><td align="center">conserved hypothetical protein</td><td align="center">NCOG</td><td align="center">-</td><td align="center">LI/LM</td><td align="center">0.170</td><td align="center">∞</td><td align="center">0.009</td><td align="center">-</td></tr><tr><td align="left">Lmo0139</td><td align="center">conserved hypothetical protein</td><td align="center">NCOG</td><td align="center">-</td><td align="center">LM2A</td><td align="center">0.160</td><td align="center">56.74</td><td align="center">0.054</td><td align="center">95</td></tr><tr><td align="left">Lmo0297</td><td align="center">PRD/PTS system IIA 2 domain protein</td><td align="center">K; G; T</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM2A</td><td align="center">0.164</td><td align="center">∞</td><td align="center">0.012</td><td align="center">499</td></tr><tr><td align="left">Lmo0397</td><td align="center">conserved hypothetical protein</td><td align="center">S</td><td align="center">GCV; MAX</td><td align="center">LM2A</td><td align="center">0.1264</td><td align="center">∞</td><td align="center">0.043</td><td align="center">-</td></tr><tr><td align="left">Lmo0429</td><td align="center">glycosyl hydrolase, family 38</td><td align="center">G</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM2A</td><td align="center">0.0325</td><td align="center">977.51</td><td align="center">0.008</td><td align="center">667</td></tr><tr><td align="left">Lmo0455</td><td align="center">conserved hypothetical protein</td><td align="center">T; Q</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM2A</td><td align="center">0.1214</td><td align="center">∞</td><td align="center">0.004</td><td align="center">-</td></tr><tr><td align="left">Lmo0653</td><td align="center">conserved hypothetical protein</td><td align="center">S</td><td align="center">MAX</td><td align="center">LM2A</td><td align="center">0.033</td><td align="center">∞</td><td align="center">0.012</td><td align="center">306</td></tr><tr><td align="left">Lmo0658</td><td align="center">endonuclease III domain protein</td><td align="center">L</td><td align="center">MAX; NSS</td><td align="center">LM2A</td><td align="center">0.108</td><td align="center">∞</td><td align="center">0.023</td><td align="center">209</td></tr><tr><td align="left">Lmo0692</td><td align="center">chemotaxis protein CheA (<italic>cheA</italic>)</td><td align="center">T; N</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM1A</td><td align="center">0.156</td><td align="center">∞</td><td align="center">0.002</td><td align="center">-</td></tr><tr><td align="left">Lmo0693</td><td align="center">flagellar motor switch domain protein</td><td align="center">N; U</td><td align="center">-</td><td align="center">LM2A</td><td align="center">0.007</td><td align="center">∞</td><td align="center">0.022</td><td align="center">-</td></tr><tr><td align="left">Lmo0695</td><td align="center">conserved hypothetical protein</td><td align="center">NCOG</td><td align="center">MAX; NSS; PHI</td><td align="center">LM1A</td><td align="center">0.175</td><td align="center">∞</td><td align="center">0.014</td><td align="center">-</td></tr><tr><td align="left">Lmo0732</td><td align="center">cell wall surface anchor family protein</td><td align="center">NCOG</td><td align="center">GCV</td><td align="center">LM1A</td><td align="center">0.185</td><td align="center">6.92</td><td align="center">0.076</td><td align="center">-</td></tr><tr><td align="left">Lmo0782</td><td align="center">PTS system, mannose/fructose/sorbose family, IIC component (<italic>mpoD</italic>)</td><td align="center">NCOG</td><td align="center">GCV; MAX</td><td align="center">Overall;</td><td align="center">0.146;</td><td align="center">88.23;</td><td align="center">0.008;</td><td align="center">-</td></tr><tr><td/><td/><td/><td/><td align="center">LM1A;</td><td align="center">0.052;</td><td align="center">∞;</td><td align="center">0.0001;</td><td align="center">-</td></tr><tr><td/><td/><td/><td/><td align="center">LM2A</td><td align="center">0.098</td><td align="center">∞</td><td align="center">0.004</td><td align="center">-</td></tr><tr><td align="left">Lmo0785</td><td align="center">sigma-54 dependent Kal regulator (<italic>manR</italic>)</td><td align="center">K; T</td><td align="center">GCV; MAX</td><td align="center">LM2A</td><td align="center">0.129</td><td align="center">1.00</td><td align="center">0.000</td><td align="center">-</td></tr><tr><td align="left">Lmo0872</td><td align="center">major facilitator family transporter</td><td align="center">G</td><td align="center">GCV; MAX; NSS</td><td align="center">LM2A</td><td align="center">0.137</td><td align="center">∞</td><td align="center">0.008</td><td align="center">-</td></tr><tr><td align="left">Lmo0910</td><td align="center">putative membrane protein</td><td align="center">R</td><td align="center">MAX; NSS</td><td align="center">LM1A</td><td align="center">0.019</td><td align="center">∞</td><td align="center">0.012</td><td align="center">-</td></tr><tr><td align="left">Lmo1146</td><td align="center">conserved hypothetical protein</td><td align="center">NCOG</td><td align="center">GCV; MAX</td><td align="center">LI/LM</td><td align="center">0.046</td><td align="center">223.37</td><td align="center">0.023</td><td align="center">169</td></tr><tr><td align="left">Lmo1164</td><td align="center">PduO protein (<italic>pduO</italic>)</td><td align="center">S; R</td><td align="center">GCV</td><td align="center">LI/LM</td><td align="center">0.170</td><td align="center">195.87</td><td align="center">0.038</td><td align="center">-</td></tr><tr><td align="left">Lmo1412</td><td align="center">DNA topology modulation protein FlaR (<italic>flaR</italic>)</td><td align="center">F</td><td align="center">-</td><td align="center">LM2A</td><td align="center">0.160</td><td align="center">8.56</td><td align="center">0.139</td><td align="center">12; 37; 68</td></tr><tr><td align="left">Lmo1424</td><td align="center">transporter, NRAMP family (<italic>mntH</italic>)</td><td align="center">P</td><td align="center">GCV</td><td align="center">LM1A</td><td align="center">0.185</td><td align="center">1.00</td><td align="center">0.000</td><td align="center">-</td></tr><tr><td align="left">Lmo1523</td><td align="center">GTP pyrophosphokinase (<italic>relA</italic>)</td><td align="center">K; T</td><td align="center">-</td><td align="center">LM2A</td><td align="center">0.033</td><td align="center">1.00</td><td align="center">0.000</td><td align="center">-</td></tr><tr><td align="left">Lmo1529</td><td align="center">preprotein translocase, YajC subunit</td><td align="center">U</td><td align="center">-</td><td align="center">LI/LM</td><td align="center">0.170</td><td align="center">∞</td><td align="center">0.011</td><td align="center">-</td></tr><tr><td align="left">Lmo2102</td><td align="center">glutamine amidotransferase, SNO family (<italic>pdxT</italic>)</td><td align="center">H</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM1A</td><td align="center">0.185</td><td align="center">∞</td><td align="center">0.017</td><td align="center">66</td></tr><tr><td align="left">Lmo2121</td><td align="center">glycosyl transferase, family 65</td><td align="center">G</td><td align="center">GCV; MAX; NSS</td><td align="center">LM1A</td><td align="center">7.6E-06</td><td align="center">35.31</td><td align="center">0.031</td><td align="center">722; 723; 725; 729; 730; 744; 752;</td></tr><tr><td align="left">Lmo2178</td><td align="center">cell wall surface anchor family protein</td><td align="center">M</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LI/LM;</td><td align="center">0.193;</td><td align="center">15.92;</td><td align="center">0.001;</td><td align="center">-</td></tr><tr><td/><td/><td/><td/><td align="center">LM2A</td><td align="center">0.137</td><td align="center">5.34</td><td align="center">0.025</td><td align="center">1769</td></tr><tr><td align="left">Lmo2215</td><td align="center">ABC transporter, ATP-binding protein</td><td align="center">V</td><td align="center">MAX; NSS</td><td align="center">Overall</td><td align="center">0.188</td><td align="center">14.23</td><td align="center">0.008</td><td align="center">-</td></tr><tr><td align="left">Lmo2222</td><td align="center">Ser/Thr protein phosphatase family protein</td><td align="center">L</td><td align="center">GCV; MAX</td><td align="center">LM2A</td><td align="center">0.160</td><td align="center">∞</td><td align="center">0.021</td><td align="center">253</td></tr><tr><td align="left">Lmo2446</td><td align="center">glycosyl hydrolase, family 31</td><td align="center">G</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM2A</td><td align="center">0.003</td><td align="center">∞</td><td align="center">0.0001</td><td align="center">-</td></tr><tr><td align="left">Lmo2596</td><td align="center">ribosomal protein S9 (<italic>rpsI</italic>)</td><td align="center">NCOG</td><td align="center">-</td><td align="center">LM2A</td><td align="center">0.021</td><td align="center">∞</td><td align="center">0.010</td><td align="center">-</td></tr><tr><td align="left">Lmo2611</td><td align="center">adenylate kinase (<italic>adk</italic>)</td><td align="center">F</td><td align="center">GCV; MAX</td><td align="center">LM2A</td><td align="center">0.120</td><td align="center">24.80</td><td align="center">0.013</td><td align="center">-</td></tr><tr><td align="left">Lmo2724</td><td align="center">conserved hypothetical protein</td><td align="center">S</td><td align="center">NSS</td><td align="center">LI/LM</td><td align="center">0.170</td><td align="center">∞</td><td align="center">0.008</td><td align="center">-</td></tr><tr><td align="left">Lmo2802</td><td align="center">glucose-inhibited division protein B (<italic>girB</italic>)</td><td align="center">M</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">LM1A</td><td align="center">0.046</td><td align="center">∞</td><td align="center">0.013</td><td align="center">-</td></tr><tr><td align="left">Lmo2804</td><td align="center">conserved hypothetical protein</td><td align="center">NCOG</td><td align="center">GCV; MAX; NSS; PHI</td><td align="center">Overall</td><td align="center">6.5E-17</td><td align="center">16.53</td><td align="center">0.044</td><td align="center">-</td></tr><tr><td align="left">Lmo2824</td><td align="center">D-isomer specific 2- hydroxyacid dehydrogenase family protein</td><td align="center">E; H</td><td align="center">GCV; MAX; NSS</td><td align="center">LM2A</td><td align="center">0.160</td><td align="center">229.62</td><td align="center">0.003</td><td align="center">-</td></tr></tbody></table></table-wrap> <table-wrap id="T5" position="float"><label>Table 5</label><caption><p>Positive selection and recombination analyses of 5 genes in 45 isolates<sup>(1)</sup></p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Gene</th><th align="center">Function</th><th align="center">Recombination evidence<sup>(2) </sup>(<italic>p</italic>-value)</th><th align="center">Positive selection evidence<sup>(3) </sup>(<italic>p</italic>-value)</th><th align="center">ω<sup>(4)</sup></th><th align="center">p<sup>(5)</sup></th><th align="center">BEB<sup>(6) </sup>sites (probability)</th></tr></thead><tbody><tr><td align="left"><italic>cheA</italic></td><td align="center">Two-component sensor histidine kinase CheA, involved in chemotaxis (Dons <italic>et a</italic>l., 2004)</td><td align="center">GENECONV (&lt; 0.001), NSS (&lt; 0.001), Max χ<sup>2 </sup>(&lt; 0.001), PHI (&lt; 0.001)</td><td align="center">LIIIA/C-LI (&lt; 0.001)</td><td align="center">∞</td><td align="center">0.002</td><td align="center">140 (98%);</td></tr><tr><td align="left"><italic>lmo0693</italic></td><td align="center">Putative flagellar motor switch protein, involved in motility</td><td align="center">NSS (0.033)</td><td align="center">LII (&lt; 0.001)</td><td align="center">∞</td><td align="center">0.022</td><td align="center">17 (73%)<break/>18 (98%)</td></tr><tr><td align="left"><italic>flaR</italic></td><td align="center">Histone like-DNA topology modulator, involved in regulation of flagellin expression (Sanchez-Campillo <italic>et al</italic>., 1995)</td><td align="center">GENECONV (0.010), Max χ<sup>2</sup>(0.010), PHI (0.005)</td><td align="center">LII (0.002)</td><td align="center">14.1</td><td align="center">0.046</td><td align="center">4 (90%)<break/>12 (99%)<break/>68 (80%)</td></tr><tr><td align="left"><italic>phoP</italic></td><td align="center">Putative two-component response phosphate regulator PhoP</td><td align="center">GENECONV (0.014), NSS (0.001), Max χ<sup>2 </sup>(0.020), PHI (0.003)</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="left"><italic>lmo2537</italic></td><td align="center">Putative UDP-N- acetylglucosamine-2-epimerase, involved in teichoic acid biogenesis (Dubail <italic>et al</italic>., 2006)</td><td align="center">GENECONV (&lt; 0.001), NSS (&lt; 0.001), Max χ<sup>2 </sup>(&lt; 0.001), PHI (&lt; 0.001)</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Isolates used to confirm positive selection and recombination patterns in five selected genes.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Primers used for re-sequencing.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Alignment of the <italic>flaR </italic>sequence for seven isolates with premature stop codons.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>(1) </sup>number of coding sequences used for cluster analysis.</p></table-wrap-foot> <table-wrap-foot><p><sup>(1) </sup>"&gt;" or "&lt;" indicates the direction of the one-sided tests (i.e. column "&gt; Codon bias" shows Bonferroni-corrected <italic>p</italic>-values for associations between genes in a given COG and higher codon bias (as compared to the genes in other COGs), while the column "&lt; Codon bias" test shows Bonferroni-corrected <italic>p</italic>-values for genes in COGs with lower codon bias as compared to genes in other COGs; "ND" = Not determined (tests were not performed for COGs that showed values that were not consistent with the tested alternative hypothesis, e.g., if the average gene length for genes in a given COG was below average than we did not test for an association of this COG with increased gene length); "NS". Not significant.</p><p><sup>(2)</sup>Tests for codon bias were performed using Nc values; a lower Nc indicates increased codon bias.</p></table-wrap-foot> <table-wrap-foot><p><sup>(1) </sup>Fisher's exact test with Bonferroni correction; NS = not significant;</p><p><sup>(2) </sup>U-test with Bonferroni correction.</p></table-wrap-foot> <table-wrap-foot><p><sup>(1) </sup>NCOG: Not in COGs; K: Transcription; G: Carbohydrate transport and metabolism; T: Signal transduction mechanisms; S: Function unknown; Q: Secondary metabolites biosynthesis, transport and catabolism; L: Replication, recombination and repair; N: Cell motility; U: Intracellular trafficking and secretion; R: General function prediction only; F: Nucleotide transport and metabolism; P: Inorganic ion transport and metabolism; E: Amino acid transport and metabolism; H: Coenzyme transport and metabolism; M: Cell wall/membrane biogenesis; V: Defense mechanisms</p><p><sup>(2) </sup>GCV: GENECONV; MAX: Maximum χ<sup>2</sup>; NSS: Neighbour Similarity Score; PHI: Pairwise Homoplasy Test</p><p><sup>(3) </sup>ω = d<sub>N</sub>/d<sub>S </sub>(Number of nonsynonymous changes per nonsynonymous sites/Number of synonymous changes per synonymous sites); infinite values of ω (∞) indicate that the model did not find synonymous changes for the branches tested (d<sub>S </sub>= 0; ω ~ ∞). However, this (i.e., ω ~ ∞) does not affect the validity of the Likelihood Ratio Test, which was used to identify the genes under positive selection (Z. Yang, pers. Communication; see <ext-link ext-link-type="uri" xlink:href="http://gsf.gc.ucdavis.edu/viewtopic.php?f=1&amp;t=2079&amp;sid=c4a82e1ca334ca84a00a8b85e0f33c9d">http://gsf.gc.ucdavis.edu/viewtopic.php?f=1&amp;t=2079&amp;sid=c4a82e1ca334ca84a00a8b85e0f33c9d</ext-link> and <ext-link ext-link-type="uri" xlink:href="http://gsf.gc.ucdavis.edu/viewtopic.php?f=1&amp;t=2329&amp;sid=c4a82e1ca334ca84a00a8b85e0f33c9d">http://gsf.gc.ucdavis.edu/viewtopic.php?f=1&amp;t=2329&amp;sid=c4a82e1ca334ca84a00a8b85e0f33c9d</ext-link>;</p><p><sup>(4) </sup>Proportion of sites under positive selection</p><p><sup>(5) </sup>This column lists sites identified using Bayes Empirical Bayes (BEB) as being under positive selection; numbers identify the amino acid sites (in alignments) that are under positive selection</p></table-wrap-foot> <table-wrap-foot><p><sup>(1) </sup>Analyses were performed using an alignment of these 5 genes for the 40 <italic>L. monocytogenes i</italic>solates, for which these genes were sequenced here (see Supp. Table 1), as well as the four <italic>L. monocytogenes </italic>and one <italic>L. innocua </italic>strain for which full genome sequences were available (Table 1).</p><p><sup>(2) </sup>NSS: Neighbour Similarity Score; Max χ<sup>2</sup>: Maximum χ<sup>2</sup>; PHI: Pairwise Homoplasy Index; GENECONV performed with g-scale = 1 did not show significant inner fragments for any of the five genes, however four genes showed significant inner fragments in GENECONV performed with g-scale = 2, these <italic>p</italic>-values are reported here. The g-scale setting in GENECONV is associated with the number of polymorphisms allowed in a putative recombinant fragment; more polymorphisms are allowed as the g-scale value decreases from 3 to 1,</p><p><sup>(3) </sup>Branches where positive selection was identified. LIIIA/C-LI: ancestral branch of lineages IIIA/C and I isolates (branch B in Figure 5); LII: Ancestral branch of lineage II isolates (branch F in Figure 5); there was no evidence for positive selection in <italic>phoP </italic>and <italic>lmo2537</italic>.</p><p><sup>(4) </sup>ω = d<sub>N</sub>/d<sub>S </sub>(number of nonsynonymous changes per nonsynonymous sites/Number of synonymous changes per synonymous sites); infinite values of ω (∞) indicate that the model did not find synonymous changes for the branches tested (d<sub>S </sub>= 0; ω ~ ∞).</p><p><sup>(5) </sup>Proportion of sites under positive selection.</p><p><sup>(6) </sup>this column lists sites identified using Bayes Empirical Bayes (BEB) as being under positive selection; numbers identify the amino acid sites (in alignments) that are under positive selection; "probabilities" refer to the posterior probabilities that the respective sites evolved by positive selection.</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-233-1"/> <graphic xlink:href="1471-2148-8-233-2"/> <graphic xlink:href="1471-2148-8-233-3"/> <graphic xlink:href="1471-2148-8-233-4"/> <graphic xlink:href="1471-2148-8-233-5"/>
<media xlink:href="1471-2148-8-233-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-233-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-233-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Petran", "Zottola"], "given-names": ["RL", "EA"], "article-title": ["A study of factors affecting growth and recovery of "], "italic": ["Listeria monocytogenes "], "source": ["J Food Science"], "year": ["1989"], "volume": ["54"], "fpage": ["458"], "lpage": ["460"], "pub-id": ["10.1111/j.1365-2621.1989.tb03105.x"]}, {"article-title": ["J Craig Venter Institute"], "ext-link": ["http://cmr.jcvi.org/"]}, {"surname": ["Hall"], "given-names": ["TA"], "article-title": ["BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT"], "source": ["Nucl Acids Symp Ser"], "year": ["1999"], "volume": ["41"], "fpage": ["95"], "lpage": ["98"]}, {"surname": ["Sawyer"], "given-names": ["SA"], "source": ["GENECONV: A computer package for the statistical detection of gene conversion"], "year": ["1999"], "publisher-name": ["Distributed by the author, Department of mathematics, Washington University in St louis"], "ext-link": ["http://www.math.wustl.edu/%7Esawyer/geneconv/"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate: a practical and powerful approach to multiple testing"], "source": ["Journal of the Royal Statistical Society B"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}, {"surname": ["Severino", "Dussurget", "Vencio", "Dumas", "Garrido", "Padilla", "Piveteau", "Lemaitre", "Kunst", "Glaser", "Buchrieser"], "given-names": ["P", "O", "RZ", "E", "P", "G", "P", "JP", "F", "P", "C"], "article-title": ["Comparative transcriptome analysis of "], "italic": ["Listeria monocytogenes "], "source": ["Appl Environ Microbiol"], "year": ["2007"]}]
BMC Evol Biol. 2008 Aug 12; 8:233