Evolutionary genetics of the capsular locus of serogroup 6 pneumococci

The evolution of the capsular biosynthetic (cps) locus of serogroup 6 Streptococcus pneumoniae was investigated by analyzing sequence variation within three serotype-specific cps genes from 102 serotype 6A and 6B isolates. Sequence variation within these cps genes was related to the genetic relatedness of the isolates, determined by multilocus sequence typing, and to the inferred patterns of recent evolutionary descent, explored using the eBURST algorithm. The serotype-specific cps genes had a low percent G+C, and there was a low level of sequence diversity in this region among serotype 6A and 6B isolates. There was also little sequence divergence between these serotypes, suggesting a single introduction of an ancestral cps sequence, followed by slight divergence to create serotypes 6A and 6B. A minority of serotype 6B isolates had cps sequences (class 2 sequences) that were approximately 5% divergent from those of other serotype 6B isolates (class 1 sequences) and which may have arisen by a second, more recent introduction from a related but distinct source. Expression of a serotype 6A or 6B capsule correlated perfectly with a single nonsynonymous polymorphism within wciP, the rhamnosyl transferase gene. In addition to ample evidence of the horizontal transfer of the serotype 6A and 6B cps locus into unrelated lineages, there was evidence for relatively frequent changes from serotype 6A to 6B, and vice versa, among very closely related isolates and examples of recent recombinational events between class 1 and 2 cps serogroup 6 sequences.     

Sequence diversity and molecular evolution of the leukotoxin (lktA) gene in bovine and ovine strains of Mannheimia (Pasteurella) haemolytica

The molecular evolution of the leukotoxin structural gene (lktA) of Mannheimia (Pasteurella) haemolytica was investigated by nucleotide sequence comparison of lktA in 31 bovine and ovine strains representing the various evolutionary lineages and serotypes of the species. Eight major allelic variants (1.4 to 15.7% nucleotide divergence) were identified; these have mosaic structures of varying degrees of complexity reflecting a history of horizontal gene transfer and extensive intragenic recombination. The presence of identical alleles in strains of different genetic backgrounds suggests that assortative (entire gene) recombination has also contributed to strain diversification in M. haemolytica. Five allelic variants occur only in ovine strains and consist of recombinant segments derived from as many as four different sources. Four of these alleles consist of DNA (52.8 to 96.7%) derived from the lktA gene of the two related species Mannheimia glucosida and Pasteurella trehalosi, and four contain recombinant segments derived from an allele that is associated exclusively with bovine or bovine-like serotype A2 strains. The two major lineages of ovine serotype A2 strains possess lktA alleles that have very different evolutionary histories and encode divergent leukotoxins (5.3% amino acid divergence), but both contain segments derived from the bovine allele. Homologous segments of donor and recipient alleles are identical or nearly identical, indicating that the recombination events are relatively recent and probably postdate the domestication of cattle and sheep. Our findings suggest that host switching of bovine strains from cattle to sheep, together with inter- and intraspecies recombinational exchanges, has played an important role in generating leukotoxin diversity in ovine strains. In contrast, there is limited allelic diversity of lktA in bovine strains, suggesting that transmission of strains from sheep to cattle has been less important in leukotoxin evolution.     

Evidence of recombination among enteroviruses.

Human enteroviruses consist of more than 60 serotypes, reflecting a wide range of evolutionary divergence. They have been genetically classified into four clusters on the basis of sequence homology in the coding region of the single-stranded RNA genome. To explore further the genetic relationships between human enteroviruses and to characterize the evolutionary mechanisms responsible for variation, previously sequenced genomes were subjected to detailed comparison. Bootstrap and genetic similarity analyses were used to systematically scan the alignments of complete genomic sequences. Bootstrap analysis provided evidence from an early recombination event at the junction of the 5' noncoding and coding regions of the progenitors of the current clusters. Analysis within the genetic clusters indicated that enterovirus prototype strains include intraspecies recombinants. Recombination breakpoints were detected in all genomic regions except the capsid protein coding region. Our results suggest that recombination is a significant and relatively frequent mechanism in the evolution of enterovirus genomes.     

Recombination Detection Under Evolutionary Scenarios Relevant to Functional Divergence

Recombination can negatively impact methods designed to detect divergent gene function that rely on explicit knowledge of a gene tree. However, we know little about how recombination detection methods perform under evolutionary scenarios encountered in studies of functional molecular divergence. We use simulation to evaluate false positive rates for six recombination detection methods (GENECONV, MaxChi, Chimera, RDP, GARD-SBP, GARD-MBP) under evolutionary scenarios that might increase false positives. Broadly, these scenarios address: (i) asymmetric tree topology and sequence divergence, (ii) non-stationary codon bias and selection pressure, and (iii) positive selection. We also evaluate power to detect recombination under truly recombinant history. As with previous studies, we find that power increases with sequence divergence. However, we also find that accuracy to correctly infer the number of breakpoints is extremely low. When recombination is absent, increased sequence divergence leads to increased false positives. Furthermore, one method (GARD-SBP) is sensitive to tree shape, with higher false positive rates under an asymmetric tree topology. Somewhat surprisingly, all methods are robust to the simulated heterogeneity in codon bias, shifts in selection pressure and presence of positive selection. Based on these findings, we recommend that studies of functional divergence in systems where recombination is plausible can, and should, include a pre-test for recombination. Application of all methods to the core genome of Prochlorococcus reveals a substantial lack of concordance among results. Based on analysis of both real and simulated datasets we present some guidelines for the investigation of recombination in genes that may have experienced functional divergence.     

Distribution of the SsuDAT1I restriction-modification system among different serotypes of Streptococcus suis

The SsuDAT1I restriction-modification (R-M) system, which contains two methyltransferases and two restriction endonucleases with recognition sequence 5'-GATC-3', was first found in a field isolate of Streptococcus suis serotype 2. Isoschizomers of the R-M system were found in the same locus between purH and purD in a field isolate of serotype 1/2 and the reference strains of serotypes 3, 7, 23, and 26 among 29 strains of different serotypes examined in this study. The R-M gene sequences in serotypes 1/2, 3, 7, and 23 were very similar to those of SsuDAT1I, whereas those in serotype 26 were less similar. These results indicate intraspecies recombination among them and genetic divergence through their evolution.     

Inference of Functional Divergence Among Proteins When the Evolutionary Process is Non-stationary

Functional shifts during protein evolution are expected to yield shifts in substitution rate, and statistical methods can test for this at both codon and amino acid levels. Although methods based on models of sequence evolution serve as powerful tools for studying evolutionary processes, violating underlying assumptions can lead to false biological conclusions. It is not unusual for functional shifts to be accompanied by changes in other aspects of the evolutionary process, such as codon or amino acid frequencies. However, models used to test for functional divergence assume these frequencies remain constant over time. We employed simulation to investigate the impact of non-stationary evolution on functional divergence inference. We investigated three likelihood ratio tests based on codon models and found varying degrees of sensitivity. Joint effects of shifts in frequencies and selection pressures can be large, leading to false signals for positive selection. Amino acid-based tests (FunDi and Bivar) were also compromised when several aspects of the substitution process were not adequately modeled. We applied the same tests to a core genome “scan” for functional divergence between light-adapted ecotypes of the cyanobacteria Prochlorococcus, and carried out gene-specific simulations for ten genes. Results of those simulations illustrated how the inference of functional divergence at the genomic level can be seriously impacted by model misspecification. Although computationally costly, simulations motivated by data in hand are warranted when several aspects of the substitution process are either misspecified or not included in the models upon which the statistical tests were built.     

载脂蛋白基因家族的密码子空间分析——核苷酸替换的非随机进化选择

ＤＮＡ分子进化中,对核苷酸替换的选择可呈选择中性或选择倾向性。为研究载脂蛋白基因进化过程中对核苷酸变化的选择方式,本文建立了基因的密码子空间分析方法。密码子空间是由密码子３个位置上核苷酸出现机率所组成的矩阵。对该空间中核苷酸分布的非随机性度量可以反映进化过程中核苷酸替换的选择方式。应用该法,我们发现载脂蛋白基因密码子空间第一及第三位的核苷酸分布呈高度非随机性。进一步研究表明：这种核苷酸的非随机分布可能与腺苷酸、胸苷酸对密码子位置的非中性选择有关。此外,还研究了同义密码子的选择使用与分支种系发生的关系。结果显示：载脂蛋白分子演化中存在着同义密码子使用的分子进化钟。这些研究提示密码子空间中核苷酸替换的非随机选择可能是载脂蛋白基因进化的一种特征。     

RNA recombination plays a major role in genomic change during circulation of coxsackie B viruses

RNA recombination has been shown to occur during circulation of enteroviruses, but most studies have focused on poliovirus. To examine the role of recombination in the evolution of the coxsackie B viruses (CVB), we determined the partial sequences of four genomic intervals for multiple clinical isolates of each of the six CVB serotypes isolated from 1970 to 1996. The regions sequenced were the 5'-nontranslated region (5'-NTR) (350 nucleotides [nt]), capsid (VP4-VP2, 416 nt, and VP1, approximately 320 nt), and polymerase (3D, 491 nt). Phylogenetic trees were constructed for each genome region, using the clinical isolate sequences and those of the prototype strains of all 65 enterovirus serotypes. The partial VP1 sequences of each CVB serotype were monophyletic with respect to serotype, as were the VP4-VP2 sequences, in agreement with previously published studies. In some cases, however, incongruent tree topologies suggested that intraserotypic recombination had occurred between the sequenced portions of VP2 and VP1. Outside the capsid region, however, isolates of the same serotype were not monophyletic, indicating that recombination had occurred between the 5'-NTR and capsid, the capsid and 3D, or both. Almost all clinical isolates were recombinant relative to the prototype strain of the same serotype. All of the recombination partners appear to be members of human enterovirus species B. These results suggest that recombination is a frequent event during enterovirus evolution but that there are genetic restrictions that may influence recombinational compatibility.     

Frequency and dynamics of recombination within different species of human enteroviruses

Enteroviruses are members of the family Picornaviridae that cause widespread infections in human and other mammalian populations. Enteroviruses are genetically and antigenically highly variable, and recombination within and between serotypes contributes to their genetic diversity. To investigate the dynamics of the recombination process, sequence phylogenies between three regions of the genome (VP4, VP1, and 3Dpol) were compared among species A and B enterovirus variants detected in a human population-based survey in Scotland between 2000 and 2001, along with contemporary virus isolates collected in the same geographical region. This analysis used novel bioinformatic methods to quantify phylogenetic compatibility and correlations with serotype assignments of evolutionary trees constructed for different regions of the enterovirus genome. Species B enteroviruses showed much more frequent, time-correlated recombination events than those found for species A, despite the equivalence in population sampling, concordant with a linkage analysis of previously characterized enterovirus sequences obtained over longer collection periods. An analysis of recombination among complete genome sequences by computation of a phylogenetic compatibility matrix (PCM) demonstrated sharply defined boundaries between the VP2/VP3/VP1 block and sequences to either side in phylogenetic compatibility. The PCM also revealed equivalent or frequently greater degrees of incompatibility between different parts within the nonstructural region (2A-3D), indicating the occurrence of extensive recombination events in the past evolution of this part of the genome. Together, these findings provide new insights into the dynamics of species A and B enterovirus recombination and evolution and into the contribution of structured sampling to documenting reservoirs, emergence, and spread of novel recombinant forms in human populations.     

Porcine teschoviruses comprise at least eleven distinct serotypes: molecular and evolutionary aspects

Nucleotide sequencing and phylogenetic analysis of 10 recognized prototype strains of the porcine enterovirus (PEV) cytopathic effect (CPE) group I reveals a close relationship of the viral genomes to the previously sequenced strain F65, supporting the concept of a reclassification of this virus group into a new picornavirus genus. Also, nucleotide sequences of the polyprotein-encoding genome region or the P1 region of 28 historic strains and recent field isolates were determined. The data suggest that several closely related but antigenically and molecular distinct serotypes constitute one species within the proposed genus Teschovirus. Based on sequence data and serological data, we propose a new serotype with strain Dresden as prototype. This hitherto unrecognized serotype is closely related to porcine teschovirus 1 (PTV-1, former PEV-1), but induces type-specific neutralizing antibodies. Sequencing of field isolates collected from animals presenting with neurological disorders prove that other serotypes than PTV-1 may also cause polioencephalomyelitis of swine.     

Comparative sequence analysis of the hexon gene in the entire spectrum of human adenovirus serotypes: phylogenetic, taxonomic, and clinical implications

The adenovirus (AdV) hexon constitutes the major virus capsid protein. The epitopes located on the hexon protein are targets of neutralizing antibodies in vivo, serve in the recognition by cytotoxic T cells, and provide the basis for the classification of adenoviruses into the 51 serotypes known to date. We have sequenced the entire hexon gene from human serotypes with incomplete or no sequence information available (n = 34) and performed a comparative analysis of all sequences. The overall sequence divergence between the 51 human serotypes ranged from 0.7 to 25.4% at the protein level. The sequence information has been exploited to assess the phylogeny of the adenovirus family, and protein distances between the six AdV species (A to F) and among individual serotypes within each species were calculated. The analysis revealed that the differences among serotypes within individual species range from 0.3 to 5.4% in the conserved regions (765 amino acids [aa]) and from 1.5 to 59.6% in the variable regions (154 to 221 aa). Serotypes of different species showed an expectedly greater divergence both in the conserved (5.9 to 12.3%) and variable (49.0 to 74.7%) regions. Construction of a phylogenetic tree revealed three major clades comprising the species B+D+E, A+F, and C, respectively. For serotypes 50 and 51, the original assignment to species B and D, respectively, is not in accordance with the hexon DNA and protein sequence data, which placed serotype 50 within species D and serotype 51 within species B. Moreover, the hexon gene of serotype 16, a member of species B, was identified as the product of recombination between sequences of species B and E. In addition to providing a basis for improved molecular diagnostics and classification, the elucidation of the complete hexon gene sequence in all AdV serotypes yields information on putative epitopes for virus recognition, which may have important implications for future treatment strategies permitting efficient targeting of any AdV serotype.     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.     

Comparative evolution of GII.3 and GII.4 norovirus over a 31-year period

Noroviruses are the most common cause of epidemic gastroenteritis. Genotype II.3 is one of the most frequently detected noroviruses associated with sporadic infections. We studied the evolution of the major capsid gene from seven archival GII.3 noroviruses collected during a cross-sectional study at the Children's Hospital in Washington, DC, from 1975 through 1991, together with capsid sequence from 56 strains available in GenBank. Evolutionary analysis concluded that GII.3 viruses evolved at a rate of 4.16 × 10(-3) nucleotide substitutions/site/year (strict clock), which is similar to that described for the more prevalent GII.4 noroviruses. The analysis of the amino acid changes over the 31-year period found that GII.3 viruses evolve at a relatively steady state, maintaining 4% distance, and have a tendency to revert back to previously used residues while preserving the same carbohydrate binding profile. In contrast, GII.4 viruses demonstrate increasing rates of distance over time because of the continued integration of new amino acids and changing HBGA binding patterns. In GII.3 strains, seven sites acting under positive selection were predicted to be surface-exposed residues in the P2 domain, in contrast to GII.4 positively selected sites located primarily in the shell domain. Our study suggests that GII.3 noroviruses caused disease as early as 1975 and that they evolve via a specific pattern, responding to selective pressures induced by the host rather than presenting a nucleotide evolution rate lower than that of GII.4 noroviruses, as previously proposed. Understanding the evolutionary dynamics of prevalent noroviruses is relevant to the development of effective prevention and control strategies.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

Environmental surveillance and molecular characterization of human enteric viruses in tropical urban wastewaters

Aims: To study the prevalence and genotypes of waterborne pathogenic viruses in urban wastewaters in the tropical region. Methods and Results: Viruses in wastewaters collected at three water reclamation plants in Singapore were studied by molecular methods. Over a 6‐month sampling period, adenoviruses, astroviruses and both norovirus genogroups I (GI) and II (GII) were detected in 100% of the sewage and secondary effluent. Enteroviruses and hepatitis A viruses (HAV) were found in 94 and 78% of sewage, and 89 and 28% of secondary effluent, respectively. By using quantitative real‐time PCR, estimated concentrations of astrovirus in the sewage were 1–2 orders of magnitude higher than those for adenovirus, noroviruses GI and GII. Genotyping of environmental isolates revealed multiple genotypes of GI and GII noroviruses. Coxsackieviruses A, astrovirus type 1 and adenovirus type 41 were prevalent. Norovirus GII/4 and coxsackievirus A24 isolates in wastewaters were closely related to respective outbreak strains isolated previously in Singapore. Conclusions: This study showed the widespread occurrence of all tested enteric virus groups in urban wastewaters. Genetic diversity of astroviruses, enteroviruses and noroviruses in the tropical region was observed. Significance and Impact of the Study: The high prevalence and great genetic diversity of human enteric viruses in urban wastewaters strongly supports the need of further comprehensive studies for evaluating the public health risk associated with viral pathogens in water environments.     

人星状病毒研究进展

人星状病毒是引起婴幼儿病毒性腹泻的重要病原之一,本文回顾了其分子生物学特征、型别、所致疾病特征、致病机制、流行病学特征及检测方法等方面的文献,指出人星状病毒目前存在多种型别,尤其是近年来发现了多种新型的星状病毒,主要的流行型别是血清型1型,不仅引起肠道内感染,还可以引起其他系统疾病,所致疾病特征和其他腹泻相关病毒相似,但是其致病机制的研究尚不清楚,对于几种新发现的星状病毒,其与疾病的相关性研究较少,且未建立完善的实验室检测方法,还需要进一步的研究。     

Gene conversion occurs within the mating-type locus of Cryptococcus neoformans during sexual reproduction

Meiotic recombination of sex chromosomes is thought to be repressed in organisms with heterogametic sex determination (e.g. mammalian X/Y chromosomes), due to extensive divergence and chromosomal rearrangements between the two chromosomes. However, proper segregation of sex chromosomes during meiosis requires crossing-over occurring within the pseudoautosomal regions (PAR). Recent studies reveal that recombination, in the form of gene conversion, is widely distributed within and may have played important roles in the evolution of some chromosomal regions within which recombination was thought to be repressed, such as the centromere cores of maize. Cryptococcus neoformans, a major human pathogenic fungus, has an unusually large mating-type locus (MAT, >100 kb), and the MAT alleles from the two opposite mating-types show extensive nucleotide sequence divergence and chromosomal rearrangements, mirroring characteristics of sex chromosomes. Meiotic recombination was assumed to be repressed within the C. neoformans MAT locus. A previous study identified recombination hot spots flanking the C. neoformans MAT, and these hot spots are associated with high GC content. Here, we investigated a GC-rich intergenic region located within the MAT locus of C. neoformans to establish if this region also exhibits unique recombination behavior during meiosis. Population genetics analysis of natural C. neoformans isolates revealed signals of homogenization spanning this GC-rich intergenic region within different C. neoformans lineages, consistent with a model in which gene conversion of this region during meiosis prevents it from diversifying within each lineage. By analyzing meiotic progeny from laboratory crosses, we found that meiotic recombination (gene conversion) occurs around the GC-rich intergenic region at a frequency equal to or greater than the meiotic recombination frequency observed in other genomic regions. We discuss the implications of these findings with regards to the possible functional and evolutionary importance of gene conversion within the C. neoformans MAT locus and, more generally, in fungi.     

Evolutionary strata on the X chromosomes of the dioecious plant Silene latifolia: evidence from new sex-linked genes

Despite its recent evolutionary origin, the sex chromosome system of the plant Silene latifolia shows signs of progressive suppression of recombination having created evolutionary strata of different X-Y divergence on sex chromosomes. However, even after 8 years of effort, this result is based on analyses of five sex-linked gene sequences, and the maximum divergence (and thus the age of this plant's sex chromosome system) has remained uncertain. More genes are therefore needed. Here, by segregation analysis of intron size variants (ISVS) and single nucleotide polymorphisms (SNPs), we identify three new Y-linked genes, one being duplicated on the Y chromosome, and test for evolutionary strata. All the new genes have homologs on the X and Y chromosomes. Synonymous divergence estimated between the X and Y homolog pairs is within the range of those already reported. Genetic mapping of the new X-linked loci shows that the map is the same in all three families that have been studied so far and that X-Y divergence increases with genetic distance from the pseudoautosomal region. We can now conclude that the divergence value is saturated, confirming the cessation of X-Y recombination in the evolution of the sex chromosomes at approximately 10-20 MYA.     

Size Variation of the Nonrecombining Region on the Mating-Type Chromosomes in the Fungal Podospora anserina Species Complex

Sex chromosomes often carry large nonrecombining regions that can extend progressively over time, generating evolutionary strata of sequence divergence. However, some sex chromosomes display an incomplete suppression of recombination. Large genomic regions without recombination and evolutionary strata have also been documented around fungal mating-type loci, but have been studied in only a few fungal systems. In the model fungus Podospora anserina (Ascomycota, Sordariomycetes), the reference S strain lacks recombination across a 0.8-Mb region around the mating-type locus. The lack of recombination in this region ensures that nuclei of opposite mating types are packaged into a single ascospore (pseudohomothallic lifecycle). We found evidence for a lack of recombination around the mating-type locus in the genomes of ten P. anserina strains and six closely related pseudohomothallic Podospora species. Importantly, the size of the nonrecombining region differed between strains and species, as indicated by the heterozygosity levels around the mating-type locus and experimental selfing. The nonrecombining region is probably labile and polymorphic, differing in size and precise location within and between species, resulting in occasional, but infrequent, recombination at a given base pair. This view is also supported by the low divergence between mating types, and the lack of strong linkage disequilibrium, chromosomal rearrangements, transspecific polymorphism and genomic degeneration. We found a pattern suggestive of evolutionary strata in P. pseudocomata. The observed heterozygosity levels indicate low but nonnull outcrossing rates in nature in these pseudohomothallic fungi. This study adds to our understanding of mating-type chromosome evolution and its relationship to mating systems.     

A new perspective on Listeria monocytogenes evolution

Listeria monocytogenes is a model organism for cellular microbiology and host-pathogen interaction studies and an important food-borne pathogen widespread in the environment, thus representing an attractive model to study the evolution of virulence. The phylogenetic structure of L. monocytogenes was determined by sequencing internal portions of seven housekeeping genes (3,288 nucleotides) in 360 representative isolates. Fifty-eight of the 126 disclosed sequence types were grouped into seven well-demarcated clonal complexes (clones) that comprised almost 75% of clinical isolates. Each clone had a unique or dominant serotype (4b for clones 1, 2 and 4, 1/2b for clones 3 and 5, 1/2a for clone 7, and 1/2c for clone 9), with no association of clones with clinical forms of human listeriosis. Homologous recombination was extremely limited (r/m<1 for nucleotides), implying long-term genetic stability of multilocus genotypes over time. Bayesian analysis based on 438 SNPs recovered the three previously defined lineages, plus one unclassified isolate of mixed ancestry. The phylogenetic distribution of serotypes indicated that serotype 4b evolved once from 1/2b, the likely ancestral serotype of lineage I. Serotype 1/2c derived once from 1/2a, with reference strain EGDe (1/2a) likely representing an intermediate evolutionary state. In contrast to housekeeping genes, the virulence factor internalin (InlA) evolved by localized recombination resulting in a mosaic pattern, with convergent evolution indicative of natural selection towards a truncation of InlA protein. This work provides a reference evolutionary framework for future studies on L. monocytogenes epidemiology, ecology, and virulence.