Pyrosequencing of chaperonin-60 (cpn60) amplicons as a means of determining microbial community composition

Human Microbiomes and Disease for the Biomedical Data Scientist

The human microbiome is complex, variable from person to person, essential for health, and related to both the risk for disease and the efficacy of our treatments. There are robust techniques to describe microbiota with high-throughput sequencing, and there are hundreds of thousands of already-sequenced specimens in public archives. The promise remains to use the microbiome both as a prognostic factor and as a target for precision medicine. However, when used as an input in biomedical data science modeling, the microbiome presents unique challenges. Here, we review the most common techniques used to describe microbial communities, explore these unique challenges, and discuss the more successful approaches for biomedical data scientists seeking to use the microbiome as an input in their studies.

The Rhizobial Microbiome from the Tropical Savannah Zones in Northern Côte d'Ivoire

Over the past decade, many projects have been initiated worldwide to decipher the composition and function of the soil microbiome, including the African Soil Microbiome (AfSM) project that aims at providing new insights into the presence and distribution of key groups of soil bacteria from across the African continent. In this national study, carried out under the auspices of the AfSM project, we assessed the taxonomy, diversity and distribution of rhizobial genera in soils from the tropical savannah zones in Northern Côte d’Ivoire. Genomic DNA extracted from seven sampled soils was analyzed by sequencing the V4-V5 variable region of the 16S rDNA using Illumina’s MiSeq platform. Subsequent bioinformatic and phylogenetic analyses showed that these soils harbored 12 out of 18 genera of Proteobacteria harboring rhizobia species validly published to date and revealed for the first time that the Bradyrhizobium genus dominates in tropical savannah soils, together with Microvirga and Paraburkholderia. In silico comparisons of different 16S rRNA gene variable regions suggested that the V5-V7 region could be suitable for differentiating rhizobia at the genus level, possibly replacing the use of the V4-V5 region. These data could serve as indicators for future rhizobial microbiome explorations and for land-use decision-making.

Ultrahigh-Throughput Multiplexing and Sequencing of >500-Base-Pair Amplicon Regions on the Illumina HiSeq 2500 Platform

Amplicon sequencing has become a popular and widespread tool for surveying microbial communities. Lower overall costs associated with high-throughput sequencing have made it a widely adopted approach, especially for projects that necessitate sample multiplexing to eliminate batch effect and reduced time to acquire data. The method for amplicon sequencing on the Illumina HiSeq 2500 platform described here provides improved multiplexing capabilities while simultaneously producing greater quality sequence data and lower per-sample cost relative to those of the Illumina MiSeq platform without sacrificing amplicon length. To make this method more flexible for various amplicon-targeted regions as well as improve amplification from low-biomass samples, we also present and validate a 2-step PCR library preparation method.

Amplification, sequencing, and analysis of the 16S rRNA gene affords characterization of microbial community composition. As this tool has become more popular and amplicon-sequencing applications have grown in the total number of samples, growth in sample multiplexing is becoming necessary while maintaining high sequence quality and sequencing depth. Here, modifications to the Illumina HiSeq 2500 platform are described which produce greater multiplexing capabilities and 300-bp paired-end reads of higher quality than those produced by the current Illumina MiSeq platform. To improve the feasibility and flexibility of this method, a 2-step PCR amplification protocol is also described that allows for targeting of different amplicon regions, and enhances amplification success from samples with low bacterial bioburden.

IMPORTANCE Amplicon sequencing has become a popular and widespread tool for surveying microbial communities. Lower overall costs associated with high-throughput sequencing have made it a widely adopted approach, especially for projects that necessitate sample multiplexing to eliminate batch effect and reduced time to acquire data. The method for amplicon sequencing on the Illumina HiSeq 2500 platform described here provides improved multiplexing capabilities while simultaneously producing greater quality sequence data and lower per-sample cost relative to those of the Illumina MiSeq platform without sacrificing amplicon length. To make this method more flexible for various amplicon-targeted regions as well as improve amplification from low-biomass samples, we also present and validate a 2-step PCR library preparation method.

Detection and Enumeration of Spore-Forming Bacteria in Powdered Dairy Products

With the abolition of milk quotas in the European Union in 2015, several member states including Ireland, Luxembourg, and Belgium have seen year on year bi-monthly milk deliveries to dairies increase by up to 35%. Milk production has also increased outside of Europe in the past number of years. Unsurprisingly, there has been a corresponding increased focus on the production of dried milk products for improved shelf life. These powders are used in a wide variety of products, including confectionery, infant formula, sports dietary supplements and supplements for health recovery. To ensure quality and safety standards in the dairy sector, strict controls are in place with respect to the acceptable quantity and species of microorganisms present in these products. A particular emphasis on spore-forming bacteria is necessary due to their inherent ability to survive extreme processing conditions. Traditional microbiological detection methods used in industry have limitations in terms of time, efficiency, accuracy, and sensitivity. The following review will explore the common spore-forming bacterial contaminants of milk powders, will review the guidelines with respect to the acceptable limits of these microorganisms and will provide an insight into recent advances in methods for detecting these microbes. The various advantages and limitations with respect to the application of these diagnostics approaches for dairy food will be provided. It is anticipated that the optimization and application of these methods in appropriate ways can ensure that the enhanced pressures associated with increased production will not result in any lessening of safety and quality standards.

High Diversity and Variability in the Vaginal Microbiome in Women following Preterm Premature Rupture of Membranes (PPROM): A Prospective Cohort Study

Objective

To characterize the vaginal microbiota of women following preterm premature rupture of membranes (PPROM), and determine if microbiome composition predicts latency duration and perinatal outcomes.

Design

A prospective cohort study

Setting

Canada

Population

Women with PPROM between 24⁺⁰ and 33⁺⁶ weeks gestational age (GA).

Methods

Microbiome profiles, based on pyrosequencing of the cpn60 universal target, were generated from vaginal samples at time of presentation with PPROM, weekly thereafter, and at delivery.

Main Outcome Measures

Vaginal microbiome composition, latency duration, gestational age at delivery, perinatal outcomes.

Results

Microbiome profiles were generated from 70 samples from 36 women. Mean GA at PPROM was 28.8 wk (mean latency 2.7 wk). Microbiome profiles were highly diverse but sequences representing Megasphaera type 1 and Prevotella spp. were detected in all vaginal samples. Only 13/70 samples were dominated by Lactobacillus spp. Microbiome profiles at the time of membrane rupture did not cluster by gestational age at PPROM, latency duration, presence of chorioamnionitis or by infant outcomes. Mycoplasma and/or Ureaplasma were detected by PCR in 81% (29/36) of women, and these women had significantly lower GA at delivery and correspondingly lower birth weight infants than Mycoplasma and/or Ureaplasma negative women.

Conclusion

Women with PPROM had mixed, abnormal vaginal microbiota but the microbiome profile at PPROM did not correlate with latency duration. Prevotella spp. and Megasphaera type I were ubiquitous. The presence of Mollicutes in the vaginal microbiome was associated with lower GA at delivery. The microbiome was remarkably unstable during the latency period.

Comparative Genomics of cpn60-Defined Enterococcus hirae Ecotypes and Relationship of Gene Content Differences to Competitive Fitness

Natural microbial communities undergo selection-driven succession with changes in environmental conditions and available nutrients. In a previous study of the pig faecal Enterococcus community, we demonstrated that cpn60 universal target (UT) sequences could resolve phenotypically and genotypically distinct ecotypes of Enterococcus spp. that emerged over time in the faecal microbiome of growing pigs. In this study, we characterized genomic diversity in the identified Enterococcus hirae ecotypes in order to define further the nature and degree of genome content differences between taxa resolved by cpn60 UT sequences. Genome sequences for six representative isolates (two from each of three ecotypes) were compared. Differences in phosphotransferase systems and amino acid metabolism pathways for glutamine, proline and selenocysteine were observed. Differences in the lac family phosphotransferase system corresponded to lactose utilization phenotypes of the isolates. Competitive fitness of the E. hirae ecotypes was evaluated by in vitro growth competition assays in pig faecal extract medium. Isolates from E. hirae-1 and E. hirae-2 ecotypes were able to out-compete isolates from the E. hirae-3 ecotype, consistent with the relatively low abundance of E. hirae-3 relative to E. hirae-1 and E. hirae-2 previously observed in the pig faecal microbiome, and with observed differences between the ecotypes in gene content related to biosynthetic capacity. Results of this study provide a genomic basis for the definition of ecotypes within E. hirae and confirm the utility of the cpn60 UT sequence for high-resolution profiling of complex microbial communities.

A comprehensive method for amplicon-based and metagenomic characterization of viruses, bacteria, and eukaryotes in freshwater samples

BACKGROUND

Studies of environmental microbiota typically target only specific groups of microorganisms, with most focusing on bacteria through taxonomic classification of 16S rRNA gene sequences. For a more holistic understanding of a microbiome, a strategy to characterize the viral, bacterial, and eukaryotic components is necessary.

RESULTS

We developed a method for metagenomic and amplicon-based analysis of freshwater samples involving the concentration and size-based separation of eukaryotic, bacterial, and viral fractions. Next-generation sequencing and culture-independent approaches were used to describe and quantify microbial communities in watersheds with different land use in British Columbia. Deep amplicon sequencing was used to investigate the distribution of certain viruses (g23 and RdRp), bacteria (16S rRNA and cpn60), and eukaryotes (18S rRNA and ITS). Metagenomic sequencing was used to further characterize the gene content of the bacterial and viral fractions at both taxonomic and functional levels.

CONCLUSION

This study provides a systematic approach to separate and characterize eukaryotic-, bacterial-, and viral-sized particles. Methodologies described in this research have been applied in temporal and spatial studies to study the impact of land use on watershed microbiomes in British Columbia.

A Study of the Vaginal Microbiome in Healthy Canadian Women Utilizing cpn60-Based Molecular Profiling Reveals Distinct Gardnerella Subgroup Community State Types

The vaginal microbiota is important in women’s reproductive and overall health. However, the relationships between the structure, function and dynamics of this complex microbial community and health outcomes remain elusive. The objective of this study was to determine the phylogenetic range and abundance of prokaryotes in the vaginal microbiota of healthy, non-pregnant, ethnically diverse, reproductive-aged Canadian women. Socio-demographic, behavioural and clinical data were collected and vaginal swabs were analyzed from 310 women. Detailed profiles of their vaginal microbiomes were generated by pyrosequencing of the chaperonin-60 universal target. Six community state types (CST) were delineated by hierarchical clustering, including three Lactobacillus-dominated CST (L. crispatus, L. iners, L. jensenii), two Gardnerella-dominated (subgroups A and C) and an “intermediate” CST which included a small number of women with microbiomes dominated by seven other species or with no dominant species but minority populations of Streptococcus, Staphylococcus, Peptoniphilus, E. coli and various Proteobacteria in co-dominant communities. The striking correspondence between Nugent score and deep sequencing CST continues to reinforce the basic premise provided by the simpler Gram stain method, while additional analyses reveal detailed cpn60-based phylogeny and estimated abundance in microbial communities from vaginal samples. Ethnicity was the only demographic or clinical characteristic predicting CST, with differences in Asian and White women (p = 0.05). In conclusion, this study confirms previous work describing four cpn60-based subgroups of Gardnerella, revealing previously undescribed CST. The data describe the range of bacterial communities seen in Canadian women presenting with no specific vaginal health concerns, and provides an important baseline for future investigations of clinically important cohorts.

Optimizing a PCR protocol for cpn60-based microbiome profiling of samples variously contaminated with host genomic DNA

BACKGROUND

The current recommended protocol for chaperonin-60 (cpn60) universal target based microbiome profiling includes universal PCR of microbiome samples across an annealing temperature gradient to maximize the diversity of sequences amplified. However, the value of including this gradient approach has not been formally evaluated since the optimization of a modified universal PCR primer cocktail for cpn60 PCR. PCR conditions that maximize representation of the microbiome while minimizing PCR-associated distortion of the community structure, especially in samples containing large amounts of host genomic DNA are critical. The goal of this study was to measure the effects of PCR annealing temperature and the ratio of host to bacterial DNA on the outcome of microbiota analysis, using pig microbiota as a model environment.

FINDINGS

Six samples were chosen with an anticipated range of ratios of pig to bacterial genomic DNA, and universal cpn60 PCR amplification with an annealing temperature gradient was used to create libraries for pyrosequencing, resulting in 426,477 sequences from the six samples. The sequences obtained were classified as target (cpn60) or non-target based on the percent identity of their closest match to the cpnDB reference database, and target sequences were further processed to create microbiome profiles for each sample at each annealing temperature. Annealing temperature affected the amount of PCR product generated, with more product generated at higher temperatures. Samples containing proportionally more host genomic DNA yielded more non-target reads, especially at lower annealing temperatures. However, microbiome composition for each sample across the annealing temperature gradient remained consistent at both the phylum and operational taxonomic unit levels. Although some microbial sequences were detected at only one annealing temperature, these sequences accounted for a minority of the total microbiome.

CONCLUSIONS

These results indicate that PCR annealing temperature does have an affect on cpn60 based microbiome profiles, but that most of the differences are due to differences in detection of low abundance sequences. Higher annealing temperatures resulted in larger amounts of PCR product and lower amounts of non-target sequence amplification, especially in samples containing proportionally large amounts of host DNA. Taken together these results provide important information to guide decisions about experimental design for cpn60 based microbiome studies.

Characterization of the fecal microbiota of pigs before and after inoculation with "Brachyspira hampsonii"

“Brachyspira hampsonii” causes disease indistinguishable from swine dysentery, and the structure of the intestinal microbiome likely plays a role in determining susceptibility of individual pigs to infection and development of clinical disease. The objectives of the current study were to determine if the pre-inoculation fecal microbiota differed between inoculated pigs that did (INOC MH) or did not (INOC non-MH) develop mucohaemorrhagic diarrhea following challenge with “B. hampsonii”, and to quantify changes in the structure of the microbiome following development of clinical disease. Fecal microbiota profiles were generated based on amplification and sequencing of the cpn60 universal target sequence from 89 samples from 18 pigs collected at −8, −5, −3 and 0 days post-inoculation, and at termination. No significant differences in richness, diversity or taxonomic composition distinguished the pre-inoculation microbiomes of INOC MH and INOC non-MH pigs. However, the development of bloody diarrhea in inoculated pigs was associated with perturbation of the microbiota relative to INOC non-MH or sham-inoculated control pigs. Specifically, the fecal microbiota of INOC MH pigs was less dense (fewer total 16S rRNA copies per gram of feces), and had a lower Bacteroidetes:Firmicutes ratio. Further investigation of the potential long-term effects of Brachyspira disease on intestinal health and performance is warranted.

Simultaneous profiling of seed-associated bacteria and fungi reveals antagonistic interactions between microorganisms within a shared epiphytic microbiome on Triticum and Brassica seeds

In order to address the hypothesis that seeds from ecologically and geographically diverse plants harbor characteristic epiphytic microbiota, we characterized the bacterial and fungal microbiota associated with Triticum and Brassica seed surfaces. The total microbial complement was determined by amplification and sequencing of a fragment of chaperonin 60 (cpn60). Specific microorganisms were quantified by qPCR. Bacteria and fungi corresponding to operational taxonomic units (OTU) that were identified in the sequencing study were isolated and their interactions examined. A total of 5477 OTU were observed from seed washes. Neither total epiphytic bacterial load nor community richness/evenness was significantly different between the seed types; 578 OTU were shared among all samples at a variety of abundances. Hierarchical clustering revealed that 203 were significantly different in abundance on Triticum seeds compared with Brassica. Microorganisms isolated from seeds showed 99-100% identity between the cpn60 sequences of the isolates and the OTU sequences from this shared microbiome. Bacterial strains identified as Pantoea agglomerans had antagonistic properties toward one of the fungal isolates (Alternaria sp.), providing a possible explanation for their reciprocal abundances on both Triticum and Brassica seeds. cpn60 enabled the simultaneous profiling of bacterial and fungal microbiota and revealed a core seed-associated microbiota shared between diverse plant genera.

Characterization of the vaginal microbiota of healthy Canadian women through the menstrual cycle

BACKGROUND

The vaginal microbial community plays a vital role in maintaining women's health. Understanding the precise bacterial composition is challenging because of the diverse and difficult-to-culture nature of many bacterial constituents, necessitating culture-independent methodology. During a natural menstrual cycle, physiological changes could have an impact on bacterial growth, colonization, and community structure. The objective of this study was to assess the stability of the vaginal microbiome of healthy Canadian women throughout a menstrual cycle by using cpn60-based microbiota analysis. Vaginal swabs from 27 naturally cycling reproductive-age women were collected weekly through a single menstrual cycle. Polymerase chain reaction (PCR) was performed to amplify the universal target region of the cpn60 gene and generate amplicons representative of the microbial community. Amplicons were pyrosequenced, assembled into operational taxonomic units, and analyzed. Samples were also assayed for total 16S rRNA gene content and Gardnerella vaginalis by quantitative PCR and screened for the presence of Mollicutes by using family and genus-specific PCR.

RESULTS

Overall, the vaginal microbiome of most women remained relatively stable throughout the menstrual cycle, with little variation in diversity and only modest fluctuations in species richness. Microbiomes between women were more different than were those collected consecutively from individual women. Clustering of microbial profiles revealed the expected groupings dominated by Lactobacillus crispatus, Lactobacillus iners, and Lactobacillus jensenii. Interestingly, two additional clusters were dominated by either Bifidobacterium breve or a heterogeneous mixture of nonlactobacilli. Direct G. vaginalis quantification correlated strongly with its pyrosequencing-read abundance, and Mollicutes, including Mycoplasma hominis, Ureaplasma parvum, and Ureaplasma urealyticum, were detected in most samples.

CONCLUSIONS

Our cpn60-based investigation of the vaginal microbiome demonstrated that in healthy women most vaginal microbiomes remained stable through their menstrual cycle. Of interest in these findings was the presence of Bifidobacteriales beyond just Gardnerella species. Bifidobacteriales are frequently underrepresented in 16S rRNA gene-based studies, and their detection by cpn60-based investigation suggests that their significance in the vaginal community may be underappreciated.

Interactions in the microbiome: communities of organisms and communities of genes

A central challenge in microbial community ecology is the delineation of appropriate units of biodiversity, which can be taxonomic, phylogenetic, or functional in nature. The term 'community' is applied ambiguously; in some cases, the term refers simply to a set of observed entities, while in other cases, it requires that these entities interact with one another. Microorganisms can rapidly gain and lose genes, potentially decoupling community roles from taxonomic and phylogenetic groupings. Trait-based approaches offer a useful alternative, but many traits can be defined based on gene functions, metabolic modules, and genomic properties, and the optimal set of traits to choose is often not obvious. An analysis that considers taxon assignment and traits in concert may be ideal, with the strengths of each approach offsetting the weaknesses of the other. Individual genes also merit consideration as entities in an ecological analysis, with characteristics such as diversity, turnover, and interactions modeled using genes rather than organisms as entities. We identify some promising avenues of research that are likely to yield a deeper understanding of microbial communities that shift from observation-based questions of 'Who is there?' and 'What are they doing?' to the mechanistically driven question of 'How will they respond?'

mPUMA: a computational approach to microbiota analysis by de novo assembly of operational taxonomic units based on protein-coding barcode sequences

BACKGROUND

Formation of operational taxonomic units (OTU) is a common approach to data aggregation in microbial ecology studies based on amplification and sequencing of individual gene targets. The de novo assembly of OTU sequences has been recently demonstrated as an alternative to widely used clustering methods, providing robust information from experimental data alone, without any reliance on an external reference database.

RESULTS

Here we introduce mPUMA (microbial Profiling Using Metagenomic Assembly, http://mpuma.sourceforge.net), a software package for identification and analysis of protein-coding barcode sequence data. It was developed originally for Cpn60 universal target sequences (also known as GroEL or Hsp60). Using an unattended process that is independent of external reference sequences, mPUMA forms OTUs by DNA sequence assembly and is capable of tracking OTU abundance. mPUMA processes microbial profiles both in terms of the direct DNA sequence as well as in the translated amino acid sequence for protein coding barcodes. By forming OTUs and calculating abundance through an assembly approach, mPUMA is capable of generating inputs for several popular microbiota analysis tools. Using SFF data from sequencing of a synthetic community of Cpn60 sequences derived from the human vaginal microbiome, we demonstrate that mPUMA can faithfully reconstruct all expected OTU sequences and produce compositional profiles consistent with actual community structure.

CONCLUSIONS

mPUMA enables analysis of microbial communities while empowering the discovery of novel organisms through OTU assembly.

Characterization of the upper respiratory tract microbiomes of patients with pandemic H1N1 influenza

The upper respiratory tract microbiome has an important role in respiratory health. Influenza A is a common viral infection that challenges that health, and a well-recognized sequela is bacterial pneumonia. Given this connection, we sought to characterize the upper respiratory tract microbiota of individuals suffering from the pandemic H1N1 influenza A outbreak of 2009 and determine if microbiome profiles could be correlated with patient characteristics. We determined the microbial profiles of 65 samples from H1N1 patients by cpn60 universal target amplification and sequencing. Profiles were examined at the phylum and nearest neighbor “species” levels using the characteristics of patient gender, age, originating health authority, sample type and designation (STAT/non-STAT). At the phylum level, Actinobacteria-, Firmicutes- and Proteobacteria-dominated microbiomes were observed, with none of the patient characteristics showing significant profile composition differences. At the nearest neighbor “species” level, the upper respiratory tract microbiomes were composed of 13-20 “species” and showed a trend towards increasing diversity with patient age. Interestingly, at an individual level, most patients had one to three organisms dominant in their microbiota. A limited number of discrete microbiome profiles were observed, shared among influenza patients regardless of patient status variables. To assess the validity of analyses derived from sequence read abundance, several bacterial species were quantified by quantitative PCR and compared to the abundance of cpn60 sequence read counts obtained in the study. A strong positive correlation between read abundance and absolute bacterial quantification was observed. This study represents the first examination of the upper respiratory tract microbiome using a target other than the 16S rRNA gene and to our knowledge, the first thorough examination of this microbiome during a viral infection.

The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses

16S ribosomal RNA currently represents the most important target of study in bacterial ecology. Its use for the description of bacterial diversity is, however, limited by the presence of variable copy numbers in bacterial genomes and sequence variation within closely related taxa or within a genome. Here we use the information from sequenced bacterial genomes to explore the variability of 16S rRNA sequences and copy numbers at various taxonomic levels and apply it to estimate bacterial genome and DNA abundances. In total, 7,081 16S rRNA sequences were in silico extracted from 1,690 available bacterial genomes (1-15 per genome). While there are several phyla containing low 16S rRNA copy numbers, in certain taxa, e.g., the Firmicutes and Gammaproteobacteria, the variation is large. Genome sizes are more conserved at all tested taxonomic levels than 16S rRNA copy numbers. Only a minority of bacterial genomes harbors identical 16S rRNA gene copies, and sequence diversity increases with increasing copy numbers. While certain taxa harbor dissimilar 16S rRNA genes, others contain sequences common to multiple species. Sequence identity clusters (often termed operational taxonomic units) thus provide an imperfect representation of bacterial taxa of a certain phylogenetic rank. We have demonstrated that the information on 16S rRNA copy numbers and genome sizes of genome-sequenced bacteria may be used as an estimate for the closest related taxon in an environmental dataset to calculate alternative estimates of the relative abundance of individual bacterial taxa in environmental samples. Using an example from forest soil, this procedure would increase the abundance estimates of Acidobacteria and decrease these of Firmicutes. Using the currently available information, alternative estimates of bacterial community composition may be obtained in this way if the variation of 16S rRNA copy numbers among bacteria is considered.

The chaperonin-60 universal target is a barcode for bacteria that enables de novo assembly of metagenomic sequence data

Barcoding with molecular sequences is widely used to catalogue eukaryotic biodiversity. Studies investigating the community dynamics of microbes have relied heavily on gene-centric metagenomic profiling using two genes (16S rRNA and cpn60) to identify and track Bacteria. While there have been criteria formalized for barcoding of eukaryotes, these criteria have not been used to evaluate gene targets for other domains of life. Using the framework of the International Barcode of Life we evaluated DNA barcodes for Bacteria. Candidates from the 16S rRNA gene and the protein coding cpn60 gene were evaluated. Within complete bacterial genomes in the public domain representing 983 species from 21 phyla, the largest difference between median pairwise inter- and intra-specific distances (“barcode gap”) was found from cpn60. Distribution of sequence diversity along the ∼555 bp cpn60 target region was remarkably uniform. The barcode gap of the cpn60 universal target facilitated the faithful de novo assembly of full-length operational taxonomic units from pyrosequencing data from a synthetic microbial community. Analysis supported the recognition of both 16S rRNA and cpn60 as DNA barcodes for Bacteria. The cpn60 universal target was found to have a much larger barcode gap than 16S rRNA suggesting cpn60 as a preferred barcode for Bacteria. A large barcode gap for cpn60 provided a robust target for species-level characterization of data. The assembly of consensus sequences for barcodes was shown to be a reliable method for the identification and tracking of novel microbes in metagenomic studies.

Resolution and characterization of distinct cpn60-based subgroups of Gardnerella vaginalis in the vaginal microbiota

Bacterial vaginosis (BV), characterized by a shift of the vaginal microbiota from a Lactobacillus-dominated community to a dense biofilm containing a complex mixture of organisms, is an important risk factor in poor reproductive health outcomes. The Nugent score, based on Gram stain, is used to diagnose BV and Gardnerella vaginalis abundance in the sample is one factor determining Nugent score. A high Nugent score is indicative of BV but does not always correspond to the presence of clinical symptoms. G. vaginalis is recognized as a heterogeneous group of organisms, which can also be part of the normal, healthy vaginal microbiome. In addition, asymptomatic BV and non-Gardnerella types of BV are being recognized. In an attempt to resolve the heterogeneous group of G. vaginalis, a phylogenetic tree of cpn60 universal target sequences from G. vaginalis isolates was constructed that indicates the existence of four subgroups of G. vaginalis. This subdivision, supported by whole genome similarity calculation of representative strains using JSpecies, demonstrates that these subgroups may represent different species. The cpn60 subgroupings did not correspond with the Piot biotyping scheme, but did show consistency with ARDRA genotyping and sialidase gene presence. Isolates from all four subgroups produced biofilm in vitro. We also investigated the distribution of G. vaginalis subgroups in vaginal samples from Kenyan women with Nugent scores consistent with BV, Intermediate and Normal microbiota (n = 44). All subgroups of G. vaginalis were detected in these women, with a significant difference (z = −3.372, n = 39, p = 0.001) in frequency of G. vaginalis subgroup B between BV and Normal groups. Establishment of a quantifiable relationship between G. vaginalis subgroup distribution and clinical status could have significant diagnostic implications.

A molecular enrichment strategy based on cpn60 for detection of epsilon-proteobacteria in the dog fecal microbiome

Members of the rare microbiome can be important components of complex microbial communities. For example, pet dog ownership is a known risk factor for human campylobacteriosis, and Campylobacter is commonly detected in dog feces by targeted assays. However, these organisms have not been detected by metagenomic methods. The goal of this study was to characterize fecal microbiota from healthy and diarrheic pet dogs using two different levels of molecular detection. PCR amplification and pyrosequencing of the universal cpn60 gene target was used to obtain microbial profiles from each dog. To investigate the relatively rare epsilon-proteobacteria component of the microbiome, a molecular enrichment was carried out using a PCR that first amplified the cpn10-cpn60 region from epsilon-proteobacteria, followed by universal cpn60 target amplification and pyrosequencing. From the non-enriched survey, the major finding was a significantly higher proportion of Bacteroidetes, notably Bacteroides vulgatus, in healthy dogs compared to diarrheic dogs. Epsilon-proteobacteria from the genera Helicobacter and Campylobacter were also detected at a low level in the non-enriched profiles of some dogs. Molecular enrichment increased the proportion of epsilon-proteobacteria sequences detected from each dog, as well as identified novel, presumably rare sequences not seen in the non-enriched profiles. Enriched profiles contained known species of Arcobacter, Campylobacter, Flexispira, and Helicobacter and identified two possibly novel species. These findings add to our understanding of the canine fecal microbiome in general, the epsilon-proteobacteria component specifically, and present a novel modification to traditional metagenomic approaches for study of the rare microbiome.