Human gut metagenomic mining reveals an untapped source of peptide antibiotics.

Drug-resistant bacteria are outpacing traditional antibiotic discovery efforts. Here, we computationally mined 444,054 families of putative small proteins from 1,773 human gut metagenomes, identifying 323 peptide antibiotics encoded in small open reading frames (smORFs). To test our computational predictions, 78 peptides were synthesized and screened for antimicrobial activity in vitro, with 59% displaying activity against either pathogens or commensals. Since these peptides were unique compared to previously reported antimicrobial peptides, we termed them smORF-encoded peptides (SEPs). SEPs killed bacteria by targeting their membrane, synergized with each other, and modulated gut commensals, indicating that they may play a role in reconfiguring microbiome communities in addition to counteracting pathogens. The lead candidates were anti-infective in both murine skin abscess and deep thigh infection models. Notably, prevotellin-2 from Prevotella copri presented activity comparable to the commonly used antibiotic polymyxin B. We report the discovery of hundreds of peptide sequences in the human gut.

Metagenomic compendium of 189,680 DNA viruses from the human gut microbiome.

Bacteriophages have important roles in the ecology of the human gut microbiome but are under-represented in reference databases. To address this problem, we assembled the Metagenomic Gut Virus catalogue that comprises 189,680 viral genomes from 11,810 publicly available human stool metagenomes. Over 75% of genomes represent double-stranded DNA phages that infect members of the Bacteroidia and Clostridia classes. Based on sequence clustering we identified 54,118 candidate viral species, 92% of which were not found in existing databases. The Metagenomic Gut Virus catalogue improves detection of viruses in stool metagenomes and accounts for nearly 40% of CRISPR spacers found in human gut Bacteria and Archaea. We also produced a catalogue of 459,375 viral protein clusters to explore the functional potential of the gut virome. This revealed tens of thousands of diversity-generating retroelements, which use error-prone reverse transcription to mutate target genes and may be involved in the molecular arms race between phages and their bacterial hosts.

MetaRibo-Seq measures translation in microbiomes.

No method exists to measure large-scale translation of genes in uncultured organisms in microbiomes. To overcome this limitation, we develop MetaRibo-Seq, a method for simultaneous ribosome profiling of tens to hundreds of organisms in microbiome samples. MetaRibo-Seq was benchmarked against gold-standard Ribo-Seq in a mock microbial community and applied to five different human fecal samples. Unlike RNA-Seq, Ribo-Seq signal of a predicted gene suggests it encodes a translated protein. We demonstrate two applications of this technique: First, MetaRibo-Seq identifies small genes, whose identification until now has been challenging. For example, MetaRibo-Seq identifies 2,091 translated, previously unannotated small protein families from five fecal samples, more than doubling the number of small proteins predicted to exist in this niche. Second, the combined application of RNA-Seq and MetaRibo-Seq identifies differences in the translation of transcripts. In summary, MetaRibo-Seq enables comprehensive translational profiling in microbiomes and identifies previously unannotated small proteins.

Large-Scale Analyses of Human Microbiomes Reveal Thousands of Small, Novel Genes.

Small proteins are traditionally overlooked due to computational and experimental difficulties in detecting them. To systematically identify small proteins, we carried out a comparative genomics study on 1,773 human-associated metagenomes from four different body sites. We describe >4,000 conserved protein families, the majority of which are novel; ∼30% of these protein families are predicted to be secreted or transmembrane. Over 90% of the small protein families have no known domain and almost half are not represented in reference genomes. We identify putative housekeeping, mammalian-specific, defense-related, and protein families that are likely to be horizontally transferred. We provide evidence of transcription and translation for a subset of these families. Our study suggests that small proteins are highly abundant and those of the human microbiome, in particular, may perform diverse functions that have not been previously reported.

DISARM is a widespread bacterial defence system with broad anti-phage activities.

The evolutionary pressure imposed by phage predation on bacteria and archaea has resulted in the development of effective anti-phage defence mechanisms, including restriction-modification and CRISPR-Cas systems. Here, we report on a new defence system, DISARM (defence island system associated with restriction-modification), which is widespread in bacteria and archaea. DISARM is composed of five genes, including a DNA methylase and four other genes annotated as a helicase domain, a phospholipase D (PLD) domain, a DUF1998 domain and a gene of unknown function. Engineering the Bacillus paralicheniformis 9945a DISARM system into Bacillus subtilis has rendered the engineered bacteria protected against phages from all three major families of tailed double-stranded DNA phages. Using a series of gene deletions, we show that four of the five genes are essential for DISARM-mediated defence, with the fifth (PLD) being redundant for defence against some of the phages. We further show that DISARM restricts incoming phage DNA and that the B. paralicheniformis DISARM methylase modifies host CCWGG motifs as a marker of self DNA akin to restriction-modification systems. Our results suggest that DISARM is a new type of multi-gene restriction-modification module, expanding the arsenal of defence systems known to be at the disposal of prokaryotes against their viruses.

BREX is a novel phage resistance system widespread in microbial genomes.

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR-Cas and restriction-modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six-gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon-like protease, an alkaline phosphatase domain protein, a putative RNA-binding protein, a DNA methylase, an ATPase-domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non-palindromic TAGGAG motifs in the bacterial genome guides self/non-self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction-modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX-like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti-BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.

Discovery of functional toxin/antitoxin systems in bacteria by shotgun cloning.

Toxin-antitoxin (TA) modules, composed of a toxic protein and a counteracting antitoxin, play important roles in bacterial physiology. We examined the experimental insertion of 1.5 million genes from 388 microbial genomes into an Escherichia coli host using more than 8.5 million random clones. This revealed hundreds of genes (toxins) that could only be cloned when the neighboring gene (antitoxin) was present on the same clone. Clustering of these genes revealed TA families widespread in bacterial genomes, some of which deviate from the classical characteristics previously described for such modules. Introduction of these genes into E. coli validated that the toxin toxicity is mitigated by the antitoxin. Infection experiments with T7 phage showed that two of the new modules can provide resistance against phage. Moreover, our experiments revealed an "antidefense" protein in phage T7 that neutralizes phage resistance. Our results expose active fronts in the arms race between bacteria and phage.

A vast collection of microbial genes that are toxic to bacteria.

In the process of clone-based genome sequencing, initial assemblies frequently contain cloning gaps that can be resolved using cloning-independent methods, but the reason for their occurrence is largely unknown. By analyzing 9,328,693 sequencing clones from 393 microbial genomes, we systematically mapped more than 15,000 genes residing in cloning gaps and experimentally showed that their expression products are toxic to the Escherichia coli host. A subset of these toxic sequences was further evaluated through a series of functional assays exploring the mechanisms of their toxicity. Among these genes, our assays revealed novel toxins and restriction enzymes, and new classes of small, non-coding toxic RNAs that reproducibly inhibit E. coli growth. Further analyses also revealed abundant, short, toxic DNA fragments that were predicted to suppress E. coli growth by interacting with the replication initiator DnaA. Our results show that cloning gaps, once considered the result of technical problems, actually serve as a rich source for the discovery of biotechnologically valuable functions, and suggest new modes of antimicrobial interventions.

Implications of Cheyne-Stokes breathing in advanced systolic heart failure.

BACKGROUND: Cheyne-Stokes breathing (CSB) has been associated with heart failure (HF) patients for many years; however, its true prevalence and its prognostic implications are still obscure. HYPOTHESIS: The goal of this study was to investigate the prevalence and the possible prognostic implications of nocturnal CSB in advanced heart failure patients. METHODS: We performed single night full polysomonography ambulatory sleep studies in 71 HF patients. We analyzed the patients' sleep studies, clinical and laboratory data, and 6 month mortality. RESULTS: A total of 71 chronic systolic HF patients were analyzed, 60 males, 11 females, age 65 +/- 13 years. Mean left ventricular ejection fraction was 27% +/- 11%. Short episodes of CSB (at least 3 min duration) were present in all patients, and mean CSB duration was 1 hour. CSB duration was associated significantly with both high serum levels of N-terminal prohormone brain natriuretic peptide (NT-proBNP) as well as with 6 month mortality. Log CSB time had a significant correlation with log NT-proBNP (r = 0.5, P<.0001). Based on median CSB duration, the Kaplan-Meier survival curve analysis showed significant association with 6 month mortality (P = .03). CONCLUSIONS: CSB prevalence in advanced HF patients is higher than previously reported and is associated with increased serum levels of NT-proBNP and higher 6 month mortality.

Pre-steady-state decoding of the Bicoid morphogen gradient.

Morphogen gradients are established by the localized production and subsequent diffusion of signaling molecules. It is generally assumed that cell fates are induced only after morphogen profiles have reached their steady state. Yet, patterning processes during early development occur rapidly, and tissue patterning may precede the convergence of the gradient to its steady state. Here we consider the implications of pre-steady-state decoding of the Bicoid morphogen gradient for patterning of the anterior-posterior axis of the Drosophila embryo. Quantitative analysis of the shift in the expression domains of several Bicoid targets (gap genes) upon alteration of bcd dosage, as well as a temporal analysis of a reporter for Bicoid activity, suggest that a transient decoding mechanism is employed in this setting. We show that decoding the pre-steady-state morphogen profile can reduce patterning errors caused by fluctuations in the rate of morphogen production. This can explain the surprisingly small shifts in gap and pair-rule gene expression domains observed in response to alterations in bcd dosage.

Just-in-time transcription program in metabolic pathways.

A primary goal of systems biology is to understand the design principles of the transcription networks that govern the timing of gene expression. Here we measured promoter activity for approximately 100 genes in parallel from living cells at a resolution of minutes and accuracy of 10%, based on GFP and Lux reporter libraries. Focusing on the amino-acid biosynthesis systems of Escherichia coli, we identified a previously unknown temporal expression program and expression hierarchy that matches the enzyme order in unbranched pathways. We identified two design principles: the closer the enzyme is to the beginning of the pathway, the shorter the response time of the activation of its promoter and the higher its maximal promoter activity. Mathematical analysis suggests that this 'just-in-time' (ref. 5) transcription program is optimal under constraints of rapidly reaching a production goal with minimal total enzyme production. Our findings suggest that metabolic regulation networks are designed to generate precision promoter timing and activity programs that can be understood using the engineering principles of production pipelines.