The pan-genome of Lactobacillus reuteri strains originating from the pig gastrointestinal tract.

BACKGROUND: Lactobacillus reuteri is a gut symbiont of a wide variety of vertebrate species that has diversified into distinct phylogenetic clades which are to a large degree host-specific. Previous work demonstrated host specificity in mice and begun to determine the mechanisms by which gut colonisation and host restriction is achieved. However, how L. reuteri strains colonise the gastrointestinal (GI) tract of pigs is unknown. RESULTS: To gain insight into the ecology of L. reuteri in the pig gut, the genome sequence of the porcine small intestinal isolate L. reuteri ATCC 53608 was completed and consisted of a chromosome of 1.94 Mbp and two plasmids of 138.5 kbp and 9.09 kbp, respectively. Furthermore, we generated draft genomes of four additional L. reuteri strains isolated from pig faeces or lower GI tract, lp167-67, pg-3b, 20-2 and 3c6, and subjected all five genomes to a comparative genomic analysis together with the previously completed genome of strain I5007. A phylogenetic analysis based on whole genomes showed that porcine L. reuteri strains fall into two distinct clades, as previously suggested by multi-locus sequence analysis. These six pig L. reuteri genomes contained a core set of 1364 orthologous gene clusters, as determined by OrthoMCL analysis, that contributed to a pan-genome totalling 3373 gene clusters. Genome comparisons of the six pig L. reuteri strains with 14 L. reuteri strains from other host origins gave a total pan-genome of 5225 gene clusters that included a core genome of 851 gene clusters but revealed that there were no pig-specific genes per se. However, genes specific for and conserved among strains of the two pig phylogenetic lineages were detected, some of which encoded cell surface proteins that could contribute to the diversification of the two lineages and their observed host specificity. CONCLUSIONS: This study extends the phylogenetic analysis of L. reuteri strains at a genome-wide level, pointing to distinct evolutionary trajectories of porcine L. reuteri lineages, and providing new insights into the genomic events in L. reuteri that occurred during specialisation to their hosts. The occurrence of two distinct pig-derived clades may reflect differences in host genotype, environmental factors such as dietary components or to evolution from ancestral strains of human and rodent origin following contact with pig populations.

Genomic and physiological variability within Group II (non-proteolytic) Clostridium botulinum.

BACKGROUND: Clostridium botulinum is a group of four physiologically and phylogenetically distinct bacteria that produce botulinum neurotoxin. While studies have characterised variability between strains of Group I (proteolytic) C. botulinum, the genetic and physiological variability and relationships between strains within Group II (non-proteolytic) C. botulinum are not well understood. In this study the genome of Group II strain C. botulinum Eklund 17B (NRP) was sequenced and used to construct a whole genome DNA microarray. This was used in a comparative genomic indexing study to compare the relatedness of 43 strains of Group II C. botulinum (14 type B, 24 type E and 5 type F). These results were compared with characteristics determined from physiological tests. RESULTS: Whole genome indexing showed that strains of Group II C. botulinum isolated from a wide variety of environments over more than 75 years clustered together indicating the genetic background of Group II C. botulinum is stable. Further analysis showed that strains forming type B or type F toxin are closely related with only toxin cluster genes targets being unique to either type. Strains producing type E toxin formed a separate subset. Carbohydrate fermentation tests supported the observation that type B and F strains form a separate subset to type E strains. All the type F strains and most of type B strains produced acid from amylopectin, amylose and glycogen whereas type E strains did not. However, these two subsets did not differ strongly in minimum growth temperature or maximum NaCl concentration for growth. No relationship was found between tellurite resistance and toxin type despite all the tested type B and type F strains carrying tehB, while the sequence was absent or diverged in all type E strains. CONCLUSIONS: Although Group II C. botulinum form a tight genetic group, genomic and physiological analysis indicates there are two distinct subsets within this group. All type B strains and type F strains are in one subset and all type E strains in the other.

Complete genome sequence of the proteolytic Clostridium botulinum type A5 (B3') strain H04402 065.

H04402 065 is one of a very small group of strains of proteolytic Clostridium botulinum that form type A5 neurotoxin. Here, we report the complete 3.9-Mb genome sequence and annotation of strain H04402 065, which was isolated from a botulism patient in the United Kingdom in 2004.

Laboratory Stock Variants of the Archetype Silver Resistance Plasmid pMG101 Demonstrate Plasmid Fusion, Loss of Transmissibility, and Transposition of Tn7/pco/sil Into the Host Chromosome.

Salmonella Typhimurium carrying the multidrug resistance (MDR) plasmid pMG101 was isolated from three burns patients in Boston United States in 1973. pMG101 was transferrable into other Salmonella spp. and Escherichia coli hosts and carried what was a novel and unusual combination of AMR genes and silver resistance. Previously published short-read DNA sequence of pMG101 showed that it was a 183.5Kb IncHI plasmid, where a Tn7-mediated transposition of pco/sil resistance genes into the chromosome of the E. coli K-12 J53 host strain had occurred. We noticed differences in streptomycin resistance and plasmid size between two stocks of E. coli K-12 J53 pMG101 we possessed, which had been obtained from two different laboratories (pMG101-A and pMG101-B). Long-read sequencing (PacBio) of the two strains unexpectedly revealed plasmid and chromosomal rearrangements in both. pMG101-A is a non-transmissible 383Kb closed-circular plasmid consisting of an IncHI2 plasmid sequence fused to an IncFI/FIIA plasmid. pMG101-B is a mobile closed-circular 154 Kb IncFI/FIIA plasmid. Sequence identity of pMG101-B with the fused IncFI/IncFIIA region of pMG101-A was >99%. Assembled host sequence reads of pMG101-B showed Tn7-mediated transposition of pco/sil into the E. coli J53 chromosome between yhiM and yhiN. Long read sequence data in combination with laboratory experiments have demonstrated large scale changes in pMG101. Loss of conjugation function and movement of resistance genes into the chromosome suggest that even under long-term laboratory storage, mobile genetic elements such as transposons and insertion sequences can drive the evolution of plasmids and host. This study emphasises the importance of utilising long read sequencing technologies of plasmids and host strains at the earliest opportunity.

A commensal gone bad: complete genome sequence of the prototypical enterotoxigenic Escherichia coli strain H10407.

In most cases, Escherichia coli exists as a harmless commensal organism, but it may on occasion cause intestinal and/or extraintestinal disease. Enterotoxigenic E. coli (ETEC) is the predominant cause of E. coli-mediated diarrhea in the developing world and is responsible for a significant portion of pediatric deaths. In this study, we determined the complete genomic sequence of E. coli H10407, a prototypical strain of enterotoxigenic E. coli, which reproducibly elicits diarrhea in human volunteer studies. We performed genomic and phylogenetic comparisons with other E. coli strains, revealing that the chromosome is closely related to that of the nonpathogenic commensal strain E. coli HS and to those of the laboratory strains E. coli K-12 and C. Furthermore, these analyses demonstrated that there were no chromosomally encoded factors unique to any sequenced ETEC strains. Comparison of the E. coli H10407 plasmids with those from several ETEC strains revealed that the plasmids had a mosaic structure but that several loci were conserved among ETEC strains. This study provides a genetic context for the vast amount of experimental and epidemiological data that have been published.

Identification and characterisation of enteroaggregative Escherichia coli subtypes associated with human disease.

Enteroaggregative E. coli (EAEC) are a major cause of diarrhoea worldwide. Due to their heterogeneity and carriage in healthy individuals, identification of diagnostic virulence markers for pathogenic strains has been difficult. In this study, we have determined phenotypic and genotypic differences between EAEC strains of sequence types (STs) epidemiologically associated with asymptomatic carriage (ST31) and diarrhoeal disease (ST40). ST40 strains demonstrated significantly enhanced intestinal adherence, biofilm formation, and pro-inflammatory interleukin-8 secretion compared with ST31 isolates. This was independent of whether strains were derived from diarrhoea patients or healthy controls. Whole genome sequencing revealed differences in putative virulence genes encoding aggregative adherence fimbriae, E. coli common pilus, flagellin and EAEC heat-stable enterotoxin 1. Our results indicate that ST40 strains have a higher intrinsic potential of human pathogenesis due to a specific combination of virulence-related factors which promote host cell colonization and inflammation. These findings may contribute to the development of genotypic and/or phenotypic markers for EAEC strains of high virulence.

The genome of Burkholderia cenocepacia J2315, an epidemic pathogen of cystic fibrosis patients.

Bacterial infections of the lungs of cystic fibrosis (CF) patients cause major complications in the treatment of this common genetic disease. Burkholderia cenocepacia infection is particularly problematic since this organism has high levels of antibiotic resistance, making it difficult to eradicate; the resulting chronic infections are associated with severe declines in lung function and increased mortality rates. B. cenocepacia strain J2315 was isolated from a CF patient and is a member of the epidemic ET12 lineage that originated in Canada or the United Kingdom and spread to Europe. The 8.06-Mb genome of this highly transmissible pathogen comprises three circular chromosomes and a plasmid and encodes a broad array of functions typical of this metabolically versatile genus, as well as numerous virulence and drug resistance functions. Although B. cenocepacia strains can be isolated from soil and can be pathogenic to both plants and man, J2315 is representative of a lineage of B. cenocepacia rarely isolated from the environment and which spreads between CF patients. Comparative analysis revealed that ca. 21% of the genome is unique in comparison to other strains of B. cenocepacia, highlighting the genomic plasticity of this species. Pseudogenes in virulence determinants suggest that the pathogenic response of J2315 may have been recently selected to promote persistence in the CF lung. The J2315 genome contains evidence that its unique and highly adapted genetic content has played a significant role in its success as an epidemic CF pathogen.

The mucin-degradation strategy of Ruminococcus gnavus: The importance of intramolecular trans-sialidases.

We previously identified and characterized an intramolecular trans-sialidase (IT-sialidase) in the gut symbiont Ruminococcus gnavus ATCC 29149, which is associated to the ability of the strain to grow on mucins. In this work we have obtained and analyzed the draft genome sequence of another R. gnavus mucin-degrader, ATCC 35913, isolated from a healthy individual. Transcriptomics analyses of both ATCC 29149 and ATCC 35913 strains confirmed that the strategy utilized by R. gnavus for mucin-degradation is focused on the utilization of terminal mucin glycans. R. gnavus ATCC 35913 also encodes a predicted IT-sialidase and harbors a Nan cluster dedicated to sialic acid utilization. We showed that the Nan cluster was upregulated when the strains were grown in presence of mucin. In addition we demonstrated that both R. gnavus strains were able to grow on 2,7-anyhydro-Neu5Ac, the IT-sialidase transglycosylation product, as a sole carbon source. Taken together these data further support the hypothesis that IT-sialidase expressing gut microbes, provide commensal bacteria such as R. gnavus with a nutritional competitive advantage, by accessing and transforming a source of nutrient to their own benefit.

Bacterial antimicrobial metal ion resistance.

Metals such as mercury, arsenic, copper and silver have been used in various forms as antimicrobials for thousands of years with until recently, little understanding of their mode of action. The discovery of antibiotics and new organic antimicrobial compounds during the twentieth century saw a general decline in the clinical use of antimicrobial metal compounds, with the exception of the rediscovery of the use of silver for burns treatments and niche uses for other metal compounds. Antibiotics and new antimicrobials were regarded as being safer for the patient and more effective than the metal-based compounds they supplanted. Bacterial metal ion resistances were first discovered in the second half of the twentieth century. The detailed mechanisms of resistance have now been characterized in a wide range of bacteria. As the use of antimicrobial metals is limited, it is legitimate to ask: are antimicrobial metal resistances in pathogenic and commensal bacteria important now? This review details the new, rediscovered and 'never went away' uses of antimicrobial metals; examines the prevalence and linkage of antimicrobial metal resistance genes to other antimicrobial resistance genes; and examines the evidence for horizontal transfer of these genes between bacteria. Finally, we discuss the possible implications of the widespread dissemination of these resistances on re-emergent uses of antimicrobial metals and how this could impact upon the antibiotic resistance problem.

Draft Genome Sequence of a Novel Lactobacillus salivarius Strain Isolated from Piglet.

Lactobacillus salivarius is part of the vertebrate indigenous microbiota of the gastrointestinal tract, oral cavity, and milk. The properties associated with some L. salivarius strains have led to their use as probiotics. Here we describe the draft genome of the pig isolate L. salivarius cp400, providing insights into host-niche specialization.

Complete Genome Sequence of the Campylobacter coli Clinical Isolate 15-537360.

Campylobacter coli strain 15-537360 was originally isolated in 2001 from a 42-year-old patient with gastroenteritis. Here, we report its complete genome sequence, which comprises a 1.7-Mbp chromosome and a 29-kbp conjugative cryptic plasmid. This is the first complete genome sequence of a clinical isolate of C. coli.

A small predatory core genome in the divergent marine Bacteriovorax marinus SJ and the terrestrial Bdellovibrio bacteriovorus.

Bacteriovorax marinus SJ is a predatory delta-proteobacterium isolated from a marine environment. The genome sequence of this strain provides an interesting contrast to that of the terrestrial predatory bacterium Bdellovibrio bacteriovorus HD100. Based on their predatory lifestyle, Bacteriovorax were originally designated as members of the genus Bdellovibrio but subsequently were re-assigned to a new genus and family based on genetic and phenotypic differences. B. marinus attaches to gram-negative bacteria, penetrates through the cell wall to form a bdelloplast, in which it replicates, as shown using microscopy. Bacteriovorax is distinct, as it shares only 30% of its gene products with its closest sequenced relatives. Remarkably, 34% of predicted genes over 500 nt in length were completely unique with no significant matches in the databases. As expected, Bacteriovorax shares several characteristic loci with the other delta-proteobacteria. A geneset shared between Bacteriovorax and Bdellovibrio that is not conserved among other delta-proteobacteria such as Myxobacteria (which destroy prey bacteria externally via lysis), or the non-predatory Desulfo-bacteria and Geobacter species was identified. These 291 gene orthologues common to both Bacteriovorax and Bdellovibrio may be the key indicators of host-interaction predatory-specific processes required for prey entry. The locus from Bdellovibrio bacteriovorus is implicated in the switch from predatory to prey/host-independent growth. Although the locus is conserved in B. marinus, the sequence has only limited similarity. The results of this study advance understanding of both the similarities and differences between Bdellovibrio and Bacteriovorax and confirm the distant relationship between the two and their separation into different families.

The complete genome, comparative and functional analysis of Stenotrophomonas maltophilia reveals an organism heavily shielded by drug resistance determinants.

BACKGROUND: Stenotrophomonas maltophilia is a nosocomial opportunistic pathogen of the Xanthomonadaceae. The organism has been isolated from both clinical and soil environments in addition to the sputum of cystic fibrosis patients and the immunocompromised. Whilst relatively distant phylogenetically, the closest sequenced relatives of S. maltophilia are the plant pathogenic xanthomonads. RESULTS: The genome of the bacteremia-associated isolate S. maltophilia K279a is 4,851,126 bp and of high G+C content. The sequence reveals an organism with a remarkable capacity for drug and heavy metal resistance. In addition to a number of genes conferring resistance to antimicrobial drugs of different classes via alternative mechanisms, nine resistance-nodulation-division (RND)-type putative antimicrobial efflux systems are present. Functional genomic analysis confirms a role in drug resistance for several of the novel RND efflux pumps. S. maltophilia possesses potentially mobile regions of DNA and encodes a number of pili and fimbriae likely to be involved in adhesion and biofilm formation that may also contribute to increased antimicrobial drug resistance. CONCLUSION: The panoply of antimicrobial drug resistance genes and mobile genetic elements found suggests that the organism can act as a reservoir of antimicrobial drug resistance determinants in a clinical environment, which is an issue of considerable concern.

A common genomic framework for a diverse assembly of plasmids in the symbiotic nitrogen fixing bacteria.

This work centres on the genomic comparisons of two closely-related nitrogen-fixing symbiotic bacteria, Rhizobium leguminosarum biovar viciae 3841 and Rhizobium etli CFN42. These strains maintain a stable genomic core that is also common to other rhizobia species plus a very variable and significant accessory component. The chromosomes are highly syntenic, whereas plasmids are related by fewer syntenic blocks and have mosaic structures. The pairs of plasmids p42f-pRL12, p42e-pRL11 and p42b-pRL9 as well large parts of p42c with pRL10 are shown to be similar, whereas the symbiotic plasmids (p42d and pRL10) are structurally unrelated and seem to follow distinct evolutionary paths. Even though purifying selection is acting on the whole genome, the accessory component is evolving more rapidly. This component is constituted largely for proteins for transport of diverse metabolites and elements of external origin. The present analysis allows us to conclude that a heterogeneous and quickly diversifying group of plasmids co-exists in a common genomic framework.

Sequence-based analysis of pQBR103; a representative of a unique, transfer-proficient mega plasmid resident in the microbial community of sugar beet.

The plasmid pQBR103 was found within Pseudomonas populations colonizing the leaf and root surfaces of sugar beet plants growing at Wytham, Oxfordshire, UK. At 425 kb it is the largest self-transmissible plasmid yet sequenced from the phytosphere. It is known to enhance the competitive fitness of its host, and parts of the plasmid are known to be actively transcribed in the plant environment. Analysis of the complete sequence of this plasmid predicts a coding sequence (CDS)-rich genome containing 478 CDSs and an exceptional degree of genetic novelty; 80% of predicted coding sequences cannot be ascribed a function and 60% are orphans. Of those to which function could be assigned, 40% bore greatest similarity to sequences from Pseudomonas spp, and the majority of the remainder showed similarity to other gamma-proteobacterial genera and plasmids. pQBR103 has identifiable regions presumed responsible for replication and partitioning, but despite being tra+ lacks the full complement of any previously described conjugal transfer functions. The DNA sequence provided few insights into the functional significance of plant-induced transcriptional regions, but suggests that 14% of CDSs may be expressed (11 CDSs with functional annotation and 54 without), further highlighting the ecological importance of these novel CDSs. Comparative analysis indicates that pQBR103 shares significant regions of sequence with other plasmids isolated from sugar beet plants grown at the same geographic location. These plasmid sequences indicate there is more novelty in the mobile DNA pool accessible to phytosphere pseudomonas than is currently appreciated or understood.

Cell-cell signaling in Xanthomonas campestris involves an HD-GYP domain protein that functions in cyclic di-GMP turnover.

HD-GYP is a protein domain of unknown biochemical function implicated in bacterial signaling and regulation. In the plant pathogen Xanthomonas campestris pv. campestris, the synthesis of virulence factors and dispersal of biofilms are positively controlled by a two-component signal transduction system comprising the HD-GYP domain regulatory protein RpfG and cognate sensor RpfC and by cell-cell signaling mediated by the diffusible signal molecule DSF (diffusible signal factor). The RpfG/RpfC two-component system has been implicated in DSF perception and signal transduction. Here we show that the role of RpfG is to degrade the unusual nucleotide cyclic di-GMP, an activity associated with the HD-GYP domain. Mutation of the conserved H and D residues of the isolated HD-GYP domain resulted in loss of both the enzymatic activity against cyclic di-GMP and the regulatory activity in virulence factor synthesis. Two other protein domains, GGDEF and EAL, are already implicated in the synthesis and degradation respectively of cyclic di-GMP. As with GGDEF and EAL domains, the HD-GYP domain is widely distributed in free-living bacteria and occurs in plant and animal pathogens, as well as beneficial symbionts and organisms associated with a range of environmental niches. Identification of the role of the HD-GYP domain thus increases our understanding of a signaling network whose importance to the lifestyle of diverse bacteria is now emerging.

The genome of Rhizobium leguminosarum has recognizable core and accessory components.

BACKGROUND: Rhizobium leguminosarum is an alpha-proteobacterial N2-fixing symbiont of legumes that has been the subject of more than a thousand publications. Genes for the symbiotic interaction with plants are well studied, but the adaptations that allow survival and growth in the soil environment are poorly understood. We have sequenced the genome of R. leguminosarum biovar viciae strain 3841. RESULTS: The 7.75 Mb genome comprises a circular chromosome and six circular plasmids, with 61% G+C overall. All three rRNA operons and 52 tRNA genes are on the chromosome; essential protein-encoding genes are largely chromosomal, but most functional classes occur on plasmids as well. Of the 7,263 protein-encoding genes, 2,056 had orthologs in each of three related genomes (Agrobacterium tumefaciens, Sinorhizobium meliloti, and Mesorhizobium loti), and these genes were over-represented in the chromosome and had above average G+C. Most supported the rRNA-based phylogeny, confirming A. tumefaciens to be the closest among these relatives, but 347 genes were incompatible with this phylogeny; these were scattered throughout the genome but were over-represented on the plasmids. An unexpectedly large number of genes were shared by all three rhizobia but were missing from A. tumefaciens. CONCLUSION: Overall, the genome can be considered to have two main components: a 'core', which is higher in G+C, is mostly chromosomal, is shared with related organisms, and has a consistent phylogeny; and an 'accessory' component, which is sporadic in distribution, lower in G+C, and located on the plasmids and chromosomal islands. The accessory genome has a different nucleotide composition from the core despite a long history of coexistence.

Extensive DNA inversions in the B. fragilis genome control variable gene expression.

The obligately anaerobic bacterium Bacteroides fragilis, an opportunistic pathogen and inhabitant of the normal human colonic microbiota, exhibits considerable within-strain phase and antigenic variation of surface components. The complete genome sequence has revealed an unusual breadth (in number and in effect) of DNA inversion events that potentially control expression of many different components, including surface and secreted components, regulatory molecules, and restriction-modification proteins. Invertible promoters of two different types (12 group 1 and 11 group 2) were identified. One group has inversion crossover (fix) sites similar to the hix sites of Salmonella typhimurium. There are also four independent intergenic shufflons that potentially alter the expression and function of varied genes. The composition of the 10 different polysaccharide biosynthesis gene clusters identified (7 with associated invertible promoters) suggests a mechanism of synthesis similar to the O-antigen capsules of Escherichia coli.

Heterologous expression of heterotrophic nitrification genes.

Paracoccus denitrificans is a heterotrophic organism capable of oxidizing ammonia to nitrite during growth on an organic carbon and energy source. This pathway, termed heterotrophic nitrification, requires the concerted action of an ammonia monooxygenase (AMO) and hydroxylamine oxidase (HAO). The genes required for heterotrophic nitrification have been isolated by introducing a Pa. denitrificans genomic library into Pseudomonas putida and screening for the accumulation of nitrite. In contrast to the situation in chemolithoautotrophic ammonia oxidizers, the genes encoding AMO and HAO are present in single linked copies in the genome of Pa. denitrificans. AMO from Pa. denitrificans expressed in Ps. putida is capable of oxidizing ethene (ethylene) to epoxyethane (ethylene oxide), which is indicative of a relaxed substrate specificity. Further, when expressed in the methylotroph Methylobacterium extorquens AM1, the AMO endows on this organism the ability to grow on ethene and methane. Thus, the Pa. denitrificans AMO is capable of oxidizing methane to methanol, as is the case for the AMO from Nitrosomonas europaea. The heterotrophic nitrification genes are moderately toxic in M. extorquens, more toxic in Ps. putida, and non-toxic in Escherichia coli. Toxicity is due to the activity of the gene products in M. extorquens, and both expression and activity in Ps. putida. This is the first time that the genes encoding an active AMO have been expressed in a heterologous host.