Resolution of phenotypically distinct strains of Enterococcus spp. in a complex microbial community using cpn60 universal target sequencing

isolateR: an R package for generating microbial libraries from Sanger sequencing data

MOTIVATION

Sanger sequencing of taxonomic marker genes (e.g. 16S/18S/ITS/rpoB/cpn60) represents the leading method for identifying a wide range of microorganisms including bacteria, archaea, and fungi. However, the manual processing of sequence data and limitations associated with conventional BLAST searches impede the efficient generation of strain libraries essential for cataloging microbial diversity and discovering novel species.

RESULTS

isolateR addresses these challenges by implementing a standardized and scalable three-step pipeline that includes: (1) automated batch processing of Sanger sequence files, (2) taxonomic classification via global alignment to type strain databases in accordance with the latest international nomenclature standards, and (3) straightforward creation of strain libraries and handling of clonal isolates, with the ability to set customizable sequence dereplication thresholds and combine data from multiple sequencing runs into a single library. The tool's user-friendly design also features interactive HTML outputs that simplify data exploration and analysis. Additionally, in silico benchmarking done on two comprehensive human gut genome catalogues (IMGG and Hadza hunter-gather populations) showcase the proficiency of isolateR in uncovering and cataloging the nuanced spectrum of microbial diversity, advocating for a more targeted and granular exploration within individual hosts to achieve the highest strain-level resolution possible when generating culture collections.

AVAILABILITY AND IMPLEMENTATION

isolateR is available at: https://github.com/bdaisley/isolateR.

Update on cpnDB: a reference database of chaperonin sequences

cpnDB was established in 2004 to provide a manually curated database of type I (60 kDa chaperonin, CPN60, also known as GroEL or HSP60) and type II (CCT, TRiC, thermosome) chaperonin sequences and to support chaperonin sequence-based applications including microbial species identification, detection and quantification, phylogenetic investigations and microbial community profiling. Since its establishment, cpnDB has grown to over 25 000 sequence records including over 4 000 records from bacterial type strains. The updated cpnDB webpage (www.cpndb.ca) provides tools for text- or sequence-based searches and links to protocols, and selected reference data sets are available for download. Here we present an updated description of the contents and taxonomic coverage of cpnDB and an analysis of cpn60 sequence diversity.

Identification of Enterococci, Staphylococci, and Enterobacteriaceae from Slurries and Air in and around Two Pork Farms

In this study, we investigated the airborne dissemination of bacteria from the inside of two very different pork farms (an intensively confined farm and an open-range farm) to the immediate environment. Samples were taken from the slurry, from the air inside the farms (area 0), and from their immediate surroundings at a distance of 50, 100, and 150 m in four directions (north, south, east, and west). A control sample in the air of a zone far away from human or animal activity was also taken. Identification of isolates was made by means of the matrix-assisted laser desorption-ionization time of flight system. A total of 1,063 isolates were obtained, of which a mere 7 came from the air of the control area. Staphylococci, enterococci, and Enterobacteriaceae were selectively targeted for isolation and represented 48.6, 27.2, and 21.6% of the isolates, respectively. The species identified from the air of surrounding areas ( Enterococcus faecalis, Enterococcus hirae, and Staphylococcus arlettae, mainly) were also present inside the farms studied. The results suggest that air is involved in bacterial dissemination, and pork farms should be considered a potential source of foodborne bacteria that might contaminate surrounding areas, including vegetable orchards. Wind direction appears as a factor involved in bacterial dispersion through the air, but its effect may be conditioned by existing vegetation and orographic conditions.

The vaginal microbiome of pregnant women is less rich and diverse, with lower prevalence of Mollicutes, compared to non-pregnant women

The vaginal microbiome plays an important role in maternal and neonatal health. Imbalances in this microbiota (dysbiosis) during pregnancy are associated with negative reproductive outcomes, such as pregnancy loss and preterm birth, but the underlying mechanisms remain poorly understood. Consequently a comprehensive understanding of the baseline microbiome in healthy pregnancy is needed. We characterized the vaginal microbiomes of healthy pregnant women at 11–16 weeks of gestational age (n = 182) and compared them to those of non-pregnant women (n = 310). Profiles were created by pyrosequencing of the cpn60 universal target region. Microbiome profiles of pregnant women clustered into six Community State Types: I, II, III, IVC, IVD and V. Overall microbiome profiles could not be distinguished based on pregnancy status. However, the vaginal microbiomes of women with healthy ongoing pregnancies had lower richness and diversity, lower prevalence of Mycoplasma and Ureaplasma and higher bacterial load when compared to non-pregnant women. Lactobacillus abundance was also greater in the microbiomes of pregnant women with Lactobacillus-dominated CSTs in comparison with non-pregnant women. This study provides further information regarding characteristics of the vaginal microbiome of low-risk pregnant women, providing a baseline for forthcoming studies investigating the diagnostic potential of the microbiome for prediction of adverse pregnancy outcomes.

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 Study of the Infant Nasal Microbiome Development over the First Year of Life and in Relation to Their Primary Adult Caregivers Using cpn60 Universal Target (UT) as a Phylogenetic Marker

Whereas the infant gut microbiome is the subject of intense study, relatively little is known regarding the nares microbiome in newborns and during early life. This study aimed to survey the typical composition and diversity of human anterior nare microflora for developing infants over time, and to explore how these correlate to their primary caregivers. Single nare swabs were collected at five time points over a one-year period for each subject from infant-caregiver pairs. Our study comprised of 50 infants (recruited at 2 weeks, post delivery) and their 50 primary caregivers. Applying the chaperonin-60 (cpn60) universal target (UT) amplicon as our molecular barcoding marker to census survey the microbial communities, we longitudinally surveyed infant nares microbiota at 5 time points over the course of the first year of life. The inter- and intra-subject diversity was catalogued and compared, both longitudinally and relative to their adult primary caregivers. Although within-subject variability over time and inter-subject variability were both observed, the assessment detected only one or two predominant genera for individual infant samples, belonging mainly to phyla Actinobacteria, Firmicutes, and Proteobacteria. Consistent with previously observed microbial population dynamics in other body sites, the diversity of nares microflora increased over the first year of life and infants showed differential operational taxonomic units (OTUs) relative to their matched primary caregiver. The collected evidence also support that both temporal and seasonal changes occur with respect to carriage of potentially pathogenic bacteria (PPBs), which may influence host predisposition to infection. This pilot study surveying paired infant/caregiver nare microbiomes provides novel longitudinal diversity information that is pertinent to better understanding nare microbiome development in infants.

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.

Use of MALDI-TOF mass spectrometry fingerprinting to characterize Enterococcus spp. and Escherichia coli isolates

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is a faster and more accurate method to identify intact bacteria than conventional microbiology and/or molecular biology methods. The MALDI-TOF MS method is potentially applicable in diagnostic laboratories to characterize commensal bacterial species, some of which are major pathogens, from human or animal gastrointestinal tracts. The aim of this study was to analyze at the cluster and statistical level the capacity of MALDI-TOF MS to distinguish between previously characterized enterococci and Escherichia coli isolated from wild birds of the Azores archipelago. Soluble proteins were extracted from intact cell cultures of 60 isolates of Enterococcus spp. and 60 isolates of E. coli by an expedient method. MALDI-TOF MS was used to obtain 1200 mass spectra that were statistically analyzed and compared. A total of 215 distinct mass-to-charge (m/z) peaks were obtained, including a peak at m/z 4428±3, which is exclusively found in spectra from Enterococcus isolates, and peaks at m/z 5379±3 and m/z 6253±3, which are only detected in spectra from E. coli isolates. By processing mass spectra and analyzing them statistically with MassUp software, including principal component analysis (PCA) and clustering, it was possible to correctly distinguish between isolates of Enterococcus and Escherichia genera. The results of the proteomic analysis confirm that these tools could be used to characterize whole bacterial cells. In the future, with an optimized protocol for obtaining plasmid-associated proteins and the further development of bioinformatics methods, it is likely that mass peak characteristic of antimicrobial resistance will be detected, increasing the potential usefulness of MALDI-TOF in routine clinical assays.

BIOLOGICAL SIGNIFICANCE

This study highlights the importance of MALDI-TOF MS in the rapid and reliable identification of bacteria by whole-cell analysis. The mass spectrometry approach performed in this study further contributes for the microbial biomarker discovery culminating in a preferable bacteria identification and antimicrobial resistance tool for the future of clinical microbiology. This article is part of a Special Issue entitled: HUPO 2014.

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.

CRISPRs of Enterococcus faecalis and E. hirae isolates from pig feces have species-specific repeats but share some common spacer sequences

Clustered regularly interspaced short palindromic repeats (CRISPR) are currently a topic of interest in microbiology due to their role as a prokaryotic immune system. Investigations of CRISPR distribution and characterization to date have focused on pathogenic bacteria, while less is known about CRISPR in commensal bacteria, where they may have a significant role in the ecology of the microbiota of humans and other animals, and act as a recorder of interactions between bacteria and viruses. A combination of PCR and sequencing was used to determine prevalence and distribution of CRISPR arrays in Enterococcus faecalis and Enterococcus hirae isolates from the feces of healthy pigs. Both type II CRISPR-Cas and Orphan CRISPR (without Cas genes) were detected in the 195 isolates examined. CRISPR-Cas was detected in 52 (46/88) and 42 % (45/107) E. faecalis and E. hirae isolates, respectively. The prevalence of Orphan CRISPR arrays was higher in E. faecalis isolates (95 %, 84/88) compared with E. hirae isolates (49 %, 53/107). Species-specific repeat sequences were identified in Orphan CRISPR arrays, and 42 unique spacer sequences were identified. Only two spacers matched previously characterized pig virome sequences, and many were apparently derived from chromosomal sequences of enterococci. Surprisingly, 17 (40 %) of the spacers were detected in both species. Shared spacer sequences are evidence of a lack of species specificity in the agents and mechanisms responsible for integration of spacers, and the abundance of spacer sequences corresponding to bacterial chromosomal sequences reflects interspecific interactions within the intestinal microbiota.

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.

Development of quantitative PCR assays targeting the 16S rRNA genes of Enterococcus spp. and their application to the identification of enterococcus species in environmental samples

The detection of environmental enterococci has been determined primarily by using culture-based techniques that might exclude some enterococcal species as well as those that are nonculturable. To address this, the relative abundances of enterococci were examined by challenging fecal and water samples against a currently available genus-specific assay (Entero1). To determine the diversity of enterococcal species, 16S rRNA gene-based group-specific quantitative PCR (qPCR) assays were developed and evaluated against eight of the most common environmental enterococcal species. Partial 16S rRNA gene sequences of 439 presumptive environmental enterococcal strains were analyzed to study further the diversity of enterococci and to confirm the specificities of group-specific assays. The group-specific qPCR assays showed relatively high amplification rates with targeted species (>98%), although some assays cross-amplified with nontargeted species (1.3 to 6.5%). The results with the group-specific assays also showed that different enterococcal species co-occurred in most fecal samples. The most abundant enterococci in water and fecal samples were Enterococcus faecalis and Enterococcus faecium, although we identified more water isolates as Enterococcus casseliflavus than as any of the other species. The prevalence of the Entero1 marker was in agreement with the combined number of positive signals determined by the group-specific assays in most fecal samples, except in gull feces. On the other hand, the number of group-specific assay signals was lower in all water samples tested, suggesting that other enterococcal species are present in these samples. While the results highlight the value of genus- and group-specific assays for detecting the major enterococcal groups in environmental water samples, additional studies are needed to determine further the diversity, distributions, and relative abundances of all enterococcal species found in water.

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.

Development and validation of a real-time TaqMan assay for the detection and enumeration of Pseudomonas fluorescens ATCC 13525 used as a challenge organism in testing of food equipments

Pseudomonas fluorescens ATCC 13525 is used as the challenge organism to evaluate the efficacy of the clean-in-place (CIP) process of food equipment (automatic ice-maker) as per NSF/ANSI Standard 12. Traditional culturing methodology is presently used to determine the concentration of the challenge organism, which takes 48 h to confirm the cell density. Storage of the challenge preparation in the refrigerator might alter the cell density as P. fluorescens is capable of growing at 4 °C. Also, background organism can grow on the Pseudomonas F agar (PFA) used for the recovery of P. fluorescens thus affecting the results of the test. Real-time TaqMan assay targeting the cpn60 gene was developed for the enumeration and the identification of P. fluorescens because of its specificity, accuracy, and shorter turnaround time. The TaqMan primer-probe pair developed using the Allele ID® 7.0 probe design software was highly specific and sensitive for the target organism. The sensitivity of the assay was 10 colony forming units (CFU)/mL. The assay was also successful in determining the concentration of the challenge preparation within 2 h. Based on these observations, TaqMan assay targeting the cpn60 gene can be efficiently used for strain level identification and enumeration of bacteria.

PRACTICAL APPLICATION

Pseudomonas fluorescens ATCC 13525 is used as a challenge organism in the efficacy testing of clean-in-place process of food equipments. Currently, culturing technique is used for its identification and estimation, which is not only time-consuming but also prone to error. Real-time TaqMan assay is more specific, sensitive, and accurate along with a shorter turnaround time compared to culturing techniques, thereby increasing the overall quality of the testing methodology to evaluate the clean-in-place process critical for the food industry to protect public health and safety.

Isolates of Thermoanaerobacter thermohydrosulfuricus from decaying wood compost display genetic and phenotypic microdiversity

In this study, 12 strains of Thermoanaerobacter were isolated from a single decaying wood compost sample and subjected to genetic and phenotypic profiling. The 16S rRNA encoding gene sequences suggested that the isolates were most similar to strains of either Thermoanaerobacter pseudethanolicus or Thermoanaerobacter thermohydrosulfuricus. Examination of the lesser conserved chaperonin-60 (cpn60) universal target showed that some isolates shared the highest sequence identity with T. thermohydrosulfuricus; however, others to Thermoanaerobacter wiegelii and Thermoanaerobacter sp. Rt8.G4 (formerly Thermoanaerobacter brockii Rt8.G4). BOX-PCR fingerprinting profiles identified differences in the banding patterns not only between the isolates and the reference strains, but also among the isolates themselves. To evaluate the extent these genetic differences were manifested phenotypically, the utilization patterns of 30 carbon substrates were examined and the niche overlap indices (NOI) calculated. Despite showing a high NOI (> 0.9), significant differences existed in the substrate utilization capabilities of the isolates suggesting that either a high degree of niche specialization or mechanisms allowing for non-competitive co-existence, were present within this ecological context. Growth studies showed that the isolates were physiologically distinct in both growth rate and the fermentation product ratios. Our data indicate that phenotypic diversity exists within genetically microdiverse Thermoanaerobacter isolates from a common environment.

A 'universal' type II chaperonin PCR detection system for the investigation of Archaea in complex microbial communities

Bacteria and Archaea are evolutionarily and biochemically distinct domains found together in many environments. Robust 'universal' PCR primer sets targeting both the bacterial 16S rRNA gene and the type I chaperonin gene have been established. However, 'universal' PCR primers for Archaea are currently limited to the 16S rRNA gene. We investigated the type II chaperonin (known as the thermosome, TF55, CCT or TCP-1) as a potential universal target (UT) for Archaea. Reproducible amplification of thermosome gene sequences from all major phyla tested was achieved through the application of a mixture or 'cocktail' of two forward and two reverse primers. Phylogenies based on the ∼750-bp thermosome UT were congruent with 16S rRNA gene phylogenies while exhibiting longer branch lengths, improving resolution of closely related taxa. 'Universal' thermosome primers were applied to profiling the archaeal community of dairy cow rumen and results compared with profiles based on the 16S rRNA gene and methyl co-enzyme M reductase (methanogen-specific) gene. Clone libraries generated from each target gene, as well as a pyrosequencing profile of one thermosome rumen library, revealed that all three targets consistently detected Methanobrevibacter smithii, Methanobrevibacter ruminantium and Methanosphaera stadtmanae as the dominant constituents; however, thermosome gene sequences were more diverse than either of the other targets providing a higher resolution description of the archaeal community. These findings demonstrate that a 'universal' thermosome PCR protocol is a powerful metagenomic tool for detecting and characterizing Archaea and archaeal communities.