What Can We Learn about the Bias of Microbiome Studies from Analyzing Data from Mock Communities?

It is known that data from both 16S and shotgun metagenomics studies are subject to biases that cause the observed relative abundances of taxa to differ from their true values. Model community analyses, in which the relative abundances of all taxa in the sample are known by construction, seem to offer the hope that these biases can be measured. However, it is unclear whether the bias we measure in a mock community analysis is the same as we measure in a sample in which taxa are spiked in at known relative abundance, or if the biases we measure in spike-in samples is the same as the bias we would measure in a real (e.g., biological) sample. Here, we consider these questions in the context of 16S rRNA measurements on three sets of samples: the commercially available Zymo cells model community; the Zymo model community mixed with Swedish Snus, a smokeless tobacco product that is virtually bacteria-free; and a set of commercially available smokeless tobacco products. Each set of samples was subject to four different extraction protocols. The goal of our analysis is to determine whether the patterns of bias observed in each set of samples are the same, i.e., can we learn about the bias in the commercially available smokeless tobacco products by studying the Zymo cells model community?

Associations between microbial communities and key chemical constituents in U.S. domestic moist snuff.

BACKGROUND: Smokeless tobacco (ST) products are widely used throughout the world and contribute to morbidity and mortality in users through an increased risk of cancers and oral diseases. Bacterial populations in ST contribute to taste, but their presence can also create carcinogenic, Tobacco-Specific N-nitrosamines (TSNAs). Previous studies of microbial communities in tobacco products lacked chemistry data (e.g. nicotine, TSNAs) to characterize the products and identify associations between carcinogen levels and taxonomic groups. This study uses statistical analysis to identify potential associations between microbial and chemical constituents in moist snuff products. METHODS: We quantitatively analyzed 38 smokeless tobacco products for TSNAs using liquid chromatography with tandem mass spectrometry (LC-MS/MS), and nicotine using gas chromatography with mass spectrometry (GC-MS). Moisture content determinations (by weight loss on drying), and pH measurements were also performed. We used 16S rRNA gene sequencing to characterize the microbial composition, and additionally measured total 16S bacterial counts using a quantitative PCR assay. RESULTS: Our findings link chemical constituents to their associated bacterial populations. We found core taxonomic groups often varied between manufacturers. When manufacturer and flavor were controlled for as confounding variables, the genus Lactobacillus was found to be positively associated with TSNAs. while the genera Enteractinococcus and Brevibacterium were negatively associated. Three genera (Corynebacterium, Brachybacterium, and Xanthomonas) were found to be negatively associated with nicotine concentrations. Associations were also investigated separately for products from each manufacturer. Products from one manufacturer had a positive association between TSNAs and bacteria in the genus Marinilactibacillus. Additionally, we found that TSNA levels in many products were lower compared with previously published chemical surveys. Finally, we observed consistent results when either relative or absolute abundance data were analyzed, while results from analyses of log-ratio-transformed abundances were divergent.

Restoring the Duality between Principal Components of a Distance Matrix and Linear Combinations of Predictors, with Application to Studies of the Microbiome.

Appreciation of the importance of the microbiome is increasing, as sequencing technology has made it possible to ascertain the microbial content of a variety of samples. Studies that sequence the 16S rRNA gene, ubiquitous in and nearly exclusive to bacteria, have proliferated in the medical literature. After sequences are binned into operational taxonomic units (OTUs) or species, data from these studies are summarized in a data matrix with the observed counts from each OTU for each sample. Analysis often reduces these data further to a matrix of pairwise distances or dissimilarities; plotting the first two or three principal components (PCs) of this distance matrix often reveals meaningful groupings in the data. However, once the distance matrix is calculated, it is no longer clear which OTUs or species are important to the observed clustering; further, the PCs are hard to interpret and cannot be calculated for subsequent observations. We show how to construct approximate decompositions of the data matrix that pair PCs with linear combinations of OTU or species frequencies, and show how these decompositions can be used to construct biplots, select important OTUs and partition the variability in the data matrix into contributions corresponding to PCs of an arbitrary distance or dissimilarity matrix. To illustrate our approach, we conduct an analysis of the bacteria found in 45 smokeless tobacco samples.

Characterization of Bacterial Communities in Selected Smokeless Tobacco Products Using 16S rDNA Analysis.

The bacterial communities present in smokeless tobacco (ST) products have not previously reported. In this study, we used Next Generation Sequencing to study the bacteria present in U.S.-made dry snuff, moist snuff and Sudanese toombak. Sample diversity and taxonomic abundances were investigated in these products. A total of 33 bacterial families from four phyla, Actinobacteria, Firmicutes, Proteobacteria and Bacteroidetes, were identified. U.S.-produced dry snuff products contained a diverse distribution of all four phyla. Moist snuff products were dominated by Firmicutes. Toombak samples contained mainly Actinobacteria and Firmicutes (Aerococcaceae, Enterococcaceae, and Staphylococcaceae). The program PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) was used to impute the prevalence of genes encoding selected bacterial toxins, antibiotic resistance genes and other pro-inflammatory molecules. PICRUSt also predicted the presence of specific nitrate reductase genes, whose products can contribute to the formation of carcinogenic nitrosamines. Characterization of microbial community abundances and their associated genomes gives us an indication of the presence or absence of pathways of interest and can be used as a foundation for further investigation into the unique microbiological and chemical environments of smokeless tobacco products.

Differences between the normal vaginal bacterial community of baboons and that of humans.

Humans and baboons (Papio spp.) share considerable anatomical and physiological similarities in their reproductive tracts. Given the similarities, it is reasonable to expect that the normal vaginal microbial composition (microbiota) of baboons would be similar to that of humans. We have used a 16S rRNA phylogenetic approach to assess the composition of the baboon vaginal microbiota in a set of nine animals from a captive facility and six from the wild. Results show that although Gram-positive bacteria dominate in baboons as they do in humans, there are major differences between the vaginal microbiota of baboons and that of humans. In contrast to humans, the species of Gram-positive bacteria (Firmicutes) were taxa other than Lactobacillus species. In addition, some groups of Gram-negative bacteria that are not normally abundant in humans were found in the baboon samples. A further level of difference was also seen even within the same bacterial phylogenetic group, as baboon strains tended to be more phylogenetically distinct from human strains than human strains were with each other. Finally, results of our analysis suggests that co-evolution of microbes and their hosts cannot account for the major differences between the microbiota of baboons and that of humans because divergences between the major bacterial genera were too ancient to have occurred since primates evolved. Instead, the primate vaginal tracts appear to have acquired discrete subsets of bacteria from the vast diversity of bacteria available in the environment and established a community responsive to and compatible with host species physiology.

Baboon vaginal microbiota: an overlooked aspect of primate physiology.

The bacterial population of the primate vaginal canal is an infant primate's first exposure to the microbial population inhabiting the outside world. Yet, little is known about this population and the effect it might have on the development and survival of the infant primate. As a first step toward characterizing the vaginal microbiota of a nonhuman primate, we used denaturing gradient gel electrophoresis to evaluate variations in the vaginal microbiota of a group of 35 baboons (Papio hamadryas), which were housed in a facility where they shared the same diet and the same environmental conditions. We found that, despite the uniform environment, there were appreciable differences in the composition of the microbiota from one individual to another. Our results also indicate that a simple swab test is sufficient for sampling the vaginal microbiota in the field, a finding that should help make more detailed characterization of the microbiota of wild primates feasible in the future.

Antibiotic resistance genes in the vaginal microbiota of primates not normally exposed to antibiotics.

Previous studies of resistance gene ecology have focused primarily on populations such as hospital patients and farm animals that are regularly exposed to antibiotics. Also, these studies have tended to focus on numerically minor populations such as enterics or enterococci. We report here a cultivation-independent approach that allowed us to assess the presence of antibiotic resistance genes in the numerically predominant populations of the vaginal microbiota of two populations of primates that are seldom or never exposed to antibiotics: baboons and mangabeys. Most of these animals were part of a captive colony in Texas that is used for scientific studies of female physiology and physical anthropology topics. Samples from some wild baboons were also tested. Vaginal swab samples, obtained in connection with a study designed to define the normal microbiota of the female vaginal canal, were tested for the presence of two types of antibiotic resistance genes: tetracycline resistance (tet) genes and erythromycin resistance (erm) genes. These genes are frequently found in human isolates of the two types of bacteria that were a substantial part of the normal microbiota of primates (Firmicutes and Bacteroidetes). Since cultivation was not feasible, polymerase chain reaction and DNA sequencing were used to detect and characterize these resistance genes. The tet(M) and tet(W) genes were found most commonly, and the tet(Q) gene was found in over a third of the samples from baboons. The ermB and ermF genes were found only in a minority of the samples. The ermG gene was not found in any of the specimens tested. Polymerase chain reaction analysis showed that at least some tet(M) and tet(Q) genes were genetically linked to DNA from known conjugative transposons (CTns), Tn916 and CTnDOT. Our results raise questions about the extent to which extensive exposure to antibiotics is the only pressure necessary to maintain resistance genes in natural settings.