A Metagenomic Approach to Characterization of the Vaginal Microbiome Signature in Pregnancy

While current major national research efforts (i.e., the NIH Human Microbiome Project) will enable comprehensive metagenomic characterization of the adult human microbiota, how and when these diverse microbial communities take up residence in the host and during reproductive life are unexplored at a population level. Because microbial abundance and diversity might differ in pregnancy, we sought to generate comparative metagenomic signatures across gestational age strata. DNA was isolated from the vagina (introitus, posterior fornix, midvagina) and the V5V3 region of bacterial 16S rRNA genes were sequenced (454FLX Titanium platform). Sixty-eight samples from 24 healthy gravidae (18 to 40 confirmed weeks) were compared with 301 non-pregnant controls (60 subjects). Generated sequence data were quality filtered, taxonomically binned, normalized, and organized by phylogeny and into operational taxonomic units (OTU); principal coordinates analysis (PCoA) of the resultant beta diversity measures were used for visualization and analysis in association with sample clinical metadata. Altogether, 1.4 gigabytes of data containing >2.5 million reads (averaging 6,837 sequences/sample of 493 nt in length) were generated for computational analyses. Although gravidae were not excluded by virtue of a posterior fornix pH >4.5 at the time of screening, unique vaginal microbiome signature encompassing several specific OTUs and higher-level clades was nevertheless observed and confirmed using a combination of phylogenetic, non-phylogenetic, supervised, and unsupervised approaches. Both overall diversity and richness were reduced in pregnancy, with dominance of Lactobacillus species (L. iners crispatus, jensenii and johnsonii, and the orders Lactobacillales (and Lactobacillaceae family), Clostridiales, Bacteroidales, and Actinomycetales. This intergroup comparison using rigorous standardized sampling protocols and analytical methodologies provides robust initial evidence that the vaginal microbial 16S rRNA gene catalogue uniquely differs in pregnancy, with variance of taxa across vaginal subsite and gestational age.

[1]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[2]  R. Knight,et al.  Quantitative and Qualitative β Diversity Measures Lead to Different Insights into Factors That Structure Microbial Communities , 2007, Applied and Environmental Microbiology.

[3]  Hongzhe Li,et al.  Disordered Microbial Communities in the Upper Respiratory Tract of Cigarette Smokers , 2010, PloS one.

[4]  R. Knight,et al.  Rapid denoising of pyrosequencing amplicon data: exploiting the rank-abundance distribution , 2010, Nature Methods.

[5]  William A. Walters,et al.  QIIME allows analysis of high-throughput community sequencing data , 2010, Nature Methods.

[6]  R. Knight,et al.  Worlds within worlds: evolution of the vertebrate gut microbiota , 2008, Nature Reviews Microbiology.

[7]  Haifeng Lu,et al.  Symbiotic gut microbes modulate human metabolic phenotypes , 2008, Proceedings of the National Academy of Sciences.

[8]  R. Knight,et al.  Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. , 2009, Genome research.

[9]  D. Gavin,et al.  A statistical approach to evaluating distance metrics and analog assignments for pollen records , 2003, Quaternary Research.

[10]  D. Gavin,et al.  Pollen‐vegetation calibration for tundra communities in the Arctic Foothills, northern Alaska , 2003 .

[11]  R. Knight,et al.  Development of the human gastrointestinal microbiota and insights from high-throughput sequencing. , 2011, Gastroenterology.

[12]  R. Whittaker Evolution and measurement of species diversity , 1972 .

[13]  Adam P. Arkin,et al.  FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix , 2009, Molecular biology and evolution.

[14]  P. Bork,et al.  A human gut microbial gene catalogue established by metagenomic sequencing , 2010, Nature.

[15]  R. Knight,et al.  Supervised classification of human microbiota. , 2011, FEMS microbiology reviews.

[16]  C. Davis,et al.  Preliminary characterization of the normal microbiota of the human vulva using cultivation-independent methods. , 2007, Journal of medical microbiology.

[17]  E. Purdom,et al.  Diversity of the Human Intestinal Microbial Flora , 2005, Science.

[18]  Witold R. Rudnicki,et al.  Feature Selection with the Boruta Package , 2010 .

[19]  B. Roe,et al.  A core gut microbiome in obese and lean twins , 2008, Nature.

[20]  W. D. de Vos,et al.  Molecular Monitoring of Succession of Bacterial Communities in Human Neonates , 2002, Applied and Environmental Microbiology.

[21]  James R. Cole,et al.  The Ribosomal Database Project: improved alignments and new tools for rRNA analysis , 2008, Nucleic Acids Res..

[22]  Rodolphe Barrangou,et al.  The genome sequence of the probiotic intestinal bacterium Lactobacillus johnsonii NCC 533. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[23]  R. Knight,et al.  Forensic identification using skin bacterial communities , 2010, Proceedings of the National Academy of Sciences.

[24]  J. Gordon,et al.  Molecular analysis of commensal host-microbial relationships in the intestine. , 2001, Science.

[25]  F. Bäckhed,et al.  Intestinal microbiota during infancy and its implications for obesity. , 2009, Journal of pediatric gastroenterology and nutrition.

[26]  Rob Knight,et al.  Supervised classification of microbiota mitigates mislabeling errors , 2011, The ISME Journal.

[27]  R. Knight,et al.  Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns , 2010, Proceedings of the National Academy of Sciences.

[28]  R. Knight,et al.  Evolution of Mammals and Their Gut Microbes , 2008, Science.

[29]  Peer Bork,et al.  Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation , 2007, Bioinform..

[30]  M. Pop,et al.  Metagenomic Analysis of the Human Distal Gut Microbiome , 2006, Science.

[31]  P. A. van den Brandt,et al.  Factors Influencing the Composition of the Intestinal Microbiota in Early Infancy , 2006, Pediatrics.

[32]  Hiroshi Mori,et al.  Comparative Metagenomics Revealed Commonly Enriched Gene Sets in Human Gut Microbiomes , 2007, DNA research : an international journal for rapid publication of reports on genes and genomes.

[33]  B. Haas,et al.  Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. , 2011, Genome research.

[34]  Roeland Kindt,et al.  Tree diversity analysis: a manual and software for common statistical methods for ecological and biodiversity studies , 2006 .

[35]  P. Kero,et al.  Fecal microflora in healthy infants born by different methods of delivery: permanent changes in intestinal flora after cesarean delivery. , 1999, Journal of pediatric gastroenterology and nutrition.

[36]  J. Tiedje,et al.  Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy , 2007, Applied and Environmental Microbiology.

[37]  R. Knight,et al.  UniFrac: a New Phylogenetic Method for Comparing Microbial Communities , 2005, Applied and Environmental Microbiology.

[38]  P. Gajer,et al.  Recent Advances in Understanding the Microbiology of the Female Reproductive Tract and the Causes of Premature Birth , 2010, Infectious diseases in obstetrics and gynecology.

[39]  Sang Joon Kim,et al.  A Mathematical Theory of Communication , 2006 .

[40]  W. D. de Vos,et al.  Development of bacterial and bifidobacterial communities in feces of newborn babies. , 2003, Anaerobe.

[41]  R. Knight,et al.  Bacterial Community Variation in Human Body Habitats Across Space and Time , 2009, Science.

[42]  Jacques Ravel,et al.  Comparison of Self-Collected and Physician-Collected Vaginal Swabs for Microbiome Analysis , 2010, Journal of Clinical Microbiology.

[43]  E. Zoetendal,et al.  Temperature Gradient Gel Electrophoresis Analysis of 16S rRNA from Human Fecal Samples Reveals Stable and Host-Specific Communities of Active Bacteria , 1998, Applied and Environmental Microbiology.

[44]  P. Legendre,et al.  vegan : Community Ecology Package. R package version 1.8-5 , 2007 .

[45]  C. Huttenhower,et al.  Metagenomic biomarker discovery and explanation , 2011, Genome Biology.

[46]  G. Olsen,et al.  Heterogeneity of Vaginal Microbial Communities within Individuals , 2009, Journal of Clinical Microbiology.

[47]  James A Foster,et al.  The vaginal flora of healthy women is not always dominated by Lactobacillus species. , 2006, The Journal of infectious diseases.

[48]  P. Gajer,et al.  Vaginal microbiome of reproductive-age women , 2010, Proceedings of the National Academy of Sciences.

[49]  R. Knight,et al.  The Human Microbiome Project , 2007, Nature.

[50]  A. Chao Nonparametric estimation of the number of classes in a population , 1984 .

[51]  R. Mackie,et al.  Developmental microbial ecology of the neonatal gastrointestinal tract. , 1999, The American journal of clinical nutrition.

[52]  R. Ley,et al.  Ecological and Evolutionary Forces Shaping Microbial Diversity in the Human Intestine , 2006, Cell.

[53]  C. Quince,et al.  Accurate determination of microbial diversity from 454 pyrosequencing data , 2009, Nature Methods.

[54]  Béla Tóthmérész,et al.  Comparison of different methods for diversity ordering , 1995 .

[55]  B. Snel,et al.  Toward Automatic Reconstruction of a Highly Resolved Tree of Life , 2006, Science.

[56]  Les Dethlefsen,et al.  The Pervasive Effects of an Antibiotic on the Human Gut Microbiota, as Revealed by Deep 16S rRNA Sequencing , 2008, PLoS biology.

[57]  Jun Ma,et al.  The Genboree Microbiome Toolset and the analysis of 16S rRNA microbial sequences , 2012, BMC Bioinformatics.

[58]  Yuzhen Ye,et al.  Identification and quantification of abundant species from pyrosequences of 16S rRNA by consensus alignment , 2010, 2010 IEEE International Conference on Bioinformatics and Biomedicine (BIBM).