Intracellular Salmonella Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular S. Typhimurium.

Although Salmonella Typhimurium (STM) and Salmonella Paratyphi A (SPA) belong to the same phylogenetic species, share large portions of their genome and express many common virulence factors, they differ vastly in their host specificity, the immune response they elicit, and the clinical manifestations they cause. In this work, we compared their intracellular transcriptomic architecture and cellular phenotypes during human epithelial cell infection. While transcription induction of many metal transport systems, purines, biotin, PhoPQ and SPI-2 regulons was similar in both intracellular SPA and STM, we identified 234 differentially expressed genes that showed distinct expression patterns in intracellular SPA vs. STM. Surprisingly, clear expression differences were found in SPI-1, motility and chemotaxis, and carbon (mainly citrate, galactonate and ethanolamine) utilization pathways, indicating that these pathways are regulated differently during their intracellular phase. Concurring, on the cellular level, we show that while the majority of STM are non-motile and reside within Salmonella-Containing Vacuoles (SCV), a significant proportion of intracellular SPA cells are motile and compartmentalized in the cytosol. Moreover, we found that the elevated expression of SPI-1 and motility genes by intracellular SPA results in increased invasiveness of SPA, following exit from host cells. These findings demonstrate unexpected flagellum-dependent intracellular motility of a typhoidal Salmonella serovar and intriguing differences in intracellular localization between typhoidal and non-typhoidal salmonellae. We propose that these differences facilitate new cycles of host cell infection by SPA and may contribute to the ability of SPA to disseminate beyond the intestinal lamina propria of the human host during enteric fever.

A comparative transcriptomic approach to understanding the formation of cork.

KEY MESSAGE: The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.

Unraveling the evolution and coevolution of small regulatory RNAs and coding genes in Listeria.

BACKGROUND: Small regulatory RNAs (sRNAs) are widely found in bacteria and play key roles in many important physiological and adaptation processes. Studying their evolution and screening for events of coevolution with other genomic features is a powerful way to better understand their origin and assess a common functional or adaptive relationship between them. However, evolution and coevolution of sRNAs with coding genes have been sparsely investigated in bacterial pathogens. RESULTS: We designed a robust and generic phylogenomics approach that detects correlated evolution between sRNAs and protein-coding genes using their observed and inferred patterns of presence-absence in a set of annotated genomes. We applied this approach on 79 complete genomes of the Listeria genus and identified fifty-two accessory sRNAs, of which most were present in the Listeria common ancestor and lost during Listeria evolution. We detected significant coevolution between 23 sRNA and 52 coding genes and inferred the Listeria sRNA-coding genes coevolution network. We characterized a main hub of 12 sRNAs that coevolved with genes encoding cell wall proteins and virulence factors. Among them, an sRNA specific to L. monocytogenes species, rli133, coevolved with genes involved either in pathogenicity or in interaction with host cells, possibly acting as a direct negative post-transcriptional regulation. CONCLUSIONS: Our approach allowed the identification of candidate sRNAs potentially involved in pathogenicity and host interaction, consistent with recent findings on known pathogenicity actors. We highlight four sRNAs coevolving with seven internalin genes, some of which being important virulence factors in Listeria.

Links between Transcription, Environmental Adaptation and Gene Variability in Escherichia coli: Correlations between Gene Expression and Gene Variability Reflect Growth Efficiencies.

Gene expression is known to be the principle factor explaining how fast genes evolve. Highly transcribed genes evolve slowly because any negative impact caused by a particular mutation is magnified by protein abundance. However, gene expression is a phenotype that depends both on the environment and on the strains or species. We studied this phenotypic plasticity by analyzing the transcriptome profiles of four Escherichia coli strains grown in three different culture media, and explored how expression variability was linked to gene allelic diversity. Genes whose expression changed according to the media and not to the strains were less polymorphic than other genes. Genes for which transcription depended predominantly on the strain were more polymorphic than other genes and were involved in sensing and responding to environmental changes, with an overrepresentation of two-component system genes. Surprisingly, we found that the correlation between transcription and gene diversity was highly variable among growth conditions and could be used to quantify growth efficiency of a strain in a medium. Genetic variability was found to increase with gene expression in poor growth conditions. As such conditions are also characterized by down-regulation of all DNA repair systems, including transcription-coupled repair, we suggest that gene expression under stressful conditions may be mutagenic and thus leads to a variability in mutation rate among genes in the genome which contributes to the pattern of protein evolution.

The BioMart community portal: an innovative alternative to large, centralized data repositories.

The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

Development of an in vitro biofilm model for the study of the impact of fluoroquinolones on sewer biofilm microbiota.

Sewer biofilms are likely to constitute hotspots for selecting and accumulating antibiotic-resistant bacteria (ARB) and antibiotic resistance genes (ARGs). This study aimed to optimize culture conditions to obtain in vitro biofilms, mimicking the biofilm collected in sewers, to study the impact of fluoroquinolones (FQs) on sewer biofilm microbiota. Biofilms were grown on coupons in CDC Biofilm Reactors®, continuously fed with nutrients and inoculum (1/100 diluted wastewater). Different culture conditions were tested: (i) initial inoculum: diluted wastewater with or without sewer biofilm, (ii) coupon material: concrete vs. polycarbonate, and (iii) time of culture: 7 versus 14 days. This study found that the biomass was highest when in vitro biofilms were formed on concrete coupons. The biofilm taxonomic diversity was not affected by adding sewer biofilm to the initial inoculum nor by the coupon material. Pseudomonadales, Burkholderiales and Enterobacterales dominated in the sewer biofilm composition, whereas in vitro biofilms were mainly composed of Enterobacterales. The relative abundance of qnrA, B, D and S genes was higher in in vitro biofilms than sewer biofilm. The resistome of sewer biofilm showed the highest Shannon diversity index compared to wastewater and in vitro biofilms. A PCoA analysis showed differentiation of samples according to the nature of the sample, and a Procrustes analysis showed that the ARG changes observed were linked to changes in the microbial community. The following growing conditions were selected for in vitro biofilms: concrete coupons, initial inoculation with sewer biofilm, and a culture duration of 14 days. Then, biofilms were established under high and low concentrations of FQs to validate our in vitro biofilm model. Fluoroquinolone exposure had no significant impact on the abundance of qnr genes, but high concentration exposure increased the proportion of mutations in gyrA (codons S83L and D87N) and parC (codon S80I). In conclusion, this study allowed the determination of the culture conditions to develop an in vitro model of sewer biofilm; and was successfully used to investigate the impact of FQs on sewer microbiota. In the future, this setup could be used to clarify the role of sewer biofilms in disseminating resistance to FQs in the environment.

Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths.

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the approximately 18,000 families of orthologous genes, we found approximately 2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

Impact of environmental micropollutants and diet composition on the gut microbiota of wild european eels (Anguilla anguilla).

In fish, the gut microbiome plays a crucial role in homeostasis and health and is affected by several organic and inorganic environmental contaminants. Amphidromous fish are sentinel species, particularly exposed to these stressors. We used whole metagenome sequencing to characterize the gut microbiome of wild European eels (Anguilla anguilla) at a juvenile stage captured from three sites with contrasted pollution levels in term of heavy metals and persistent organic pollutants. The objectives were to identify what parameters could alter the gut microbiome of this catadromous fish and to explore the potential use of microbiota as bioindicators of environment quality. We identified a total of 1079 microbial genera. Overall, gut microbiome was dominated by Proteobacteria, Firmicutes and Actinobacteria. Alpha and beta diversity were different amongst sites and could be explained by a reduced number of environmental and biological factors, specifically the relative abundance of fish preys in eels' diet, PCB101, γHCH (lindane), transnonachlor and arsenic. Furthermore, we identified a series of indicator taxa with differential abundance between the three sites. Changes in the microbial communities in the gut caused by environmental pollutants were previously undocumented in European eels. Our results indicate that microbiota might represent another route by which pollutants affect the health of these aquatic sentinel organisms.

Phylodynamic study of the conserved RNA structure encompassing the hemagglutinin cleavage site encoding region of H5 and H7 low pathogenic avian influenza viruses.

Highly pathogenic avian influenza viruses (HPAIVs) evolve from low pathogenic avian influenza viruses (LPAIVs) of the H5 and H7 subtypes. This evolution is characterized by the acquisition of a multi-basic cleavage site (MBCS) motif in the hemagglutinin (HA) that leads to an extended viral tropism and severe disease in poultry. One key unanswered question is whether the risk of transition to HPAIVs is similar for all LPAIVs H5 or H7 strains, or whether specific determinants in the HA sequence of some H5 or H7 LPAIV strains correlate with a higher risk of transition to HPAIVs. Here, we determined if specific features of the conserved RNA stem-loop located at the HA cleavage site-encoding region could be detected along the LPAIV to HPAIV evolutionary pathway. Analysis of the thermodynamic stability of the predicted RNA structures showed no specific patterns common to HA sequences leading to HPAIVs and distinct from those remaining LPAIVs. However, RNA structure clustering analysis revealed that most of the American lineage ancestors leading to H7 emergences via recombination shared the same viral RNA (vRNA) structure topology at the HA1/HA2 boundary region. Our study thus identified predicted secondary RNA structures present in the HA of H7 viruses, which could promote genetic recombination and acquisition of a multibasic cleavage site motif (MBCS).

Widely rhythmic transcriptome in Calanus finmarchicus during the high Arctic summer solstice period.

Solar light/dark cycles and seasonal photoperiods underpin daily and annual rhythms of life on Earth. Yet, the Arctic is characterized by several months of permanent illumination ("midnight sun"). To determine the persistence of 24h rhythms during the midnight sun, we investigated transcriptomic dynamics in the copepod Calanus finmarchicus during the summer solstice period in the Arctic, with the lowest diel oscillation and the highest altitude of the sun's position. Here we reveal that in these extreme photic conditions, a widely rhythmic daily transcriptome exists, showing that very weak solar cues are sufficient to entrain organisms. Furthermore, at extremely high latitudes and under sea-ice, gene oscillations become re-organized to include <24h rhythms. Environmental synchronization may therefore be modulated to include non-photic signals (i.e. tidal cycles). The ability of zooplankton to be synchronized by extremely weak diel and potentially tidal cycles, may confer an adaptive temporal reorganization of biological processes at high latitudes.

Daily transcriptomes of the copepod Calanus finmarchicus during the summer solstice at high Arctic latitudes.

The zooplankter Calanus finmarchicus is a member of the so-called "Calanus Complex", a group of copepods that constitutes a key element of the Arctic polar marine ecosystem, providing a crucial link between primary production and higher trophic levels. Climate change induces the shift of C. finmarchicus to higher latitudes with currently unknown impacts on its endogenous timing. Here we generated a daily transcriptome of C. finmarchicus at two high Arctic stations, during the more extreme time of Midnight Sun, the summer solstice. While the southern station (74.5 °N) was sea ice-free, the northern one (82.5 °N) was sea ice-covered. The mRNAs of the 42 samples have been sequenced with an average of 126 ± 5 million reads (mean ± SE) per sample, and aligned to the reference transcriptome. We detail the quality assessment of the datasets and the complete annotation procedure, providing the possibility to investigate daily gene expression of this ecologically important species at high Arctic latitudes, and to compare gene expression according to latitude and sea ice-coverage.

A module located at a chromosomal integration hot spot is responsible for the multidrug resistance of a reference strain from Escherichia coli clonal group A.

Escherichia coli clonal group A (CGA) commonly exhibits a distinctive multidrug antimicrobial resistance phenotype-i.e., resistance to ampicillin, chloramphenicol, streptomycin, sulfonamides, tetracycline, and trimethoprim (ACSSuTTp)-and has accounted for up to 50% of trimethoprim-sulfamethoxazole-resistant E. coli urinary tract infections in some locales. Annotation of the whole-genome sequencing of UMN026, a reference CGA strain, clarified the genetic basis for this strain's ACSSuTTp antimicrobial resistance phenotype. Most of the responsible genes were clustered in a unique 23-kbp chromosomal region, designated the genomic resistance module (GRM), which occurred within a 105-kbp genomic island situated at the leuX tRNA. The GRM is characterized by numerous remnants of mobilization and rearrangement events suggesting multiple horizontal transfers. Additionally, comparative genomic analysis of the leuX tRNA genomic island in 14 sequenced E. coli genomes showed that this region is a hot spot of integration, with the presence/absence of specific subregions being uncorrelated with either the phylogenetic group or the pathotype. Our data illustrate the importance of whole-genome sequencing in the detection of genetic elements involved in antimicrobial resistance. Additionally, this is the first documentation of the bla(TEM) and dhfrVII genes in a chromosomal location in E. coli strains.

aes, the gene encoding the esterase B in Escherichia coli, is a powerful phylogenetic marker of the species.

BACKGROUND: Previous studies have established a correlation between electrophoretic polymorphism of esterase B, and virulence and phylogeny of Escherichia coli. Strains belonging to the phylogenetic group B2 are more frequently implicated in extraintestinal infections and include esterase B2 variants, whereas phylogenetic groups A, B1 and D contain less virulent strains and include esterase B1 variants. We investigated esterase B as a marker of phylogeny and/or virulence, in a thorough analysis of the esterase B-encoding gene. RESULTS: We identified the gene encoding esterase B as the acetyl-esterase gene (aes) using gene disruption. The analysis of aes nucleotide sequences in a panel of 78 reference strains, including the E. coli reference (ECOR) strains, demonstrated that the gene is under purifying selection. The phylogenetic tree reconstructed from aes sequences showed a strong correlation with the species phylogenetic history, based on multi-locus sequence typing using six housekeeping genes. The unambiguous distinction between variants B1 and B2 by electrophoresis was consistent with Aes amino-acid sequence analysis and protein modelling, which showed that substituted amino acids in the two esterase B variants occurred mostly at different sites on the protein surface. Studies in an experimental mouse model of septicaemia using mutant strains did not reveal a direct link between aes and extraintestinal virulence. Moreover, we did not find any genes in the chromosomal region of aes to be associated with virulence. CONCLUSION: Our findings suggest that aes does not play a direct role in the virulence of E. coli extraintestinal infection. However, this gene acts as a powerful marker of phylogeny, illustrating the extensive divergence of B2 phylogenetic group strains from the rest of the species.

Selection acts on DNA secondary structures to decrease transcriptional mutagenesis.

Single-stranded DNA is more subject to mutation than double stranded. During transcription, DNA is transiently single stranded and therefore subject to higher mutagenesis. However, if local intra-strand secondary structures are formed, some bases will be paired and therefore less sensitive to mutation than unpaired bases. Using complete genome sequences of Escherichia coli, we show that local intra-strand secondary structures can, as a consequence, be used to define an index of transcription-driven mutability. At gene level, we show that natural selection has favoured a reduced transcription-driven mutagenesis via the higher than expected frequency of occurrence of intra-strand secondary structures. Such selection is stronger in highly expressed genes and suggests a sequence-dependent way to control mutation rates and a novel form of selection affecting the evolution of synonymous mutations.

Human Adipose-Derived Stem Cells Expanded Under Ambient Oxygen Concentration Accumulate Oxidative DNA Lesions and Experience Procarcinogenic DNA Replication Stress.

Adipose-derived stem cells (ADSCs) have led to growing interest in cell-based therapy because they can be easily harvested from an abundant tissue. ADSCs must be expanded in vitro before transplantation. This essential step causes concerns about the safety of adult stem cells in terms of potential transformation. Tumorigenesis is driven in its earliest step by DNA replication stress, which is characterized by the accumulation of stalled DNA replication forks and activation of the DNA damage response. Thus, to evaluate the safety of ADSCs during ex vivo expansion, we monitored DNA replication under atmospheric (21%) or physiologic (1%) oxygen concentration. Here, by combining immunofluorescence and DNA combing, we show that ADSCs cultured under 21% oxygen accumulate endogenous oxidative DNA lesions, which interfere with DNA replication by increasing fork stalling events, thereby leading to incomplete DNA replication and fork collapse. Moreover, we found by RNA sequencing (RNA-seq) that culture of ADSCs under atmospheric oxygen concentration leads to misexpression of cell cycle and DNA replication genes, which could contribute to DNA replication stress. Finally, analysis of acquired small nucleotide polymorphism shows that expansion of ADSCs under 21% oxygen induces a mutational bias toward deleterious transversions. Overall, our results suggest that expanding ADSCs at a low oxygen concentration could reduce the risk for DNA replication stress-associated transformation, as occurs in neoplastic tissues. Stem Cells Translational Medicine 2017;6:68-76.

Listeriomics: an Interactive Web Platform for Systems Biology of Listeria.

As for many model organisms, the amount of Listeria omics data produced has recently increased exponentially. There are now >80 published complete Listeria genomes, around 350 different transcriptomic data sets, and 25 proteomic data sets available. The analysis of these data sets through a systems biology approach and the generation of tools for biologists to browse these various data are a challenge for bioinformaticians. We have developed a web-based platform, named Listeriomics, that integrates different tools for omics data analyses, i.e., (i) an interactive genome viewer to display gene expression arrays, tiling arrays, and sequencing data sets along with proteomics and genomics data sets; (ii) an expression and protein atlas that connects every gene, small RNA, antisense RNA, or protein with the most relevant omics data; (iii) a specific tool for exploring protein conservation through the Listeria phylogenomic tree; and (iv) a coexpression network tool for the discovery of potential new regulations. Our platform integrates all the complete Listeria species genomes, transcriptomes, and proteomes published to date. This website allows navigation among all these data sets with enriched metadata in a user-friendly format and can be used as a central database for systems biology analysis. IMPORTANCE In the last decades, Listeria has become a key model organism for the study of host-pathogen interactions, noncoding RNA regulation, and bacterial adaptation to stress. To study these mechanisms, several genomics, transcriptomics, and proteomics data sets have been produced. We have developed Listeriomics, an interactive web platform to browse and correlate these heterogeneous sources of information. Our website will allow listeriologists and microbiologists to decipher key regulation mechanism by using a systems biology approach.

Remodeling of the cycling transcriptome of the oyster Crassostrea gigas by the harmful algae Alexandrium minutum.

As a marine organism, the oyster Crassostrea gigas inhabits a complex biotope governed by interactions between the moon and the sun cycles. We used next-generation sequencing to investigate temporal regulation of oysters under light/dark entrainment and the impact of harmful algal exposure. We found that ≈6% of the gills' transcriptome exhibits circadian expression, characterized by a nocturnal and bimodal pattern. Surprisingly, a higher number of ultradian transcripts were also detected under solely circadian entrainment. The results showed that a bloom of Alexandrium minutum generated a remodeling of the bivalve's temporal structure, characterized by a loss of oscillations, a genesis of de novo oscillating transcripts, and a switch in the period of oscillations. These findings provide unprecedented insights into the diurnal landscape of the oyster's transcriptome and pleiotropic remodeling due to toxic algae exposure, revealing the intrinsic plasticity of the cycling transcriptome in oysters.

RNAbrowse: RNA-Seq de novo assembly results browser.

Transcriptome analysis based on a de novo assembly of next generation RNA sequences is now performed routinely in many laboratories. The generated results, including contig sequences, quantification figures, functional annotations and variation discovery outputs are usually bulky and quite diverse. This article presents a user oriented storage and visualisation environment permitting to explore the data in a top-down manner, going from general graphical views to all possible details. The software package is based on biomart, easy to install and populate with local data. The software package is available under the GNU General Public License (GPL) at http://bioinfo.genotoul.fr/RNAbrowse.

PASTEC: an automatic transposable element classification tool.

SUMMARY: The classification of transposable elements (TEs) is key step towards deciphering their potential impact on the genome. However, this process is often based on manual sequence inspection by TE experts. With the wealth of genomic sequences now available, this task requires automation, making it accessible to most scientists. We propose a new tool, PASTEC, which classifies TEs by searching for structural features and similarities. This tool outperforms currently available software for TE classification. The main innovation of PASTEC is the search for HMM profiles, which is useful for inferring the classification of unknown TE on the basis of conserved functional domains of the proteins. In addition, PASTEC is the only tool providing an exhaustive spectrum of possible classifications to the order level of the Wicker hierarchical TE classification system. It can also automatically classify other repeated elements, such as SSR (Simple Sequence Repeats), rDNA or potential repeated host genes. Finally, the output of this new tool is designed to facilitate manual curation by providing to biologists with all the evidence accumulated for each TE consensus. AVAILABILITY: PASTEC is available as a REPET module or standalone software (http://urgi.versailles.inra.fr/download/repet/REPET_linux-x64-2.2.tar.gz). It requires a Unix-like system. There are two standalone versions: one of which is parallelized (requiring Sun grid Engine or Torque), and the other of which is not.

Animal and human pathogenic Escherichia coli strains share common genetic backgrounds.

Escherichia coli is a versatile species encompassing both commensals of the digestive tracts of many vertebrates, including humans, and pathogenic strains causing various intra- and extraintestinal infections. Despite extensive gene flow between strains, the E. coli species has a globally clonal population structure, consisting of distinct phylogenetic groups. Little is known about the relationships between phylogenetic groups and host specificity. We therefore used multilocus sequence typing (MLST) to investigate phylogenetic relationships and evaluated the virulence gene content of 35 E. coli strains representative of the diverse diseases encountered in domestic animals. We compared these strains with a panel of 101 human pathogenic and 98 non-human and human commensal strains representative of the phylogenetic and pathovar diversity of this species. A global factorial analysis of correspondence indicated that extraintestinal infections were caused mostly by phylogenetic group B2 strains, whereas intraintestinal infections were caused mostly by phylogenetic group A/B1/E strains, with strains responsible from extraintestinal or intraintestinal infections having specific virulence factors. It was not possible to distinguish between strains of human and animal origin. A detailed phylogenetic analysis of the MLST data showed that numerous pathogenic animal and human strains are very closely related, and had a number of virulence genes in common. However, a set of specific adhesins was identified in animal non-B2 group strains of all pathotypes. In conclusion, human and animal pathogenic strains share common genetic backgrounds, but non-B2 strains of different origins seem to have different sets of adhesins that could be involved in host specificity.