Despite the critical role of the human microbiota in health, our understanding of microbiota compositional dynamics during and after pregnancy is incomplete. We conducted a case-control study of 49 pregnant women, 15 of whom delivered preterm. From 40 of these women, we analyzed bacterial taxonomic composition of 3,767 specimens collected prospectively and weekly during gestation and monthly after delivery from the vagina, distal gut, saliva, and tooth/ gum. Linear mixed-effects modeling, medoid-based clustering, and Markov chain modeling were used to analyze community temporal trends, community structure, and vaginal community state transitions. Microbiota community taxonomic composition and diversity remained remarkably stable at all four body sites during pregnancy (P > 0.05 for trends over time). Prevalence of a Lactobacillus-poor vaginal community state type (CST 4) was inversely correlated with gestational age at delivery (P = 0.0039). Risk for preterm birth was more pronounced for subjects with CST 4 accompanied by elevated Gardnerella or Ureaplasma abundances. This finding was validated with a set of 246 vaginal specimens from nine women (four of whom delivered preterm). Most women experienced a postdelivery disturbance in the vaginal community characterized by a decrease in Lactobacillus species and an increase in diverse anaerobes such as Peptoniphilus, Prevotella, and Anaerococcus species. This disturbance was unrelated to gestational age at delivery and persisted for up to 1 y. These findings have important implications for predicting premature labor, a major global health problem, and for understanding the potential impact of a persistent, altered postpartum microbiota on maternal health, including outcomes of pregnancies following short interpregnancy intervals.16S rRNA gene | pregnancy | preterm birth | microbiome | premature labor

Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. Previous studies have suggested that the maternal vaginal microbiota contributes to the pathophysiology of PTB, but conflicting results in recent years have raised doubts. We conducted a study of PTB compared with term birth in two cohorts of pregnant women: one predominantly Caucasian (n = 39) at low risk for PTB, the second predominantly African American and at high-risk (n = 96). We profiled the taxonomic composition of 2,179 vaginal swabs collected prospectively and weekly during gestation using 16S rRNA gene sequencing. Previously proposed associations between PTB and lower Lactobacillus and higher Gardnerella abundances replicated in the low-risk cohort, but not in the high-risk cohort. High-resolution bioinformatics enabled taxonomic assignment to the species and subspecies levels, revealing that Lactobacillus crispatus was associated with low risk of PTB in both cohorts, while Lactobacillus iners was not, and that a subspecies clade of Gardnerella vaginalis explained the genus association with PTB. Patterns of cooccurrence between L. crispatus and Gardnerella were highly exclusive, while Gardnerella and L. iners often coexisted at high frequencies. We argue that the vaginal microbiota is better represented by the quantitative frequencies of these key taxa than by classifying communities into five community state types. Our findings extend and corroborate the association between the vaginal microbiota and PTB, demonstrate the benefits of high-resolution statistical bioinformatics in clinical microbiome studies, and suggest that previous conflicting results may reflect the different risk profile of women of black race.pregnancy | prematurity | vaginal microbiota | Lactobacillus | Gardnerella P reterm birth (PTB; delivery at <37 gestational wk) affects ≈12% of US births and is the leading cause of neonatal death and morbidity worldwide. Multiple lines of evidence support a role for the indigenous microbial communities of the mother (the maternal microbiota) in the pathophysiology of PTB. Microbial invasion of the amniotic cavity is one of the most frequent causes of spontaneous PTB (1), and the most common invading taxa are consistent with maternal origin (2-4). Bacterial vaginosis (BV), a condition involving an altered vaginal microbiota, has been consistently identified as a risk factor for PTB (5, 6). Multiple studies have also found chronic periodontitis, another condition associated with an altered microbiota, to be a risk factor for PTB (7,8).High-throughput sequencing methods have facilitated new lines of investigation into the microbial etiology of PTB (9, 10). Amplification and high-throughput sequencing of the 16S rRNA gene (metabarcoding) simultaneously measures the presence and relative abundance of thousands of bacterial taxa (composition), and resolves differences to the level of genus and sometimes species or subspecies. To date, metabarcoding studies of the relationship between the vaginal microbiota and PTB have y...

Bacteriophages typically have small genomes 1 and depend on their bacterial hosts for replication 2 . Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.Phages-viruses that infect bacteria-are considered distinct from cellular life owing to their inability to carry out most biological processes required for reproduction. They are agents of ecosystem change because they prey on specific bacterial populations, mediate lateral gene transfer, alter host metabolism and redistribute bacterially derived compounds through cell lysis 2-4 . They spread antibiotic resistance 5 and disperse pathogenicity factors that cause disease in humans and animals 6,7 . Most knowledge about phages is based on laboratorystudied examples, the vast majority of which have genomes that are a few tens of kb in length. Widely used isolation-based methods select against large phage particles, and they can be excluded from phage concentrates obtained by passage through 100-nm or 200-nm filters 1 . In 2017, only 93 isolated phages with genomes that were more than 200 kb in length were published 1 . Sequencing of whole-community DNA can uncover phage-derived fragments; however, large genomes can still escape detection owing to fragmentation 8 . A new clade of human-and animal-associated megaphages was recently described on the basis of genomes that were manually curated to completion from metagenomic datasets 9 . This finding prompted us to carry out a more-comprehensive analysis of microbial communities to evaluate the prevalence, diversity and ecosystem distribution of phages with large genomes. Previously, phages with genomes of more than 200 kb have been referred to as 'jumbophages' 1 or, in the case of phages with genomes of more than 500 kb, as megaphages 9 . As the set reconstructed here span both size ranges we refer to them simply as 'huge phage...

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems provide adaptive immunity against phage via spacer-encoded CRISPR RNAs that are complementary to invasive nucleic acids. Here, we challenge Streptococcus thermophilus with a bacteriophage, and used PCR-based metagenomics to monitor phage-derived spacers daily for 15 days in two experiments. Spacers that target the host chromosome are infrequent and strongly selected against, suggesting autoimmunity is lethal. In experiments that recover over half a million spacers, we observe early dominance by a few spacer sub-populations and rapid oscillations in sub-population abundances. In two CRISPR systems and in replicate experiments, a few spacers account for the majority of spacer sequences. Nearly all phage locations targeted by the acquired spacers have a proto-spacer adjacent motif (PAM), indicating PAMs are involved in spacer acquisition. We detect a strong and reproducible bias in the phage genome locations from which spacers derive. This may reflect selection for specific spacers based on location and effectiveness.

Current understanding of microorganism–virus interactions, which shape the evolution and functioning of Earth's ecosystems, is based primarily on cultivated organisms. Here we investigate thousands of viral and microbial genomes recovered using a cultivation-independent approach to study the frequency, variety and taxonomic distribution of viral defence mechanisms. CRISPR-Cas systems that confer microorganisms with immunity to viruses are present in only 10% of 1,724 sampled microorganisms, compared with previous reports of 40% occurrence in bacteria and 81% in archaea. We attribute this large difference to the lack of CRISPR-Cas systems across major bacterial lineages that have no cultivated representatives. We correlate absence of CRISPR-Cas with lack of nucleotide biosynthesis capacity and a symbiotic lifestyle. Restriction systems are well represented in these lineages and might provide both non-specific viral defence and access to nucleotides.

Well-studied innate immune systems exist throughout bacteria and archaea, but a more recently discovered genomic locus may offer prokaryotes surprising immunological adaptability. Mediated by a cassette-like genomic locus termed Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), the microbial adaptive immune system differs from its eukaryotic immune analogues by incorporating new immunities unidirectionally. CRISPR thus stores genomically recoverable timelines of virus-host coevolution in natural organisms refractory to laboratory cultivation. Here we combined a population genetic mathematical model of CRISPR-virus coevolution with six years of metagenomic sequencing to link the recoverable genomic dynamics of CRISPR loci to the unknown population dynamics of virus and host in natural communities. Metagenomic reconstructions in an acid-mine drainage system document CRISPR loci conserving ancestral immune elements to the base-pair across thousands of microbial generations. This ‘trailer-end conservation’ occurs despite rapid viral mutation and despite rapid prokaryotic genomic deletion. The trailer-ends of many reconstructed CRISPR loci are also largely identical across a population. ‘Trailer-end clonality’ occurs despite predictions of host immunological diversity due to negative frequency dependent selection (kill the winner dynamics). Statistical clustering and model simulations explain this lack of diversity by capturing rapid selective sweeps by highly immune CRISPR lineages. Potentially explaining ‘trailer-end conservation,’ we record the first example of a viral bloom overwhelming a CRISPR system. The polyclonal viruses bloom even though they share sequences previously targeted by host CRISPR loci. Simulations show how increasing random genomic deletions in CRISPR loci purges immunological controls on long-lived viral sequences, allowing polyclonal viruses to bloom and depressing host fitness. Our results thus link documented patterns of genomic conservation in CRISPR loci to an evolutionary advantage against persistent viruses. By maintaining old immunities, selection may be tuning CRISPR-mediated immunity against viruses reemerging from lysogeny or migration.

We analyzed near-complete population (composite) genomic sequences for coexisting acidophilic ironoxidizing Leptospirillum group II and III bacteria (phylum Nitrospirae) and an extrachromosomal plasmid from a Richmond Mine, Iron Mountain, CA, acid mine drainage biofilm. Community proteomic analysis of the genomically characterized sample and two other biofilms identified 64.6% and 44.9% of the predicted proteins of Leptospirillum groups II and III, respectively, and 20% of the predicted plasmid proteins. The bacteria share 92% 16S rRNA gene sequence identity and >60% of their genes, including integrated plasmid-like regions. The extrachromosomal plasmid carries conjugation genes with detectable sequence similarity to genes in the integrated conjugative plasmid, but only those on the extrachromosomal element were identified by proteomics. Both bacterial groups have genes for community-essential functions, including carbon fixation and biosynthesis of vitamins, fatty acids, and biopolymers (including cellulose); proteomic analyses reveal these activities. Both Leptospirillum types have multiple pathways for osmotic protection. Although both are motile, signal transduction and methyl-accepting chemotaxis proteins are more abundant in Leptospirillum group III, consistent with its distribution in gradients within biofilms. Interestingly, Leptospirillum group II uses a methyldependent and Leptospirillum group III a methyl-independent response pathway. Although only Leptospirillum group III can fix nitrogen, these proteins were not identified by proteomics. The abundances of core proteins are similar in all communities, but the abundance levels of unique and shared proteins of unknown function vary. Some proteins unique to one organism were highly expressed and may be key to the functional and ecological differentiation of Leptospirillum groups II and III.To understand how microorganisms contribute to biogeochemical cycling, it is necessary to determine the roles of uncultivated as well as cultivated groups and to establish how these roles vary during ecological succession and when environmental conditions change. Shotgun genomic sequencing (metagenomics) has opened new opportunities for cultureindependent studies of microbial communities. Examples include investigations of acid mine drainage (AMD) biofilm communities (4, 43, 75), symbiosis in a marine worm involving sulfur-oxidizing and sulfate-reducing bacteria (85), and enhanced biological phosphorous removal by sludge communities (32). From these genomic data sets, it has been possible to reconstruct aspects of the metabolism of individual organisms (32) and coexisting community members (29, 75) and to identify which organisms contribute community-essential functions (75). An interesting question relates to how differences in metabolic potential between organisms from the same lineage allow them to occupy distinct niches. Identification of potentially adaptive traits in closely related organisms is also important from an evolutionary perspective.Genomic data do not...

Northern acidic peatlands are important sources of atmospheric methane, yet the methanogens in them are poorly characterized. We examined methanogenic activities and methanogen populations at different depths in two peatlands, McLean bog (MB) and Chicago bog (CB). Both have acidic (pH 3.5-4.5) peat soils, but the pH of the deeper layers of CB is near-neutral, reflecting its previous existence as a neutral-pH fen. Acetotrophic and hydrogenotrophic methanogenesis could be stimulated in upper samples from both bogs, and phylotypes of methanogens using H2/CO2 (Methanomicrobiales) or acetate (Methanosarcinales) were identified in 16S rRNA gene clone libraries and by terminal restriction fragment length polymorphism (T-RFLP) analyses using a novel primer/restriction enzyme set that we developed. Particularly dominant in the upper layers was a clade in the Methanomicrobiales, called E2 here and the R10 or fen group elsewhere, estimated by quantitative polymerase chain reaction to be present at approximately 10(8) cells per gram of dry peat. Methanogenic activity was considerably lower in deeper samples from both bogs. The methanogen populations detected by T-RFLP in deeper portions of MB were mainly E2 and the uncultured euryarchaeal rice cluster (RC)-II group, whereas populations in the less acidic CB deep layers were considerably different, and included a Methanomicrobiales clade we call E1-E1', as well as RC-I, RC-II, marine benthic group D, and a new cluster that we call the subaqueous cluster. E2 was barely detectable in the deeper samples from CB, further evidence for the associations of most organisms in this group with acidic habitats.