Two-Target Quantitative PCR To Predict Library Composition for Shallow Shotgun Sequencing

Gardnerella diversity and ecology in pregnancy and preterm birth

The vaginal microbiome has been linked to negative health outcomes including preterm birth. Specific taxa, including Gardnerella spp., have been identified as risk factors for these conditions. Historically, microbiome analysis methods have treated all Gardnerella spp. as one species, but the broad diversity of Gardnerella has become more apparent. We explore the diversity of Gardnerella clades and genomic species in the vaginal microbiome of pregnant women and their associations with microbiome composition and preterm birth. Relative abundance of Gardnerella clades and genomic species and other taxa was quantified in shotgun metagenomic sequencing data from three distinct cohorts of pregnant women. We also assessed the diversity and abundance of Gardnerella variants in 16S rRNA gene amplicon sequencing data from seven previously conducted studies in differing populations. Individual microbiomes often contained multiple Gardnerella variants, and the number of clades was associated with increased microbial load, or the ratio of non-human reads to human reads. Taxon co-occurrence patterns were largely consistent across Gardnerella clades and among cohorts. Some variants previously described as rare were prevalent in other cohorts, highlighting the importance of surveying a diverse set of populations to fully capture the diversity of Gardnerella. The diversity of Gardnerella both across populations and within individual vaginal microbiomes has long been unappreciated, as has been the intra-species diversity of many other members of the vaginal microbiome. The broad genomic diversity of Gardnerella has led to its reclassification as multiple species; here we demonstrate the diversity of Gardnerella found within and between vaginal microbiomes.IMPORTANCEThe present study shows that single microbiomes can contain all currently known species of Gardnerella and that multiple similar species can exist within the same environment. Furthermore, surveys of demographically distinct populations suggest that some species appear more commonly in certain populations. Further studies in broad and diverse populations will be necessary to fully understand the ecological roles of each Gardnerella sp., how they can co-exist, and their distinct impacts on microbial communities, preterm birth, and other health outcomes.

The airway microbiome of persons with cystic fibrosis correlates with acquisition and microbiological outcomes of incident Stenotrophomonas maltophilia infection

Rationale

Chronic infection with Stenotrophomonas maltophilia in persons with cystic fibrosis (pwCF) has been linked to an increased risk of pulmonary exacerbations and lung function decline. We sought to establish whether baseline sputum microbiome associates with risk of S. maltophilia incident infection and persistence in pwCF.

Methods

pwCF experiencing incident S. maltophilia infections attending the Calgary Adult CF Clinic from 2010-2018 were compared with S. maltophilia-negative sex, age (+/-2 years), and birth-cohort-matched controls. Infection outcomes were classified as persistent (when the pathogen was recovered in ≥50% of cultures in the subsequent year) or transient. We assessed microbial communities from prospectively biobanked sputum using V3-V4 16S ribosomal RNA (rRNA) gene sequencing, in the year preceding (Pre) (n = 57), at (At) (n = 22), and after (Post) (n = 31) incident infection. We verified relative abundance data using S. maltophilia-specific qPCR and 16S rRNA-targeted qPCR to assess bioburden. Strains were typed using pulse-field gel electrophoresis.

Results

Twenty-five pwCF with incident S. maltophilia (56% female, median 29 years, median FEV1 61%) with 33 total episodes were compared with 56 uninfected pwCF controls. Demographics and clinical characteristics were similar between cohorts. Among those with incident S. maltophilia infection, sputum communities did not cluster based on infection timeline (Pre, At, Post). Communities differed between the infection cohort and controls (n = 56) based on Shannon Diversity Index (SDI, p = 0.04) and clustered based on Aitchison distance (PERMANOVA, p = 0.01) prior to infection. At the time of incident S. maltophilia isolation, communities did not differ in SDI but clustered based on Aitchison distance (PERMANOVA, p = 0.03) in those that ultimately developed persistent infection versus those that were transient. S. maltophilia abundance within sputum was increased in samples from patients (Pre) relative to controls, measuring both relative (p = 0.004) and absolute (p = 0.001). Furthermore, S. maltophilia abundance was increased in sputum at incident infection in those who ultimately developed persistent infection relative to those with transient infection, measured relatively (p = 0.04) or absolute (p = 0.04), respectively.

Conclusion

Microbial community composition of CF sputum associates with S. maltophilia infection acquisition as well as infection outcome. Our study suggests sputum microbiome may serve as a surrogate for identifying infection risk and persistence risk.

Microbial DNA extraction of high-host content and low biomass samples: Optimized protocol for nasopharynx metagenomic studies

Introduction

Low microbial biomass and high human DNA content in nasopharyngeal aspirate samples hinder comprehensive characterization of microbiota and resistome. We obtained samples from premature infants, a group with increased risk of developing respiratory disorders and infections, and consequently frequent exposure to antibiotics. Our aim was to devise an optimal protocol for handling nasopharyngeal aspirate samples from premature infants, focusing on host DNA depletion and microbiome and resistome characterization.

Methods

Three depletion and three DNA extraction protocols were compared, using RT-PCR and whole metagenome sequencing to determine the efficiency of human DNA removal, taxonomic profiling and assignment of antibiotic resistance genes. Protocols were tested using mock communities, as well as pooled and individual patient samples.

Results

The only extraction protocol to retrieve the expected DNA yield from mock community samples was based on a lytic method to improve Gram positive recovery (MasterPure™). Host DNA content in non-depleted aliquots from pooled patient samples was 99%. Only samples depleted with MolYsis™ showed satisfactory, but varied reduction in host DNA content, in both pooled and individual patient samples, allowing for microbiome and resistome characterisation (host DNA content from 15% to 98%). Other depletion protocols either retrieved too low total DNA yields, preventing further analysis, or failed to reduce host DNA content. By using Mol_MasterPure protocol on aliquots from pooled patient samples, we increased the number of bacterial reads by 7.6 to 1,725.8-fold compared to non-depleted reference samples. PCR results were indicative of achieved microbial enrichment. Individual patient samples processed with Mol_MasterPure protocol varied greatly in total DNA yield, host DNA content (from 40% to 98%), species and antibiotic resistance gene richness.

Discussion

Despite high human DNA and low microbial biomass content in nasopharynx aspirates of preterm infants, we were able to reduce host DNA content to levels compatible with downstream shotgun metagenomic analysis, including bacterial species identification and coverage of antibiotic resistance genes. Whole metagenomic sequencing of microbes colonizing the nasopharynx may contribute to explaining the possible role of airway microbiota in respiratory conditions and reveal carriage of antibiotic resistance genes.

Performance Characteristics of Next-Generation Sequencing for the Detection of Antimicrobial Resistance Determinants in Escherichia coli Genomes and Metagenomes

Short-read sequencing can provide detection of multiple genomic determinants of antimicrobial resistance from single bacterial genomes and metagenomic samples. Despite its increasing application in human, animal, and environmental microbiology, including human clinical trials, the performance of short-read Illumina sequencing for antimicrobial resistance gene (ARG) detection, including resistance-conferring single nucleotide polymorphisms (SNPs), has not been systematically characterized. Using paired-end 2 × 150 bp (base pair) Illumina sequencing and an assembly-based method for ARG prediction, we determined sensitivity, positive predictive value (PPV), and sequencing depths required for ARG detection in an Escherichia coli isolate of sequence type (ST) 38 spiked into a synthetic microbial community at varying abundances. Approximately 300,000 reads or 15× genome coverage was sufficient to detect ARGs in E. coli ST38, with comparable sensitivity and PPV to ~100× genome coverage. Using metagenome assembly of mixed microbial communities, ARG detection at E. coli relative abundances of 1% would require assembly of approximately 30 million reads to achieve 15× target coverage. The minimum sequencing depths were validated using public data sets of 948 E. coli genomes and 10 metagenomic rectal swab samples. A read-based approach using k-mer alignment (KMA) for ARG prediction did not substantially improve minimum sequencing depths for ARG detection compared to assembly of the E. coli ST38 genome or the combined metagenomic samples. Analysis of sequencing depths from recent studies assessing ARG content in metagenomic samples demonstrated that sequencing depths had a median estimated detection frequency of 84% (interquartile range: 30%-92%) for a relative abundance of 1%. IMPORTANCE Systematically determining Illumina sequencing performance characteristics for detection of ARGs in metagenomic samples is essential to inform study design and appraisal of human, animal, and environmental metagenomic antimicrobial resistance studies. In this study, we quantified the performance characteristics of ARG detection in E. coli genomes and metagenomes and established a benchmark of ~15× coverage for ARG detection for E. coli in metagenomes. We demonstrate that for low relative abundances, sequencing depths of ~30 million reads or more may be required for adequate sensitivity for many applications.