Recombinant production of a diffusible signal factor inhibits Salmonella invasion and animal carriage.

The complex chemical environment of the intestine is defined largely by the metabolic products of the resident microbiota. Enteric pathogens, elegantly evolved to thrive in the gut, use these chemical products as signals to recognize specific niches and to promote their survival and virulence. Our previous work has shown that a specific class of quorum-sensing molecules found within the gut, termed diffusible signal factors (DSF), signals the repression of Salmonella tissue invasion, thus defining a means by which this pathogen recognizes its location and modulates virulence to optimize its survival. Here, we determined whether the recombinant production of a DSF could reduce Salmonella virulence in vitro and in vivo. We found that the most potent repressor of Salmonella invasion, cis-2-hexadecenoic acid (c2-HDA), could be recombinantly produced in E. coli by the addition of a single exogenous gene encoding a fatty acid enoyl-CoA dehydratase/thioesterase and that co-culture of the recombinant strain with Salmonella potently inhibited tissue invasion by repressing Salmonella genes required for this essential virulence function. Using the well characterized E. coli Nissle 1917 strain and a chicken infection model, we found that the recombinant DSF-producing strain could be stably maintained in the large intestine. Further, challenge studies demonstrated that this recombinant organism could significantly reduce Salmonella colonization of the cecum, the site of carriage in this animal species. These findings thus describe a plausible means by which Salmonella virulence may be affected in animals by in situ chemical manipulation of functions essential for colonization and virulence.

Despite our best efforts, infections of agricultural animals with Salmonella persist, posing threats to food safety. Few, if any, measures have proven effective in reducing Salmonella carriage in animals used for food, a major source of this pathogen. Antibiotics are ineffective at curtailing infection and have served only to exacerbate the global crisis of antimicrobial resistance. The alternative then is to seek novel means to reduce Salmonella disease and carriage by preventing its colonization of livestock and poultry. Here we describe an approach targeting invasion, a function essential for Salmonella carriage and disease in animals. We show that a potent chemical inhibitor of invasion, the diffusible signal factor cis-2 hexadecenoic acid, can be produced by recombinant E. coli strains capable of stably colonizing the animal intestine, providing a means to directly affect the virulence of Salmonella within an animal host. These studies may thus provide a route to reduce the carriage of this pathogen in production animals and thus the spread of disease to humans.

Salmonella Invasion Is Controlled by Competition among Intestinal Chemical Signals.

The intestine is a complex, ever-changing environment replete with an array of signaling molecules. To colonize such a complex organ, pathogens have adapted to utilize specific cues from the local environment to intricately regulate the expression of their virulence determinants. Salmonella preferentially colonizes the distal ileum, a niche enriched in the metabolite formic acid. Here, we show that the relatively higher concentration of this metabolite in the distal ileum prevents other signals from repressing Salmonella invasion in that region. We show that imported and unmetabolized formic acid functions as a cytoplasmic signal that competitively binds to HilD, the master transcriptional regulator of Salmonella invasion, thus preventing repressive fatty acids from binding to the protein. This results in an increased lifetime of HilD and subsequent derepression of invasion genes. This study demonstrates an important mechanism by which Salmonella utilizes competition among signals in the gut to its advantage as a pathogen. IMPORTANCE Enteric pathogens acutely sense their environment for signals to regulate their virulence functions. We demonstrate here that the enteric pathogen Salmonella utilizes the competition among certain regional intestinal constituents to modulate its virulence determinants in that region. We show that the high concentration of formic acid in the ileum outcompetes other signals and triggers the activation of virulence genes in the ileum. This study shows a delicate spatial and temporal mechanism by which enteric pathogens may utilize the competition among environmental cues to optimize their pathogenicity.

Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead.

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

AraC-type regulators HilC and RtsA are directly controlled by an intestinal fatty acid to regulate Salmonella invasion.

Invasion of the intestinal epithelium is an essential but energetically expensive survival strategy and is, therefore, tightly regulated by using specific cues from the environment. The enteric pathogen Salmonella controls its invasion machinery through the elegant coordination of three AraC-type transcription activators, HilD, HilC, and RtsA. Most environmental signals target HilD to control invasion, whereas HilC and RtsA are known only to augment these effects on HilD. Here we show that a fatty acid found in the murine colon, cis-2-hexadecenoic acid (c2-HDA), represses Salmonella invasion by directly targeting HilC and RtsA, in addition to HilD. c2-HDA directly binds each of these regulators and inhibits their attachment to DNA targets, repressing invasion even in the absence of HilD. Fatty acid binding, however, does not affect HilC and RtsA protein stability, unlike HilD. Importantly, we show that HilC and RtsA are highly effective in restoring HilD production and invasion gene expression after elimination of the repressive fatty acid c2-HDA. Together, these results illuminate a precise mechanism by which HilC and RtsA may modulate invasion as Salmonella navigates through different regions of the intestine, contributing to our understanding of how this enteric pathogen senses and adapts to a diverse intestinal environment while maintaining its virulence.

Genomic Assessment of Global Population Structure in a Highly Migratory and Habitat Versatile Apex Predator, the Tiger Shark (Galeocerdo cuvier).

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST = 0.125-0.144) and candidate outlier loci (FST = 0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared with sharks from the Indian Ocean (FST = 0.003-0.005, P < 0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Who infects whom?-Reconstructing infection chains of Mycobacterium avium ssp. paratuberculosis in an endemically infected dairy herd by use of genomic data.

Recent evidence of circulation of multiple strains within herds and mixed infections of cows marks the beginning of a rethink of our knowledge on Mycobacterium avium ssp. paratuberculosis (MAP) epidemiology. Strain typing opens new ways to investigate MAP transmission. This work presents a method for reconstructing infection chains in a setting of endemic Johne's disease on a well-managed dairy farm. By linking genomic data with demographic field data, strain-specific differences in spreading patterns could be quantified for a densely sampled dairy herd. Mixed infections of dairy cows with MAP are common, and some strains spread more successfully. Infected cows remain susceptible for co-infections with other MAP genotypes. The model suggested that cows acquired infection from 1-4 other cows and spread infection to 0-17 individuals. Reconstructed infection chains supported the hypothesis that high shedding animals that started to shed at an early age and showed a progressive infection pattern represented a greater risk for spreading MAP. Transmission of more than one genotype between animals was recorded. In this farm with a good MAP control management program, adult-to-adult contact was proposed as the most important transmission route to explain the reconstructed networks. For each isolate, at least one more likely ancestor could be inferred. Our study results help to capture underlying transmission processes and to understand the challenges of tracing MAP spread within a herd. Only the combination of precise longitudinal field data and bacterial strain type information made it possible to trace infection in such detail.

A diffusible signal factor of the intestine dictates Salmonella invasion through its direct control of the virulence activator HilD.

Successful intestinal infection by Salmonella requires optimized invasion of the gut epithelium, a function that is energetically costly. Salmonella have therefore evolved to intricately regulate the expression of their virulence determinants by utilizing specific environmental cues. Here we show that a powerful repressor of Salmonella invasion, a cis-2 unsaturated long chain fatty acid, is present in the murine large intestine. Originally identified in Xylella fastidiosa as a diffusible signal factor for quorum sensing, this fatty acid directly interacts with HilD, the master transcriptional regulator of Salmonella, and prevents hilA activation, thus inhibiting Salmonella invasion. We further identify the fatty acid binding region of HilD and show it to be selective and biased in favour of signal factors with a cis-2 unsaturation over other intestinal fatty acids. Single mutation of specific HilD amino acids to alanine prevented fatty acid binding, thereby alleviating their repressive effect on invasion. Together, these results highlight an exceedingly sensitive mechanism used by Salmonella to colonize its host by detecting and exploiting specific molecules present within the complex intestinal environment.

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Population Gene Introgression and High Genome Plasticity for the Zoonotic Pathogen Streptococcus agalactiae.

The influence that bacterial adaptation (or niche partitioning) within species has on gene spillover and transmission among bacterial populations occupying different niches is not well understood. Streptococcus agalactiae is an important bacterial pathogen that has a taxonomically diverse host range making it an excellent model system to study these processes. Here, we analyze a global set of 901 genome sequences from nine diverse host species to advance our understanding of these processes. Bayesian clustering analysis delineated 12 major populations that closely aligned with niches. Comparative genomics revealed extensive gene gain/loss among populations and a large pan genome of 9,527 genes, which remained open and was strongly partitioned among niches. As a result, the biochemical characteristics of 11 populations were highly distinctive (significantly enriched). Positive selection was detected and biochemical characteristics of the dispensable genes under selection were enriched in ten populations. Despite the strong gene partitioning, phylogenomics detected gene spillover. In particular, tetracycline resistance (which likely evolved in the human-associated population) from humans to bovine, canines, seals, and fish, demonstrating how a gene selected in one host can ultimately be transmitted into another, and biased transmission from humans to bovines was confirmed with a Bayesian migration analysis. Our findings show high bacterial genome plasticity acting in balance with selection pressure from distinct functional requirements of niches that is associated with an extensive and highly partitioned dispensable genome, likely facilitating continued and expansive adaptation.

Salmonella invasion is controlled through the secondary structure of the hilD transcript.

Virulence functions of bacterial pathogens are often energetically costly and thus are subjected to intricate regulatory mechanisms. In Salmonella, invasion of the intestinal epithelium, an essential early step in virulence, requires the production of a multi-protein type III secretion apparatus. The pathogen mitigates the overall cost of invasion by inducing it in only a fraction of its population. This constitutes a successful virulence strategy as invasion by a small number is sufficient to promote the proliferation of the non-invading majority. Such a system suggests the existence of a sensitive triggering mechanism that permits only a minority of Salmonella to reach a threshold of invasion-gene induction. We show here that the secondary structure of the invasion regulator hilD message provides such a trigger. The 5' end of the hilD mRNA is predicted to contain two mutually exclusive stem-loop structures, the first of which (SL1) overlaps the ribosome-binding site and the ORF start codon. Changes that reduce its stability enhance invasion gene expression, while those that increase stability reduce invasion. Conversely, disrupting the second stem-loop (SL2) represses invasion genes. Although SL2 is the energetically more favorable, repression through SL1 is enhanced by binding of the global regulator CsrA. This system thus alters the levels of hilD mRNA and is so sensitive that changing a single base pair within SL1, predicted to augment its stability, eliminates expression of invasion genes and significantly reduces Salmonella virulence in mice. This system thus provides a possible means to rapidly and finely tune an essential virulence function.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy.

BACKGROUND: Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy. RESULTS: Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism. CONCLUSIONS: This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Mitochondrial genome of an Atlantic white shark (Carcharodon carcharias).

Here we report the first full mitochondrial genome sequence for a white shark caught in the Atlantic Ocean. The mitochondrial genome is 16,745 bp in length and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding control region. The base composition of this mtDNA lineage is A: 30.7%, C: 26.9%, G: 13.8%, and T: 28.6%. In concordance with previous population genetic studies, the Atlantic caught individual forms a separate lineage from individuals caught on either side of the Pacific Ocean.

Genome based phylogeny and comparative genomic analysis of intra-mammary pathogenic Escherichia coli.

Escherichia coli is an important cause of bovine mastitis and can cause both severe inflammation with a short-term transient infection, as well as less severe, but more chronic inflammation and infection persistence. E. coli is a highly diverse organism that has been classified into a number of different pathotypes or pathovars, and mammary pathogenic E. coli (MPEC) has been proposed as a new such pathotype. The purpose of this study was to use genome sequence data derived from both transient and persistent MPEC isolates (two isolates of each phenotype) to construct a genome-based phylogeny that places MPEC in its phylogenetic context with other E. coli pathovars. A subsidiary goal was to conduct comparative genomic analyses of these MPEC isolates with other E. coli pathovars to provide a preliminary perspective on loci that might be correlated with the MPEC phenotype. Both concatenated and consensus tree phylogenies did not support MPEC monophyly or the monophyly of either transient or persistent phenotypes. Three of the MPEC isolates (ECA-727, ECC-Z, and ECA-O157) originated from within the predominately commensal clade of E. coli, referred to as phylogroup A. The fourth MPEC isolate, of the persistent phenotype (ECC-1470), was sister group to an isolate of ETEC, falling within the E. coli B1 clade. This suggests that the MPEC phenotype has arisen on numerous independent occasions and that this has often, although not invariably, occurred from commensal ancestry. Examination of the genes present in the MPEC strains relative to the commensal strains identified a consistent presence of the type VI secretion system (T6SS) in the MPEC strains, with only occasional representation in commensal strains, suggesting that T6SS may be associated with MPEC pathogenesis and/or as an inter-bacterial competitive attribute and therefore could represent a useful target to explore for the development of MPEC specific inhibitors.

Phylogenomics and the dynamic genome evolution of the genus Streptococcus.

The genus Streptococcus comprises important pathogens that have a severe impact on human health and are responsible for substantial economic losses to agriculture. Here, we utilize 46 Streptococcus genome sequences (44 species), including eight species sequenced here, to provide the first genomic level insight into the evolutionary history and genetic basis underlying the functional diversity of all major groups of this genus. Gene gain/loss analysis revealed a dynamic pattern of genome evolution characterized by an initial period of gene gain followed by a period of loss, as the major groups within the genus diversified. This was followed by a period of genome expansion associated with the origins of the present extant species. The pattern is concordant with an emerging view that genomes evolve through a dynamic process of expansion and streamlining. A large proportion of the pan-genome has experienced lateral gene transfer (LGT) with causative factors, such as relatedness and shared environment, operating over different evolutionary scales. Multiple gene ontology terms were significantly enriched for each group, and mapping terms onto the phylogeny showed that those corresponding to genes born on branches leading to the major groups represented approximately one-fifth of those enriched. Furthermore, despite the extensive LGT, several biochemical characteristics have been retained since group formation, suggesting genomic cohesiveness through time, and that these characteristics may be fundamental to each group. For example, proteolysis: mitis group; urea metabolism: salivarius group; carbohydrate metabolism: pyogenic group; and transcription regulation: bovis group.

Transcriptomic and genomic evidence for Streptococcus agalactiae adaptation to the bovine environment.

BACKGROUND: Streptococcus agalactiae is a major cause of bovine mastitis, which is the dominant health disorder affecting milk production within the dairy industry and is responsible for substantial financial losses to the industry worldwide. However, there is considerable evidence for host adaptation (ecotypes) within S. agalactiae, with both bovine and human sourced isolates showing a high degree of distinctiveness, suggesting differing ability to cause mastitis. Here, we (i) generate RNAseq data from three S. agalactiae isolates (two putative bovine adapted and one human) and (ii) compare publicly available whole genome shotgun sequence data from an additional 202 isolates, obtained from six host species, to elucidate possible genetic factors/adaptations likely important for S. agalactiae growth and survival in the bovine mammary gland. RESULTS: Tests for differential expression showed distinct expression profiles for the three isolates when grown in bovine milk. A key finding for the two putatively bovine adapted isolates was the up regulation of a lactose metabolism operon (Lac.2) that was strongly correlated with the bovine environment (all 36 bovine sourced isolates on GenBank possessed the operon, in contrast to only 8/151 human sourced isolates). Multi locus sequence typing of all genome sequences and phylogenetic analysis using conserved operon genes from 44 S. agalactiae isolates and 16 additional Streptococcus species provided strong evidence for acquisition of the operon via multiple lateral gene transfer events, with all Streptococcus species known to be major causes of mastitis, identified as possible donors. Furthermore, lactose fermentation tests were only positive for isolates possessing Lac.2. Combined, these findings suggest that lactose metabolism is likely an important adaptation to the bovine environment. Additional up regulation in the bovine adapted isolates included genes involved in copper homeostasis, metabolism of purine, pyrimidine, glycerol and glucose, and possibly aminoglycoside antibiotic resistance. CONCLUSION: We detected several genetic factors likely important in S. agalactiae's adaptation to the bovine environment, in particular lactose metabolism. Of concern is the up regulation of a putative antibiotic resistance gene (GCN5-related N-acetyltransferase) that might reflect an adaptation to the use of aminoglycoside antibiotics within this environment.

Comparative characterization of the virulence gene clusters (lipooligosaccharide [LOS] and capsular polysaccharide [CPS]) for Campylobacter coli, Campylobacter jejuni subsp. jejuni and related Campylobacter species.

Campylobacter jejuni subsp. jejuni and Campylobacter coli are leading causes of gastroenteritis, with virulence linked to cell surface carbohydrate diversity. Although the associated gene clusters are well studied for C. jejuni subsp. jejuni, C. coli has been largely neglected. Here we provide comparative analysis of the lipooligosaccharide (LOS) and capsular polysaccharide (CPS) gene clusters, using genome and cluster sequence data for 36 C. coli strains, 67 C. jejuni subsp. jejuni strains and ten additional Campylobacter species. Similar to C. jejuni subsp. jejuni, C. coli showed high LOS/CPS gene diversity, with each cluster delineated into eight gene content classes. This diversity was predominantly due to extensive gene gain/loss, with the lateral transfer of genes likely occurring both within and between species and also between the LOS and CPS. Additional mechanisms responsible for LOS/CPS diversity included phase-variable homopolymeric repeats, gene duplication/inactivation, and possibly host environment selection pressure. Analyses also showed that (i) strains of C. coli and Campylobacter upsaliensis possessed genes homologous to the sialic acid genes implicated in the neurological disorder Guillain-Barré syndrome (GBS), and (ii) C. coli LOS classes were differentiated between bovine and poultry hosts, potentially aiding post infection source tracking.

Genome characterization and population genetic structure of the zoonotic pathogen, Streptococcus canis.

BACKGROUND: Streptococcus canis is an important opportunistic pathogen of dogs and cats that can also infect a wide range of additional mammals including cows where it can cause mastitis. It is also an emerging human pathogen. RESULTS: Here we provide characterization of the first genome sequence for this species, strain FSL S3-227 (milk isolate from a cow with an intra-mammary infection). A diverse array of putative virulence factors was encoded by the S. canis FSL S3-227 genome. Approximately 75% of these gene sequences were homologous to known Streptococcal virulence factors involved in invasion, evasion, and colonization. Present in the genome are multiple potentially mobile genetic elements (MGEs) [plasmid, phage, integrative conjugative element (ICE)] and comparison to other species provided convincing evidence for lateral gene transfer (LGT) between S. canis and two additional bovine mastitis causing pathogens (Streptococcus agalactiae, and Streptococcus dysgalactiae subsp. dysgalactiae), with this transfer possibly contributing to host adaptation. Population structure among isolates obtained from Europe and USA [bovine = 56, canine = 26, and feline = 1] was explored. Ribotyping of all isolates and multi locus sequence typing (MLST) of a subset of the isolates (n = 45) detected significant differentiation between bovine and canine isolates (Fisher exact test: P = 0.0000 [ribotypes], P = 0.0030 [sequence types]), suggesting possible host adaptation of some genotypes. Concurrently, the ancestral clonal complex (54% of isolates) occurred in many tissue types, all hosts, and all geographic locations suggesting the possibility of a wide and diverse niche. CONCLUSION: This study provides evidence highlighting the importance of LGT in the evolution of the bacteria S. canis, specifically, its possible role in host adaptation and acquisition of virulence factors. Furthermore, recent LGT detected between S. canis and human bacteria (Streptococcus urinalis) is cause for concern, as it highlights the possibility for continued acquisition of human virulence factors for this emerging zoonotic pathogen.

Gene repertoire evolution of Streptococcus pyogenes inferred from phylogenomic analysis with Streptococcus canis and Streptococcus dysgalactiae.

Streptococcus pyogenes, is an important human pathogen classified within the pyogenic group of streptococci, exclusively adapted to the human host. Our goal was to employ a comparative evolutionary approach to better understand the genomic events concomitant with S. pyogenes human adaptation. As part of ascertaining these events, we sequenced the genome of one of the potential sister species, the agricultural pathogen S. canis, and combined it in a comparative genomics reconciliation analysis with two other closely related species, Streptococcus dysgalactiae and Streptococcus equi, to determine the genes that were gained and lost during S. pyogenes evolution. Genome wide phylogenetic analyses involving 15 Streptococcus species provided convincing support for a clade of S. equi, S. pyogenes, S. dysgalactiae, and S. canis and suggested that the most likely S. pyogenes sister species was S. dysgalactiae. The reconciliation analysis identified 113 genes that were gained on the lineage leading to S. pyogenes. Almost half (46%) of these gained genes were phage associated and 14 showed significant matches to experimentally verified bacteria virulence factors. Subsequent to the origin of S. pyogenes, over half of the phage associated genes were involved in 90 different LGT events, mostly involving different strains of S. pyogenes, but with a high proportion involving the horse specific pathogen S. equi subsp. equi, with the directionality almost exclusively (86%) in the S. pyogenes to S. equi direction. Streptococcus agalactiae appears to have played an important role in the evolution of S. pyogenes with a high proportion of LGTs originating from this species. Overall the analysis suggests that S. pyogenes adaptation to the human host was achieved in part by (i) the integration of new virulence factors (e.g. speB, and the sal locus) and (ii) the construction of new regulation networks (e.g. rgg, and to some extent speB).

Comparative genomic analysis of the Streptococcus dysgalactiae species group: gene content, molecular adaptation, and promoter evolution.

Comparative genomics of closely related bacterial species with different pathogenesis and host preference can provide a means of identifying the specifics of adaptive differences. Streptococcus dysgalactiae (SD) is comprised of two subspecies: S. dysgalactiae subsp. equisimilis is both a human commensal organism and a human pathogen, and S. dysgalactiae subsp. dysgalactiae is strictly an animal pathogen. Here, we present complete genome sequences for both taxa, with analyses involving other species of Streptococcus but focusing on adaptation in the SD species group. We found little evidence for enrichment in biochemical categories of genes carried by each SD strain, however, differences in the virulence gene repertoire were apparent. Some of the differences could be ascribed to prophage and integrative conjugative elements. We identified approximately 9% of the nonrecombinant core genome to be under positive selection, some of which involved known virulence factors in other bacteria. Analyses of proteomes by pooling data across genes, by biochemical category, clade, or branch, provided evidence for increased rates of evolution in several gene categories, as well as external branches of the tree. Promoters were primarily evolving under purifying selection but with certain categories of genes evolving faster. Many of these fast-evolving categories were the same as those associated with rapid evolution in proteins. Overall, these results suggest that adaptation to changing environments and new hosts in the SD species group has involved the acquisition of key virulence genes along with selection of orthologous protein-coding loci and operon promoters.