Evaluation of variant calling for cpn60 barcode sequence-based microbiome profiling

Birth Mode Does Not Determine the Presence of Shared Bacterial Strains between the Maternal Vaginal Microbiome and the Infant Stool Microbiome

Dysbiosis of the neonatal gut microbiome during early life has been suggested as the missing link that may explain higher rates of certain diseases in caesarean section-delivered infants. Many studies report delivery mode-related dysbiosis in infants due to a lack of maternal vaginal microbiome exposure, prompting interventions to correct the neonatal gut microbiome by transferring these missing microbes after caesarean delivery. The maternal vaginal microbiome is among the first microbial exposures that many infants experience, yet little is known about the extent of direct transmission of maternal vaginal microbes. As part of the Maternal Microbiome Legacy Project, we aimed to determine if maternal vaginal bacteria are vertically transmitted to infants. We employed cpn60 microbiome profiling, culture-based screening, molecular strain typing, and whole-genome sequencing to determine whether identical maternal vaginal strains were present in infant stool microbiomes. We identified identical cpn60 sequence variants in both halves of maternal-infant dyads in 204 of 585 Canadian women and their newborn infants (38.9%). The same species of Bifidobacterium and Enterococcus were cultured from maternal and corresponding infant samples in 33 and 13 of these mother-infant dyads, respectively. Pulsed-field gel electrophoresis and whole-genome sequencing determined that near-identical strains were detected in these dyads irrespective of delivery mode, indicating an alternative source in cases of caesarean delivery. Overall, we demonstrated that vertical transmission of maternal vaginal microbiota is likely limited and that transmission from other maternal body sites, such as the gut and breast milk, may compensate for the lack of maternal vaginal microbiome exposure during caesarean delivery. IMPORTANCE The importance of the gut microbiome in human health and disease is widely recognized, and there has been a growing appreciation that alterations in gut microbiome composition during a "critical window" of development may impact health in later life. Attempts to correct gut microbiome dysbiosis related to birth mode are underpinned by the assumption that the lack of exposure to maternal vaginal microbes during caesarean delivery is responsible for dysbiosis. Here, we demonstrate that there is limited transmission of the maternal vaginal microbiome to the neonatal gut, even in cases of vaginal delivery. Furthermore, the presence of identical strains shared between mothers and infants in early life, even in cases of caesarean delivery, highlights compensatory microbial exposures and sources for the neonatal stool microbiome other than the maternal vagina.

Rapid and accurate taxonomic classification of cpn60 amplicon sequence variants

The "universal target" region of the gene encoding the 60 kDa chaperonin protein (cpn60, also known as groEL or hsp60) is a proven sequence barcode for bacteria and a useful target for marker gene amplicon-based studies of complex microbial communities. To date, identification of cpn60 sequence variants from microbiome studies has been accomplished by alignment of queries to a reference database. Naïve Bayesian classifiers offer an alternative identification method that provides variable rank classification and shorter analysis times. We curated a set of cpn60 barcode sequences to train the RDP classifier and tested its performance on data from previous human microbiome studies. Results showed that sequences accounting for 79%, 86% and 92% of the observations (read counts) in saliva, vagina and infant stool microbiome data sets were classified to the species rank. We also trained the QIIME 2 q2-feature-classifier on cpn60 sequence data and demonstrated that it gives results consistent with the standalone RDP classifier. Successful implementation of a naïve Bayesian classifier for cpn60 sequences will facilitate future microbiome studies and open opportunities to integrate cpn60 amplicon sequence identification into existing analysis pipelines.

Evaluating cpn60 for high-resolution profiling of the mammalian skin microbiome and detection of phylosymbiosis

Despite being the most widely used phylogenetic marker for amplicon-based profiling of microbial communities, limited phylogenetic resolution of the 16S rRNA gene limits its use for studies of host-microbe co-evolution. In contrast, the cpn60 gene is a universal phylogenetic marker with greater sequence variation capable of species-level resolution. This research compared mammalian skin microbial profiles generated from cpn60 and 16S rRNA gene sequencing approaches, testing for patterns of phylosymbiosis that suggest co-evolutionary host-microbe associations. An ~560 bp fragment of the cpn60 gene was amplified with universal primers and subjected to high-throughput sequencing. Taxonomic classification of cpn60 sequences was completed using a naïve-Bayesian QIIME2 classifier created for this project, trained with an NCBI-supplemented curated cpn60 database (cpnDB_nr). The cpn60 dataset was then compared to published 16S rRNA gene amplicon data. Beta diversity comparisons of microbial community profiles generated with cpn60 and 16S rRNA gene amplicons were not significantly different, based on Procrustes analysis of Bray-Curtis and UniFrac distances. Despite similar relationships among skin microbial profiles, improved phylogenetic resolution provided by the cpn60 gene sequencing permitted observations of phylosymbiosis between microbial community profiles and their mammalian hosts that were not previously observed with 16S rRNA gene profiles. Subsequent investigation of Staphylococcaceae taxa using the cpn60 gene showed increased phylogenetic resolution compared the 16S rRNA gene profiles, revealing potential co-evolutionary host-microbe associations. Overall, our results demonstrate that 16S rRNA and cpn60 marker genes generate comparable microbial community composition patterns while cpn60 better facilitates analyses, such as phylosymbiosis, that require increased phylogenetic resolution.

The effect of wheat genotype on the microbiome is more evident in roots and varies through time

Crop breeding has traditionally ignored the plant-associated microbial communities. Considering the interactions between plant genotype and associated microbiota is of value since different genotypes of the same crop often harbor distinct microbial communities which can influence the plant phenotype. However, recent studies have reported contrasting results, which led us to hypothesize that the effect of genotype is constrained by growth stages, sampling year and plant compartment. To test this hypothesis, we sampled bulk soil, rhizosphere soil and roots of 10 field-grown wheat genotypes, twice per year, for 4 years. DNA was extracted and regions of the bacterial 16 S rRNA and CPN60 genes and the fungal ITS region were amplified and sequenced. The effect of genotype was highly contingent on the time of sampling and on the plant compartment sampled. Only for a few sampling dates, were the microbial communities significantly different across genotypes. The effect of genotype was most often significant for root microbial communities. The three marker genes used provided a highly coherent picture of the effect of genotype. Taken together, our results confirm that microbial communities in the plant environment strongly vary across compartments, growth stages, and years, and that this can mask the effect of genotype.

Maternal vaginal microbiome composition does not affect development of the infant gut microbiome in early life

Birth mode has been implicated as a major factor influencing neonatal gut microbiome development, and it has been assumed that lack of exposure to the maternal vaginal microbiome is responsible for gut dysbiosis among caesarean-delivered infants. Consequently, practices to correct dysbiotic gut microbiomes, such as vaginal seeding, have arisen while the effect of the maternal vaginal microbiome on that of the infant gut remains unknown. We conducted a longitudinal, prospective cohort study of 621 Canadian pregnant women and their newborn infants and collected pre-delivery maternal vaginal swabs and infant stool samples at 10-days and 3-months of life. Using cpn60-based amplicon sequencing, we defined vaginal and stool microbiome profiles and evaluated the effect of maternal vaginal microbiome composition and various clinical variables on the development of the infant stool microbiome. Infant stool microbiomes showed significant differences in composition by delivery mode at 10-days postpartum; however, this effect could not be explained by maternal vaginal microbiome composition and was vastly reduced by 3 months. Vaginal microbiome clusters were distributed across infant stool clusters in proportion to their frequency in the overall maternal population, indicating independence of the two communities. Intrapartum antibiotic administration was identified as a confounder of infant stool microbiome differences and was associated with lower abundances of Escherichia coli, Bacteroides vulgatus, Bifidobacterium longum and Parabacteroides distasonis. Our findings demonstrate that maternal vaginal microbiome composition at delivery does not affect infant stool microbiome composition and development, suggesting that practices to amend infant stool microbiome composition focus factors other than maternal vaginal microbes.

Gardnerella subgroup dominant microbiomes are associated with divergent cervicovaginal immune responses in a longitudinal cohort of Kenyan women

Most cervicovaginal microbiome-immunology studies to date have relied on 16S rDNA microbial profiling which does not resolve the molecular subgroups of Gardnerella, believed to be central to the pathogenesis of bacterial vaginosis (BV) and subsequent risk of HIV acquisition. Here we used the cpn60 universal target which in addition to other microbial taxa, resolves four Gardnerella subgroups, for cervicovaginal microbial profiling in a longitudinal cohort of Kenyan women to examine associations with cellular and soluble markers of inflammation and HIV susceptibility. Participants (N = 41) were sampled, contributing 362 samples for microbiome analysis. All non-Lactobacillus dominant microbial communities were associated with high pro-inflammatory cytokine levels. Divergent associations were observed among different Gardnerella subgroup dominated communities with respect to the chemokine IP-10. Specifically, Gardnerella subgroup A dominant and polymicrobial communities were associated with reduced concentrations of IP-10 in adjusted linear mixed models (p<0.0001), compared to microbial communities dominated by Lactobacillus (non-iners) species. However, these associations did not translate to significant differences in the proportion or absolute number of CCR5, HLA-DR and CD38 expressed on cervical CD4⁺ T- cells. These findings suggest that some associations between Gardnerella subgroup dominant microbiomes and mucosal immunity differ and are relevant for the study of BV-pathogenesis and understanding the mechanisms of BV-associated HIV risk.

DNA barcoding, an effective tool for species identification: a review

DNA barcoding is a powerful taxonomic tool to identify and discover species. DNA barcoding utilizes one or more standardized short DNA regions for taxon identification. With the emergence of new sequencing techniques, such as Next-generation sequencing (NGS), ONT MinION nanopore sequencing, and Pac Bio sequencing, DNA barcoding has become more accurate, fast, and reliable. Rapid species identification by DNA barcodes has been used in a variety of fields, including forensic science, control of the food supply chain, and disease understanding. The Consortium for Barcode of Life (CBOL) presents various working groups to identify the universal barcode gene, such as COI in metazoans; rbcL, matK, and ITS in plants; ITS in fungi; 16S rRNA gene in bacteria and archaea, and creating a reference DNA barcode library. In this article, an attempt has been made to analyze the various proposed DNA barcode for different organisms, strengths & limitations, recent advancements in DNA barcoding, and methods to speed up the DNA barcode reference library construction. This study concludes that constructing a reference library with high species coverage would be a major step toward identifying species by DNA barcodes. This can be achieved in a short period of time by using advanced sequencing and data analysis methods.

Life barcoded by DNA barcodes

The modern concept of DNA-based barcoding for cataloguing biodiversity was proposed in 2003 by first adopting an approximately 600 bp fragment of the mitochondrial COI gene to compare via nucleotide alignments with known sequences from specimens previously identified by taxonomists. Other standardized regions meeting barcoding criteria then are also evolving as DNA barcodes for fast, reliable and inexpensive assessment of species composition across all forms of life, including animals, plants, fungi, bacteria and other microorganisms. Consequently, global DNA barcoding campaigns have resulted in the formation of many online workbenches and databases, such as BOLD system, as barcode references, and facilitated the development of mini-barcodes and metabarcoding strategies as important extensions of barcode techniques. Here we intend to give an overview of the characteristics and features of these barcode markers and major reference libraries existing for barcoding the planet’s life, as well as to address the limitations and opportunities of DNA barcodes to an increasingly broader community of science and society.

Vaginal microbiota associated with oncogenic HPV in a cohort of HPV-vaccinated women living with HIV

Background

Women living with HIV (WLWH) experience higher rates of human papillomavirus (HPV) infection and cervical cancer than women without HIV. Changes in the vaginal microbiome have been implicated in HPV-related disease processes such as persistence of high-risk HPV infection but this has not been well defined in a population living with HIV.

Methods

Four hundred and 20 girls and WLWH, age ≥9, across 14 clinical sites in Canada were enrolled to receive three doses of quadrivalent HPV vaccine for assessment of vaccine immunogenicity. Blood, cervical cytology, and cervico-vaginal swabs were collected. Cervico-vaginal samples were tested for HPV DNA and underwent microbiota sequencing.

Results

Principal component analysis (PCA) and hierarchical clustering generated community state types (CSTs). Relationships between taxa and CSTs with HPV infection were examined using mixed-effects logistic regressions, Poisson regressions, or generalized linear mixed-effects models, as appropriate. Three hundred and fifty-six cervico-vaginal microbiota samples from 172 women were sequenced. Human papillomavirus DNA was detected in 211 (59%) samples; 110 (31%) contained oncogenic HPV. Sixty-five samples (18%) were taken concurrently with incident oncogenic HPV infection and 56 (16%) were collected from women with concurrent persistent oncogenic HPV infection.

Conclusions

No significant associations between taxa, CST, or microbial diversity and HPV-related outcomes were found. However, we observed weak associations between a dysbiotic microbiome and specific species, including Gardnerella, Porphyromonas, and Prevotella species, with incident HPV infection.

cpn60 barcode sequences accurately identify newly defined genera within the Lactobacillaceae

The cpn60 barcode sequence is established as an informative target for microbial species identification. Applications of cpn60 barcode sequencing are supported by the availability of "universal" PCR primers for its amplification and a curated reference database of cpn60 sequences, cpnDB. A recent reclassification of lactobacilli involving the definition of 23 new genera provided an opportunity to update cpnDB and to determine if the cpn60 barcode could be used for accurate identification of species consistent with the new framework. Analysis of 275 cpn60 sequences representing 258/269 of the validly named species in Lactobacillus, Paralactobacillus and the 23 newer genera showed that cpn60-based sequence relationships were generally consistent with the whole-genome-based phylogeny. Aligning or mapping full length barcode sequences or a 150 bp subsequence resulted in accurate and unambiguous species identification in almost all cases. Taken together, our results show that the combination of available reference sequence data, "universal" barcode amplification primers, and the inherent sequence diversity within the cpn60 barcode make it a useful target for the detection and identification of lactobacilli as defined by the latest taxonomic framework.

Early Neonatal Meconium Does Not Have a Demonstrable Microbiota Determined through Use of Robust Negative Controls with cpn60-Based Microbiome Profiling

Detection of bacterial DNA within meconium is often cited as evidence supporting in utero colonization. However, many studies fail to adequately control for contamination. We aimed to define the microbial content of meconium under properly controlled conditions. DNA was extracted from 141 meconium samples and subjected to cpn60-based microbiome profiling, with controls to assess contamination throughout. Total bacterial loads of neonatal meconium, infant stool, and controls were compared by 16S rRNA quantitative PCR (qPCR). Viable bacteria within meconium were cultured, and isolate clonality was assessed by pulsed-field gel electrophoresis (PFGE). Meconium samples did not differ significantly from controls with respect to read numbers or taxonomic composition. Twenty (14%) outliers with markedly higher read numbers were collected significantly later after birth and appeared more like transitional stool than meconium. Total bacterial loads were significantly higher in stool than in meconium, which did not differ from that of sequencing controls, and correlated well with read numbers. Cultured isolates were most frequently identified as Staphylococcus epidermidis, Enterococcus faecalis, or Escherichia coli, with PFGE indicating high intraspecies diversity. Our findings highlight the importance of robust controls in studies of low microbial biomass samples and argue against meaningful bacterial colonization in utero. Given that meconium microbiome profiles could not be distinguished from sequencing controls, and that viable bacteria within meconium appeared uncommon and largely consistent with postnatal skin colonization, there does not appear to be a meconium microbiota.

IMPORTANCE Much like the recent placental microbiome controversy, studies of neonatal meconium reporting bacterial communities within the fetal and neonatal gut imply that microbial colonization begins prior to birth. However, recent work has shown that placental microbiomes almost exclusively represent contamination from lab reagents and the environment. Here, we demonstrate that prior studies of neonatal meconium are impacted by the same issue, showing that the microbial content of meconium does not differ from negative controls that have never contained any biological material. Our culture findings similarly supported this notion and largely comprised bacteria normally associated with healthy skin. Overall, our work adds to the growing body of evidence against the in utero colonization hypothesis.

CaptureSeq: Hybridization-Based Enrichment of cpn60 Gene Fragments Reveals the Community Structures of Synthetic and Natural Microbial Ecosystems

Background. The molecular profiling of complex microbial communities has become the basis for examining the relationship between the microbiome composition, structure and metabolic functions of those communities. Microbial community structure can be partially assessed with “universal” PCR targeting taxonomic or functional gene markers. Increasingly, shotgun metagenomic DNA sequencing is providing more quantitative insight into microbiomes. However, both amplicon-based and shotgun sequencing approaches have shortcomings that limit the ability to study microbiome dynamics. Methods. We present a novel, amplicon-free, hybridization-based method (CaptureSeq) for profiling complex microbial communities using probes based on the chaperonin-60 gene. Molecular profiles of a commercially available synthetic microbial community standard were compared using CaptureSeq, whole metagenome sequencing, and 16S universal target amplification. Profiles were also generated for natural ecosystems including antibiotic-amended soils, manure storage tanks, and an agricultural reservoir. Results. The CaptureSeq method generated a microbial profile that encompassed all of the bacteria and eukaryotes in the panel with greater reproducibility and more accurate representation of high G/C content microorganisms compared to 16S amplification. In the natural ecosystems, CaptureSeq provided a much greater depth of coverage and sensitivity of detection compared to shotgun sequencing without prior selection. The resulting community profiles provided quantitatively reliable information about all three domains of life (Bacteria, Archaea, and Eukarya) in the different ecosystems. The applications of CaptureSeq will facilitate accurate studies of host-microbiome interactions for environmental, crop, animal and human health. Conclusions: cpn60-based hybridization enriched for taxonomically informative DNA sequences from complex mixtures. In synthetic and natural microbial ecosystems, CaptureSeq provided sequences from prokaryotes and eukaryotes simultaneously, with quantitatively reliable read abundances. CaptureSeq provides an alternative to PCR amplification of taxonomic markers with deep community coverage while minimizing amplification biases.