Jaehyoung Kim

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.

Resistant Starches Types 2 and 4 Have Differential Effects on the Composition of the Fecal Microbiota in Human Subjects

Background To systematically develop dietary strategies based on resistant starch (RS) that modulate the human gut microbiome, detailed in vivo studies that evaluate the effects of different forms of RS on the community structure and population dynamics of the gut microbiota are necessary. The aim of the present study was to gain a community wide perspective of the effects of RS types 2 (RS2) and 4 (RS4) on the fecal microbiota in human individuals. Methods and Findings Ten human subjects consumed crackers for three weeks each containing either RS2, RS4, or native starch in a double-blind, crossover design. Multiplex sequencing of 16S rRNA tags revealed that both types of RS induced several significant compositional alterations in the fecal microbial populations, with differential effects on community structure. RS4 but not RS2 induced phylum-level changes, significantly increasing Actinobacteria and Bacteroidetes while decreasing Firmicutes. At the species level, the changes evoked by RS4 were increases in Bifidobacterium adolescentis and Parabacteroides distasonis, while RS2 significantly raised the proportions of Ruminococcus bromii and Eubacterium rectale when compared to RS4. The population shifts caused by RS4 were numerically substantial for several taxa, leading for example, to a ten-fold increase in bifidobacteria in three of the subjects, enriching them to 18–30% of the fecal microbial community. The responses to RS and their magnitudes varied between individuals, and they were reversible and tightly associated with the consumption of RS. Conclusion Our results demonstrate that RS2 and RS4 show functional differences in their effect on human fecal microbiota composition, indicating that the chemical structure of RS determines its accessibility by groups of colonic bacteria. The findings imply that specific bacterial populations could be selectively targeted by well designed functional carbohydrates, but the inter-subject variations in the response to RS indicates that such strategies might benefit from more personalized approaches.

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.

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.

Ancestral Divergence, Genome Diversification, and Phylogeographic Variation in Subpopulations of Sorbitol-Negative, β-Glucuronidase-Negative Enterohemorrhagic Escherichia coli O157

The O157:H7 lineage of enterohemorrhagic Escherichia coli is a geographically disseminated complex of highly related genotypes that share common ancestry. The common clone that is found worldwide carries several markers of events in its evolution, including markers for acquisition of virulence genes and loss of physiological characteristics, such as sorbitol fermentation ability and beta-glucuronidase production. Populations of variants that are distinct with respect to motility and the sorbitol and beta-glucuronidase markers appear to have diverged at several points along the inferred evolutionary pathway. In addition to these variants, distinct subpopulations of the contemporary non-sorbitol-fermenting, beta-glucuronidase-negative O157:H7 clone were recently detected among bovine and human clinical isolates in the United States by using high-resolution genome comparison. In order to determine if these recently described subpopulations were derived from a regional or ancestral divergence event, we used octamer-based genome scanning, marker sorting, and DNA sequence analysis to examine their phylogenetic relationship to populations of non-sorbitol-fermenting, beta-glucuronidase negative O157:H7 and O157:H- strains from Australia. The inferred phylogeny is consistent with the hypothesis that subpopulations on each continent resulted from geographic spread of an ancestral divergence event and subsequent expansion of distinct subpopulations. Marker sorting and DNA sequence analyses identified sets of monophyletic markers consistent with the pattern of divergence and demonstrated that phylogeographic variation occurred through emergence of regional subclones and concentration of regional polymorphisms among distinct subpopulations. DNA sequence analysis of representative polyphyletic markers showed that genome diversity accrued through random drift and bacteriophage-mediated events.

Investigation of bacterial diversity in the feces of cattle fed different diets.

The objective of this study is to investigate individual animal variation of bovine fecal microbiota including as affected by diets. Fecal samples were collected from 426 cattle fed 1 of 3 diets typically fed to feedlot cattle: 1) 143 steers fed finishing diet (83% dry-rolled corn, 13% corn silage, and 4% supplement), 2) 147 steers fed late growing diet (66% dry-rolled corn, 26% corn silage, and 8% supplement), and 3) 136 heifers fed early growing diet (70% corn silage and 30% alfalfa haylage). Bacterial 16S rRNA gene amplicons were determined from individual fecal samples using next-generation pyrosequencing technology. A total of 2,149,008 16S rRNA gene sequences from 333 cattle with at least 2,000 sequences were analyzed. Firmicutes and Bacteroidetes were dominant phyla in all fecal samples. At the genus level, Oscillibacter, Turicibacter, Roseburia, Fecalibacterium, Coprococcus, Clostridium, Prevotella, and Succinivibrio were represented by more than 1% of total sequences. However, numerous sequences could not be assigned to a known genus. Dominant unclassified groups were unclassified Ruminococcaceae and unclassified Lachnospiraceae that could be classified to a family but not to a genus. These dominant genera and unclassified groups differed (P < 0.001) with diets. A total of 176,692 operational taxonomic units (OTU) were identified in combination across all the 333 cattle. Only 2,359 OTU were shared across 3 diet groups. UniFrac analysis showed that bacterial communities in cattle feces were greatly affected by dietary differences. This study indicates that the community structure of fecal microbiota in cattle is greatly affected by diet, particularly between forage- and concentrate-based diets.

Bile-acid-mediated decrease in endoplasmic reticulum stress: a potential contributor to the metabolic benefits of ileal interposition surgery in UCD-T2DM rats

SUMMARY Post-operative increases in circulating bile acids have been suggested to contribute to the metabolic benefits of bariatric surgery; however, their mechanistic contributions remain undefined. We have previously reported that ileal interposition (IT) surgery delays the onset of type 2 diabetes in UCD-T2DM rats and increases circulating bile acids, independently of effects on energy intake or body weight. Therefore, we investigated potential mechanisms by which post-operative increases in circulating bile acids improve glucose homeostasis after IT surgery. IT, sham or no surgery was performed on 2-month-old weight-matched male UCD-T2DM rats. Animals underwent an oral fat tolerance test (OFTT) and serial oral glucose tolerance tests (OGTT). Tissues were collected at 1.5 and 4.5 months after surgery. Cell culture models were used to investigate interactions between bile acids and ER stress. IT-operated animals exhibited marked improvements in glucose and lipid metabolism, with concurrent increases in postprandial glucagon-like peptide-1 (GLP-1) secretion during the OFTT and OGTTs, independently of food intake and body weight. Measurement of circulating bile acid profiles revealed increases in circulating total bile acids in IT-operated animals, with a preferential increase in circulating cholic acid concentrations. Gut microbial populations were assessed as potential contributors to the increases in circulating bile acid concentrations, which revealed proportional increases in Gammaproteobacteria in IT-operated animals. Furthermore, IT surgery decreased all three sub-arms of ER stress signaling in liver, adipose and pancreas tissues. Amelioration of ER stress coincided with improved insulin signaling and preservation of β-cell mass in IT-operated animals. Incubation of hepatocyte, adipocyte and β-cell lines with cholic acid decreased ER stress. These results suggest that postoperative increases in circulating cholic acid concentration contribute to improvements in glucose homeostasis after IT surgery by ameliorating ER stress.

Microbial Successions Are Associated with Changes in Chemical Profiles of a Model Refrigerated Fresh Pork Sausage during an 80-Day Shelf Life Study

ABSTRACT Fresh pork sausage is produced without a microbial kill step and therefore chilled or frozen to control microbial growth. In this report, the microbiota in a chilled fresh pork sausage model produced with or without an antimicrobial combination of sodium lactate and sodium diacetate was studied using a combination of traditional microbiological methods and deep pyrosequencing of 16S rRNA gene amplicons. In the untreated system, microbial populations rose from 102 to 106 CFU/g within 15 days of storage at 4°C, peaking at nearly 108 CFU/g by day 30. Pyrosequencing revealed a complex community at day 0, with taxa belonging to the Bacilli, Gammaproteobacteria, Betaproteobacteria, Actinobacteria, Bacteroidetes, and Clostridia. During storage at 4°C, the untreated system displayed a complex succession, with species of Weissella and Leuconostoc that dominate the product at day 0 being displaced by species of Pseudomonas (P. lini and P. psychrophila) within 15 days. By day 30, a second wave of taxa (Lactobacillus graminis, Carnobacterium divergens, Buttiauxella brennerae, Yersinia mollaretti, and a taxon of Serratia) dominated the population, and this succession coincided with significant chemical changes in the matrix. Treatment with lactate-diacetate altered the dynamics dramatically, yielding a monophasic growth curve of a single species of Lactobacillus (L. graminis), followed by a uniform selective die-off of the majority of species in the population. Of the six species of Lactobacillus that were routinely detected, L. graminis became the dominant member in all samples, and its origins were traced to the spice blend used in the formulation.

Next-Generation Sequencing of Crown and Rhizome Transcriptome from an Upland, Tetraploid Switchgrass

The crown and rhizome transcriptome of an upland tetraploid switchgrass cultivar cv Summer well adapted to the upper Midwest was investigated using the Roche 454-FLX pyrosequencing platform. Overall, approximately one million reads consisting of 216 million bases were assembled into 27,687 contigs and 43,094 singletons. Analyses of these sequences revealed minor contamination with non-plant sequences (< 0.5%), indicating that a majority were for transcripts coded by the switchgrass genome. Blast2Gos comparisons resulted in the annotation of ~65% of the contig sequences and ~40% of the singleton sequences. Contig sequences were mostly homologous to other plant sequences, dominated by matches to Sorghum bicolor genome. Singleton sequences, while displaying significant matches to S. bicolor, also contained sequences matching non-plant species. Comparisons of the 454 dataset to existing EST collections resulted in the identification of 30,177 new sequences. These new sequences coded for a number of different proteins and a selective analysis of two categories, namely, peroxidases and transcription factors, resulted in the identification of specific peroxidases and a number of low-abundance transcription factors expected to be involved in chromatin remodeling. KEGG maps for glycolysis and sugar metabolism showed high levels of transcript coding for enzymes involved in primary metabolism. The assembly provided significant insights into the status of these tissues and broadly indicated that there was active metabolism taking place in the crown and rhizomes at post-anthesis, the seed maturation stage of plant development.

Genomic instability in regions adjacent to a highly conserved pch prophage in Escherichia coli O157:H7 generates diversity in expression patterns of the LEE pathogenicity island.

The LEE pathogenicity island has been acquired on multiple occasions within the different lineages of enteropathogenic and enterohemorrhagic Escherichia coli. In each lineage, LEE expression is regulated by complex networks of pathways, including core pathways shared by all lineages and lineage-specific pathways. Within the O157:H7 lineage of enterohemorrhagic E. coli, strain-to-strain variation in LEE expression has been observed, implying that expression patterns can diversify even within highly related subpopulations. Using comparative genomics of E. coli O157:H7 subpopulations, we have identified one source of strain-level variation affecting LEE expression. The variation occurs in prophage-dense regions of the genome that lie immediately adjacent to the late regions of the pch prophage carrying pchA, pchB, pchC, and a newly identified pch gene, pchX. Genomic segments extending from the holin S region to the pchA, pchB, pchC, and pchX genes of their respective prophage are highly conserved but are nonetheless embedded within adjacent genomic segments that are extraordinarily variable, termed pch adjacent genomic regions (pch AGR). Despite the remarkable degree of variation, the pattern of variation in pch AGR is highly correlated with the distribution of phylogenetic markers on the backbone of the genome. Quantitative analysis of transcription from the LEE1 promoter further revealed that variation in the pch AGR has substantial effects on absolute levels and patterns of LEE1 transcription. Variation in the pch AGR therefore serves as a mechanism to diversify LEE expression patterns, and the lineage-specific pattern of pch AGR variation could ultimately influence ecological or virulence characteristics of subpopulations within each lineage.