Andrew K. Benson

Comparative genomics of the lactic acid bacteria

Lactic acid-producing bacteria are associated with various plant and animal niches and play a key role in the production of fermented foods and beverages. We report nine genome sequences representing the phylogenetic and functional diversity of these bacteria. The small genomes of lactic acid bacteria encode a broad repertoire of transporters for efficient carbon and nitrogen acquisition from the nutritionally rich environments they inhabit and reflect a limited range of biosynthetic capabilities that indicate both prototrophic and auxotrophic strains. Phylogenetic analyses, comparison of gene content across the group, and reconstruction of ancestral gene sets indicate a combination of extensive gene loss and key gene acquisitions via horizontal gene transfer during the coevolution of lactic acid bacteria with their habitats.

Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors

In vertebrates, including humans, individuals harbor gut microbial communities whose species composition and relative proportions of dominant microbial groups are tremendously varied. Although external and stochastic factors clearly contribute to the individuality of the microbiota, the fundamental principles dictating how environmental factors and host genetic factors combine to shape this complex ecosystem are largely unknown and require systematic study. Here we examined factors that affect microbiota composition in a large (n = 645) mouse advanced intercross line originating from a cross between C57BL/6J and an ICR-derived outbred line (HR). Quantitative pyrosequencing of the microbiota defined a core measurable microbiota (CMM) of 64 conserved taxonomic groups that varied quantitatively across most animals in the population. Although some of this variation can be explained by litter and cohort effects, individual host genotype had a measurable contribution. Testing of the CMM abundances for cosegregation with 530 fully informative SNP markers identified 18 host quantitative trait loci (QTL) that show significant or suggestive genome-wide linkage with relative abundances of specific microbial taxa. These QTL affect microbiota composition in three ways; some loci control individual microbial species, some control groups of related taxa, and some have putative pleiotropic effects on groups of distantly related organisms. These data provide clear evidence for the importance of host genetic control in shaping individual microbiome diversity in mammals, a key step toward understanding the factors that govern the assemblages of gut microbiota associated with complex diseases.

The Evolution of Host Specialization in the Vertebrate Gut Symbiont Lactobacillus reuteri

Recent research has provided mechanistic insight into the important contributions of the gut microbiota to vertebrate biology, but questions remain about the evolutionary processes that have shaped this symbiosis. In the present study, we showed in experiments with gnotobiotic mice that the evolution of Lactobacillus reuteri with rodents resulted in the emergence of host specialization. To identify genomic events marking adaptations to the murine host, we compared the genome of the rodent isolate L. reuteri 100-23 with that of the human isolate L. reuteri F275, and we identified hundreds of genes that were specific to each strain. In order to differentiate true host-specific genome content from strain-level differences, comparative genome hybridizations were performed to query 57 L. reuteri strains originating from six different vertebrate hosts in combination with genome sequence comparisons of nine strains encompassing five phylogenetic lineages of the species. This approach revealed that rodent strains, although showing a high degree of genomic plasticity, possessed a specific genome inventory that was rare or absent in strains from other vertebrate hosts. The distinct genome content of L. reuteri lineages reflected the niche characteristics in the gastrointestinal tracts of their respective hosts, and inactivation of seven out of eight representative rodent-specific genes in L. reuteri 100-23 resulted in impaired ecological performance in the gut of mice. The comparative genomic analyses suggested fundamentally different trends of genome evolution in rodent and human L. reuteri populations, with the former possessing a large and adaptable pan-genome while the latter being subjected to a process of reductive evolution. In conclusion, this study provided experimental evidence and a molecular basis for the evolution of host specificity in a vertebrate gut symbiont, and it identified genomic events that have shaped this process.

Diversification of the gut symbiont Lactobacillus reuteri as a result of host-driven evolution

The vertebrate digestive tract, including that of humans, is the habitat to trillions of bacteria that are of significant importance to host biology and health. Although these communities are often postulated to have coevolved with their hosts, evidence is lacking, yet critical for our understanding of microbial symbiosis in vertebrates. To gain insight into the evolution of a gut symbiont, we have characterized the population genetic structure and phylogeny of Lactobacillus reuteri strains isolated from six different host species (human, mouse, rat, pig, chicken and turkey) using Amplified-Fragment Length Polymorphism (AFLP) and Multi-Locus Sequence Analysis (MLSA). The results revealed considerable genetic heterogeneity within the L. reuteri population and distinct monophyletic clades reflecting host origin but not provenance. The evolutionary patterns detected indicate a long-term association of L. reuteri lineages with particular vertebrate species and host-driven diversification. Results from a competition experiment in a gnotobiotic mouse model revealed that rodent isolates showed elevated ecological performance, indicating that evolution of L. reuteri lineages was adaptive. These findings provide evidence that some vertebrate gut microbes are not promiscuous, but have diversified into host-adapted lineages by a long-term evolutionary process, allowing the development of a highly specialized symbiosis.

Genome evolution in major Escherichia coli O157:H7 lineages

BackgroundGenetic analysis of Escherichia coli O157:H7 strains has shown divergence into two distinct lineages, lineages I and II, that appear to have distinct ecological characteristics, with lineage I strains more commonly associated with human disease. In this study, microarray-based comparative genomic hybridization (CGH) was used to identify genomic differences among 31 E. coli O157:H7 strains that belong to various phage types (PTs) and different lineage-specific polymorphism assay (LSPA) types.ResultsA total of 4,084 out of 6,057 ORFs were detected in all E. coli O157:H7 strains and 1,751 were variably present or absent. Based on this data, E. coli O157:H7 strains were divided into three distinct clusters, which consisted of 15 lineage I (LSPA type 111111), four lineage I/II (designated in this study) (LSPA type 211111) and 12 lineage II strains (LSPA 222222, 222211, 222212, and 222221), respectively. Eleven different genomic regions that were dominant in lineage I strains (present in ≥80% of lineage I and absent from ≥ 92% of lineage II strains) spanned segments containing as few as two and up to 25 ORFs each. These regions were identified within E. coli Sakai S-loops # 14, 16, 69, 72, 78, 83, 85, 153 and 286, Sakai phage 10 (S-loops # 91, 92 and 93) and a genomic backbone region. All four lineage I/II strains were of PT 2 and possessed eight of these 11 lineage I-dominant loci. Several differences in virulence-associated loci were noted between lineage I and lineage II strains, including divergence within S-loop 69, which encodes Shiga toxin 2, and absence of the non-LEE encoded effector genes nleF and nleH1-2 and the perC homologue gene pchD in lineage II strains.ConclusionCGH data suggest the existence of two dominant lineages as well as LSPA type and PT-related subgroups within E. coli O157:H7. The genomic composition of these subgroups supports the phylogeny that has been inferred from other methods and further suggests that genomic divergence from an ancestral form and lateral gene transfer have contributed to their evolution. The genomic features identified in this study may contribute to apparent differences in the epidemiology and ecology of strains of different E. coli O157:H7 lineages.

Role of sigma(B) in adaptation of Listeria monocytogenes to growth at low temperature.

The activity of sigma(B) in Listeria monocytogenes is stimulated by high osmolarity and is necessary for efficient uptake of osmoprotectants. Here we demonstrate that, during cold shock, sigma(B) contributes to adaptation in a growth phase-dependent manner and is necessary for efficient accumulation of betaine and carnitine as cryoprotectants.

Next-Generation Sequencing Reveals Significant Bacterial Diversity of Botrytized Wine

While wine fermentation has long been known to involve complex microbial communities, the composition and role of bacteria other than a select set of lactic acid bacteria (LAB) has often been assumed either negligible or detrimental. This study served as a pilot study for using barcoded amplicon next-generation sequencing to profile bacterial community structure in wines and grape musts, comparing the taxonomic depth achieved by sequencing two different domains of prokaryotic 16S rDNA (V4 and V5). This study was designed to serve two goals: 1) to empirically determine the most taxonomically informative 16S rDNA target region for barcoded amplicon sequencing of wine, comparing V4 and V5 domains of bacterial 16S rDNA to terminal restriction fragment length polymorphism (TRFLP) of LAB communities; and 2) to explore the bacterial communities of wine fermentation to better understand the biodiversity of wine at a depth previously unattainable using other techniques. Analysis of amplicons from the V4 and V5 provided similar views of the bacterial communities of botrytized wine fermentations, revealing a broad diversity of low-abundance taxa not traditionally associated with wine, as well as atypical LAB communities initially detected by TRFLP. The V4 domain was determined as the more suitable read for wine ecology studies, as it provided greater taxonomic depth for profiling LAB communities. In addition, targeted enrichment was used to isolate two species of Alphaproteobacteria from a finished fermentation. Significant differences in diversity between inoculated and uninoculated samples suggest that Saccharomyces inoculation exerts selective pressure on bacterial diversity in these fermentations, most notably suppressing abundance of acetic acid bacteria. These results determine the bacterial diversity of botrytized wines to be far higher than previously realized, providing further insight into the fermentation dynamics of these wines, and demonstrate the utility of next-generation sequencing for wine ecology studies.

Murine Gut Microbiota Is Defined by Host Genetics and Modulates Variation of Metabolic Traits

The gastrointestinal tract harbors a complex and diverse microbiota that has an important role in host metabolism. Microbial diversity is influenced by a combination of environmental and host genetic factors and is associated with several polygenic diseases. In this study we combined next-generation sequencing, genetic mapping, and a set of physiological traits of the BXD mouse population to explore genetic factors that explain differences in gut microbiota and its impact on metabolic traits. Molecular profiling of the gut microbiota revealed important quantitative differences in microbial composition among BXD strains. These differences in gut microbial composition are influenced by host-genetics, which is complex and involves many loci. Linkage analysis defined Quantitative Trait Loci (QTLs) restricted to a particular taxon, branch or that influenced the variation of taxa across phyla. Gene expression within the gastrointestinal tract and sequence analysis of the parental genomes in the QTL regions uncovered candidate genes with potential to alter gut immunological profiles and impact the balance between gut microbial communities. A QTL region on Chr 4 that overlaps several interferon genes modulates the population of Bacteroides, and potentially Bacteroidetes and Firmicutes–the predominant BXD gut phyla. Irak4, a signaling molecule in the Toll-like receptor pathways is a candidate for the QTL on Chr15 that modulates Rikenellaceae, whereas Tgfb3, a cytokine modulating the barrier function of the intestine and tolerance to commensal bacteria, overlaps a QTL on Chr 12 that influence Prevotellaceae. Relationships between gut microflora, morphological and metabolic traits were uncovered, some potentially a result of common genetic sources of variation.

Identification of Common Subpopulations of Non-Sorbitol-Fermenting, β-Glucuronidase-Negative Escherichia coli O157:H7 from Bovine Production Environments and Human Clinical Samples

ABSTRACT Non-sorbitol-fermenting, β-glucuronidase-negative Escherichia coli O157:H7 strains are regarded as a clone complex, and populations from different geographical locations are believed to share a recent common ancestor. Despite their relatedness, high-resolution genotyping methods can detect significant genome variation among different populations. Phylogenetic analysis of high-resolution genotyping data from these strains has shown that subpopulations from geographically unlinked continents can be divided into two primary phylogenetic lineages, termed lineage I and lineage II, and limited studies of the distribution of these lineages suggest there could be differences in their propensity to cause disease in humans or to be transmitted to humans. Because the genotyping methods necessary to discriminate the two lineages are tedious and subjective, these methods are not particularly suited for studying the large sets of strains that are required to systematically evaluate the ecology and transmission characteristics of these lineages. To overcome this limitation, we have developed a lineage-specific polymorphism assay (LSPA) that can readily distinguish between the lineage I and lineage II subpopulations. In the studies reported here, we describe the development of a six-marker test (LSPA-6) and its validation in a side-by-side comparison with octamer-based genome scanning. Analysis of over 1,400 O157:H7 strains with the LSPA-6 demonstrated that five genotypes comprise over 91% of the strains, suggesting that these subpopulations may be widespread.

Genome Diversification in Phylogenetic Lineages I and II of Listeria monocytogenes: Identification of Segments Unique to Lineage II Populations

Thirteen different serotypes of Listeria monocytogenes can be distinguished on the basis of variation in somatic and flagellar antigens. Although the known virulence genes are present in all serotypes, greater than 90% of human cases of listeriosis are caused by serotypes 1/2a, 1/2b, and 4b and nearly all outbreaks of food-borne listeriosis have been caused by serotype 4b strains. Phylogenetic analysis of these three common clinical serotypes places them into two different lineages, with serotypes 1/2b and 4b belonging to lineage I and 1/2a belonging to lineage II. To begin examining evolution of the genome in these serotypes, DNA microarray analysis was used to identify lineage-specific and serotype-specific differences in genome content. A set of 44 strains representing serotypes 1/2a, 1/2b, and 4b was probed with a shotgun DNA microarray constructed from the serotype 1/2a strain 10403s. Clones spanning 47 different genes in 16 different contiguous segments relative to the lineage II 1/2a genome were found to be absent in all lineage I strains tested (serotype 4b and 1/2b) and an additional nine were altered exclusively in 4b strains. Southern hybridization confirmed that conserved alterations were, in all but two loci, due to absence of the segments from the genome. Genes within these contiguous segments comprise five functional categories, including genes involved in synthesis of cell surface molecules and regulation of virulence gene expression. Phylogenetic reconstruction and examination of compositional bias in the regions of difference are consistent with a model in which the ancestor of the two lineages had the 1/2 somatic serotype and the regions absent in the lineage I genome arose by loss of ancestral sequences.