Previous Article |
Table of Contents
| Next Article
MICROBIOLOGY
Molecular analysis of the bacterial microbiota in the human stomach
, 
, 






,
, 
*Department of Microbiology and Immunology, Stanford University School of Medicine, Fairchild Science Building, 299 Campus Drive, Stanford, CA 94305;
Veterans Affairs Palo Alto Health Care System, 3801 Miranda Avenue, Palo Alto, CA 94304;
Department of Medicine, Stanford University School of Medicine, 300 Pasteur Drive, Stanford, CA 94305; ¶The Institute for Genomic Research, 9712 Medical Center Drive, Rockville, MD 20850; ||Department of Statistics, Sequoia Hall, 390 Serra Mall, Stanford University, Stanford, CA 94305; and Departments of **Medicine and 
Microbiology, New York University School of Medicine, New York, NY 10016
Edited by Jeffrey I. Gordon, Washington University School of Medicine, St. Louis, MO, and approved November 23, 2005 (received for review August 3, 2005)
| Abstract |
|---|
|
|
|---|
16S ribosomal rRNA clone library | Helicobacter pylori | microbial diversity | human indigenous microbiota
The human gastric microbiota has not been extensively characterized by using 16S rDNA sequence analysis. Most PCR-based studies of the gastric biota have focused on specific genera, e.g., Helicobacter (19) and Lactobacilli (20). In 2000, a combination of broad-range bacterial PCR and temporal temperature gradient gel electrophoresis suggested the presence of a complex gastric mucosal biota (21). However, this study relied on a small number of short (200 bp) 16S rDNA fragments and lacked the resolution to define sequence differences at the species level (only 11 genera were described). It remains unclear whether the gastric mucosal microbiota varies between different anatomic sites (fundus, corpus, antrum, and pylorus), is distinct from the oral and esophageal microbiota, and whether the presence of H. pylori influences the composition of the gastric microbial community. To address these issues, we characterized the gastric microbiota of H. pylori-positive and -negative subjects by using broad-range PCR and 16S rDNA sequence analysis.
| Results |
|---|
|
|
|---|
|
|
|
Influence of H. pylori and Anatomical Location on the Composition of the Gastric Microflora. To investigate the possible impact of H. pylori on the presence of other bacterial species, we first examined the phylogenetic tree and phylotype abundance plot (Fig. 1). Subjects in group 3 demonstrated a relative lack of non-Proteobacteria phylotypes, especially Bacteroidetes phylotypes, compared with H. pylori-negative subjects. This difference could largely be attributed to the relative low number of non-H. pylori clones among these libraries. When H. pylori sequences were left out of the analysis, phylotype evenness and diversity among libraries in which H. pylori clones were found (groups 2 and 3) were higher than that of libraries from H. pylori-negative subjects (group 1) (Table 2). Evaluation of the phylum distribution of all non-H. pylori phylotypes from each of the three groups of samples and subjects based on H. pylori status (Fig. 3) revealed no gross differences in the taxonomic patterns.
Subject and clone library variability were explored with double principal coordinate analysis (DPCoA) and redundancy analysis of non-H. pylori sequences. No significant association between phylotype distribution and H. pylori status (P > 0.5) or anatomical location (P > 0.5) was found. In addition, phylotype distribution was not significantly related to gastric pH levels (P = 0.47). Hierarchical clustering of the data did not reveal distinct clusters between H. pylori-negative and -positive subjects (Fig. 5, which is published as supporting information on the PNAS web site).
|
-LIBSHUFF revealed that the majority of the libraries were distinct from any other (Table 4, which is published as supporting information on the PNAS web site). As predicted by the DPCoA results, no pattern of library relatedness between subject groups was evident in this analysis. However, the combined gastric sequence library (with or without H. pylori clones) was distinct from libraries derived from the subgingival crevice and esophageal mucosa obtained from other studies (P < 0.00005) (P. W. Lepp, K. Palm, and K. Lightfield, unpublished results; ref. 14). Caution must be taken in interpreting the latter
-LIBSHUFF results, given that the oral and esophageal libraries were not derived from the same individuals, the subjects in these studies had different clinical syndromes, and the clone libraries were constructed by using different PCR and sequencing protocols.
|
| Discussion |
|---|
|
|
|---|
The current study used sensitive molecular methods to reveal previously uncharacterized features of the gastric bacterial biota among adults with symptomatic upper gastrointestinal disease. Interestingly, H. pylori was the most abundant phylotype in the human stomach in subjects who tested positive for this organism by using conventional clinical approaches. However, H. pylori DNA also was present in the biopsy samples of seven subjects that were considered to be H. pylori-negative by conventional methods. It is possible that some of our biopsies became cross-contaminated with H. pylori DNA during tissue processing, although abundant extraction controls did not yield visible amounts of DNA on gel electrophoresis after broad-range PCR. An alternative explanation is that the number of H. pylori bacteria in these samples (group 2) was too low to be detected by conventional methods and too few to elicit host antibody production. This hypothesis is supported by the low percentage of H. pylori clones in the libraries derived from samples and subjects that were negative by using conventional H. pylori tests. Detection of bacterial DNA does not necessarily indicate the presence of live, resident organisms. Bacterial DNA in stomach biopsy samples might reflect the presence of bacterial cell remnants or transient bacteria (3).
We found that the human gastric environment contains ribosomal DNA from a wealth of bacteria, in addition to H. pylori (Fig. 1). One hundred twenty-eight phylotypes among eight bacterial phyla were detected, suggesting much greater gastric ecosystem diversity than previously described; this higher-order taxon diversity is similar to that described for the lower gastrointestinal tract (1517). Phylotypes were found that were affiliated with bacterial taxa not yet described with cultivation-based analyses of the gastric flora, including Caulobacter, Actinobacillus, Corynebacterium, Rothia, Gemella, Leptotrichia, Porphyromonas, Capnocytophaga, TM7, Flexistipes, and Deinococcus. Among the representative sequences of the 128 phylotypes found in this study, 64 (50%) were derived from uncultivated bacteria. Of these phylotypes, 43 (67%) have been previously described in specimens from the human mouth (1013). This finding is not unexpected because it is likely that the composition of the gastric community is not only determined by niche-specific factors but also by stochastic colonization from upstream components of the alimentary tract.
Ten percent of the phylotypes found in this study were previously uncharacterized, and some of these sequences display >5% sequence dissimilarity from publicly available sequences (Table 1). Because little is known about their closest relatives, the biological and clinical significance of these putative organisms is unclear. Of the previously uncharacterized phylotypes found in this study, the sequence belonging to the Deinococcus/Thermus phylum is particularly interesting. Deinococcus species (specifically D. radiodurans) have been isolated from extreme environments, such as radioactive waste disposal sites and hot springs, as well as from animal fecal samples (23). To our knowledge, this Deinococcus-related sequence is the first identified from a human.
-LIBSHUFF analysis verified a large degree of intersubject variability, similar to that observed in other 16S rDNA surveys of the gastrointestinal microbiota (1517). Given the predominance of men in this study, comparison of the gastric community of the only woman with those of the 22 male subjects would not be helpful. Based on visual inspection of the phylogenetic trees (Fig. 1), assessment of higher-order taxa distribution among samples (Fig. 2), and DPCoA results, the presence of H. pylori, gastric anatomical location, and gastric pH level did not significantly influence or alter local bacterial community composition. In addition, the presence of H. pylori did not have a negative effect on diversity and evenness of the other members of the gastric community (Table 2). Further molecular analyses are needed to understand better the impact of H. pylori infection on the gastric mucosal microbial community.
In this study, bacterial sequences in the stomach were not simply a random sampling of bacterial sequences from oral (subgingival) or esophageal communities. Given prior assumptions about limited microbial diversity in the human stomach, and the natural flow of oroesophageal secretions into the stomach, these findings are somewhat surprising and suggest the presence of distinct bacterial communities that have adapted to multiple specific environmental habitats in the stomach. These results must be interpreted with caution, because the libraries were derived from different subjects with differing medical conditions and constructed by using different DNA extraction methods and amplification protocols. Whether the stomach is home to distinct bacterial communities remains to be verified. Further studies should include a more detailed spatial (anatomical) and temporal analysis of gastric microbial community structure within subjects of different ethnic and racial backgrounds, as well as an analysis of the effect of H. pylori eradication therapy on community structure. A better understanding of the resident microbial communities at healthy and diseased sites should shed light on the pathogenesis, diagnosis, and treatment of gastrointestinal illnesses.
| Materials and Methods |
|---|
|
|
|---|
Construction of 16S rDNA Clone Libraries. Total genomic DNA was extracted from the biopsy samples by using a combination of the QIAamp DNA isolation kit (Qiagen, Valencia, CA) and a bead-beater method. Universal bacterial PCR primers were used to amplify the region corresponding to positions 8 to 806 of the Escherichia coli 16S rRNA gene by using a 25-cycle PCR. Experimental details are provided in Supporting Materials and Methods.
Amplification products were gel-purified, ligated into the pCR2.1 vector, and transformed into E. coli TOP10 cells by using the TOPO-TA cloning kit (Invitrogen) according to the manufacturer's instructions. From each gastric biopsy, 95 transformants were selected. Plasmid template DNA from each transformant was prepared by a modified alkaline lysis method. Sequencing reactions were carried out with BigDye terminator and run on ABI 3730xl sequencers (Applied Biosystems, Foster City, CA).
Phylogenetic Analysis. A total of 1,833 nonchimeric 16S rDNA sequences were aligned and analyzed by using the ARB software package (24); phylogenetic analysis details are available in Supporting Materials and Methods. A neighbor-joining tree containing all sequences was constructed by using Olsen correction and 738 nonambiguous filter positions. Phylotypes were defined as single clones or groups of sequences having 99% or more identity to a published sequence, or, in the case of a previously uncharacterized phylotype, <99% identity to any public database sequence. The reliability of phylogenetic tree topology was tested by bootstrap resampling with 100 replicates.
Estimation of Microbial Diversity. Rarefaction curves, richness estimations, and diversity indices were determined with ESTIMATES (Version 7, http://viceroy.eeb.uconn.edu/estimates). The species rarefaction curve of the entire data set was computed by using the individual-based Coleman method and the sample-based Mao Tau method. The bias-corrected Chao1 estimator of species richness was calculated after 1,000 randomizations of sampling without replacement. The percentage of coverage was calculated by Good's method with the formula [1 (n/N)] x 100, where n is the number of phylotypes in a sample represented by one clone (singletons) and N is the total number of sequences in that sample (25). Diversity of the sampled sequence set was estimated by using the Simpson and Shannon indices within the ESTIMATES application. The Shannon index of evenness was calculated with the formula E = eD/N, where D is the Shannon diversity index.
Statistical Analysis. DPCoA was implemented to investigate the relationships between phylotype dissimilarities and Rao diversity data among the gastric communities (26) by using the ade4 and dpcoa functions within the R statistical package (www.rproject.org). Redundancy analysis (1,000 permutations) applied to the results of DPCoA was used to evaluate the statistical significance of categorical explanatory variables, such as positive or negative H. pylori status (27). Cluster analysis methods were used to analyze the relationships between sequence phylotypes as detected among each subject. More details are provided in Supporting Materials and Methods.
To determine whether differences in clone libraries were due to underlying variability in the microbial communities versus artifacts of subsampling, gastric clone libraries were compared with the
-LIBSHUFF program (28, 29). This tool compares homologous and heterologous coverage curves by using the integral form of the Cramér-von Mises statistic and performs multiple pairwise comparisons among a set of libraries (details in Supporting Materials and Methods).
-LIBSHUFF was used to compare individual clone libraries and combined clone libraries from the stomach, the oral cavity, and esophagus. Regarding the latter, we compared our pooled gastric library to 1,563 subgingival 16S rDNA sequences derived from subjects with and without gingivitis or periodontitis (P. W. Lepp, K. Palm, and K. Lightfield, unpublished results) and 896 esophageal 16S rDNA sequences from subjects with normal esophageal histology (14).
| Acknowledgements |
|---|
| Footnotes |
|---|
This paper was submitted directly (Track II) to the PNAS office.
Freely available online through the PNAS open access option.
Abbreviation: DPCoA, double principal coordinate analysis.
Data deposition: The 16S rDNA sequences of clones have been deposited in the GenBank database (accession nos. AY582885AY582898).
To whom correspondence may be addressed. E-mail: eliesbik{at}stanford.edu or relman{at}stanford.edu.
© 2006 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
M. J Blaser, Y. Chen, and J. Reibman Does Helicobacter pylori protect against asthma and allergy? Gut, May 1, 2008; 57(5): 561 - 567. [Full Text] [PDF] |
||||
![]() |
C.-W. Lee, B. Rickman, A. B. Rogers, Z. Ge, T. C. Wang, and J. G. Fox Helicobacter pylori Eradication Prevents Progression of Gastric Cancer in Hypergastrinemic INS-GAS Mice Cancer Res., May 1, 2008; 68(9): 3540 - 3548. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Frank, C. I. Reich, S. Sharma, J. S. Weisbaum, B. A. Wilson, and G. J. Olsen Critical Evaluation of Two Primers Commonly Used for Amplification of Bacterial 16S rRNA Genes Appl. Envir. Microbiol., April 15, 2008; 74(8): 2461 - 2470. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. H. Liu, S. M. Lee, J. M. VanLare, D. L. Kasper, and S. K. Mazmanian From the Cover: Regulation of surface architecture by symbiotic bacteria mediates host colonization PNAS, March 11, 2008; 105(10): 3951 - 3956. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Levine, E. A. Lin, W. Emara, J. Kang, M. DiBenedetto, T. Ando, D. Falush, and M. J. Blaser Plastic cells and populations: DNA substrate characteristics in Helicobacter pylori transformation define a flexible but conservative system for genomic variation FASEB J, November 1, 2007; 21(13): 3458 - 3467. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A. Rader, S. R. Campagna, M. F. Semmelhack, B. L. Bassler, and K. Guillemin The Quorum-Sensing Molecule Autoinducer 2 Regulates Motility and Flagellar Morphogenesis in Helicobacter pylori J. Bacteriol., September 1, 2007; 189(17): 6109 - 6117. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Fraune and T. C. G. Bosch Long-term maintenance of species-specific bacterial microbiota in the basal metazoan Hydra PNAS, August 7, 2007; 104(32): 13146 - 13151. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Chen and M. J. Blaser Inverse Associations of Helicobacter pylori With Asthma and Allergy Arch Intern Med, April 23, 2007; 167(8): 821 - 827. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. O. Glazunova, T. Launay, D. Raoult, and V. Roux Prevotella timonensis sp. nov., isolated from a human breast abscess Int J Syst Evol Microbiol, April 1, 2007; 57(4): 883 - 886. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. R. Cox and M. S. Gilmore Native Microbial Colonization of Drosophila melanogaster and Its Use as a Model of Enterococcus faecalis Pathogenesis Infect. Immun., April 1, 2007; 75(4): 1565 - 1576. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Gao, C.-h. Tseng, Z. Pei, and M. J. Blaser Molecular analysis of human forearm superficial skin bacterial biota PNAS, February 20, 2007; 104(8): 2927 - 2932. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. R. Hoffman, E. Deziel, D. A. D'Argenio, F. Lepine, J. Emerson, S. McNamara, R. L. Gibson, B. W. Ramsey, and S. I. Miller Selection for Staphylococcus aureus small-colony variants due to growth in the presence of Pseudomonas aeruginosa PNAS, December 26, 2006; 103(52): 19890 - 19895. [Abstract] [Full Text] [PDF] |
||||
![]() |
W.-K. Lee, K. Ogura, J. T. Loh, T. L. Cover, and D. E. Berg Quantitative Effect of luxS Gene Inactivation on the Fitness of Helicobacter pylori Appl. Envir. Microbiol., October 1, 2006; 72(10): 6615 - 6622. [Abstract] [Full Text] [PDF] |
||||
![]() |
M.-E. Dumas, R. H. Barton, A. Toye, O. Cloarec, C. Blancher, A. Rothwell, J. Fearnside, R. Tatoud, V. Blanc, J. C. Lindon, et al. Metabolic profiling reveals a contribution of gut microbiota to fatty liver phenotype in insulin-resistant mice PNAS, August 15, 2006; 103(33): 12511 - 12516. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Oh, H. Kling-Backhed, M. Giannakis, J. Xu, R. S. Fulton, L. A. Fulton, H. S. Cordum, C. Wang, G. Elliott, J. Edwards, et al. The complete genome sequence of a chronic atrophic gastritis Helicobacter pylori strain: Evolution during disease progression PNAS, June 27, 2006; 103(26): 9999 - 10004. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||