Deep Divergence and Genomic Diversification of Gut Symbionts of Neotropical Stingless Bees.

Social bees harbor conserved gut microbiotas that may have been acquired in a common ancestor of social bees and subsequently codiversified with their hosts. However, most of this knowledge is based on studies on the gut microbiotas of honey bees and bumblebees. Much less is known about the gut microbiotas of the third and most diverse group of social bees, the stingless bees. Specifically, the absence of genomic data from their microbiotas presents an important knowledge gap in understanding the evolution and functional diversity of the social bee microbiota. Here, we combined community profiling with culturing and genome sequencing of gut bacteria from six neotropical stingless bee species from Brazil. Phylogenomic analyses show that most stingless bee gut isolates form deep-branching sister clades of core members of the honey bee and bumblebee gut microbiota with conserved functional capabilities, confirming the common ancestry and ecology of their microbiota. However, our bacterial phylogenies were not congruent with those of the host, indicating that the evolution of the social bee gut microbiota was not driven by strict codiversification but included host switches and independent symbiont gain and losses. Finally, as reported for the honey bee and bumblebee microbiotas, we found substantial genomic divergence among strains of stingless bee gut bacteria, suggesting adaptation to different host species and glycan niches. Our study offers first insights into the genomic diversity of the stingless bee microbiota and highlights the need for broader samplings to understand the evolution of the social bee gut microbiota. IMPORTANCE Stingless bees are the most diverse group of the corbiculate bees and represent important pollinator species throughout the tropics and subtropics. They harbor specialized microbial communities in their gut that are related to those found in honey bees and bumblebees and that are likely important for bee health. Few bacteria have been cultured from the gut of stingless bees, which has prevented characterization of their genomic diversity and functional potential. Here, we established cultures of major members of the gut microbiotas of six stingless bee species and sequenced their genomes. We found that most stingless bee isolates belong to novel bacterial species distantly related to those found in honey bees and bumblebees and encoding similar functional capabilities. Our study offers a new perspective on the evolution of the social bee gut microbiota and presents a basis for characterizing the symbiotic relationships between gut bacteria and stingless bees.

Microbiome-mediated fructose depletion restricts murine gut colonization by vancomycin-resistant Enterococcus.

Multidrug-resistant organisms (MDRO) are a major threat to public health. MDRO infections, including those caused by vancomycin-resistant Enterococcus (VRE), frequently begin by colonization of the intestinal tract, a crucial step that is impaired by the intestinal microbiota. However, the specific members of the microbiota that suppress MDRO colonization and the mechanisms of such protection are largely unknown. Here, using metagenomics and mouse models that mimic the patients' exposure to antibiotics, we identified commensal bacteria associated with protection against VRE colonization. We further found a consortium of five strains that was sufficient to restrict VRE gut colonization in antibiotic treated mice. Transcriptomics in combination with targeted metabolomics and in vivo assays indicated that the bacterial consortium inhibits VRE growth through nutrient depletion, specifically by reducing the levels of fructose, a carbohydrate that boosts VRE growth in vivo. Finally, in vivo RNA-seq analysis of each strain of the consortium in combination with ex vivo and in vivo assays demonstrated that a single bacterium (Olsenella sp.) could recapitulate the effect of the consortium. Our results indicate that nutrient depletion by specific commensals can reduce VRE intestinal colonization, which represents a novel non-antibiotic based strategy to prevent infections caused by this multidrug-resistant organism.

Lactobacillus supports Clostridiales to restrict gut colonization by multidrug-resistant Enterobacteriaceae.

Infections by multidrug-resistant Enterobacteriaceae (MRE) are life-threatening to patients. The intestinal microbiome protects against MRE colonization, but antibiotics cause collateral damage to commensals and open the way to colonization and subsequent infection. Despite the significance of this problem, the specific commensals and mechanisms that restrict MRE colonization remain largely unknown. Here, by performing a multi-omic prospective study of hospitalized patients combined with mice experiments, we find that Lactobacillus is key, though not sufficient, to restrict MRE gut colonization. Lactobacillus rhamnosus and murinus increase the levels of Clostridiales bacteria, which induces a hostile environment for MRE growth through increased butyrate levels and reduced nutrient sources. This mechanism of colonization resistance, an interaction between Lactobacillus spp. and Clostridiales involving cooperation between microbiota members, is conserved in mice and patients. These results stress the importance of exploiting microbiome interactions for developing effective probiotics that prevent infections in hospitalized patients.

Listeria monocytogenes faecal carriage is common and depends on the gut microbiota.

Listeria genus comprises two pathogenic species, L. monocytogenes (Lm) and L. ivanovii, and non-pathogenic species. All can thrive as saprophytes, whereas only pathogenic species cause systemic infections. Identifying Listeria species' respective biotopes is critical to understand the ecological contribution of Listeria virulence. In order to investigate the prevalence and abundance of Listeria species in various sources, we retrieved and analyzed 16S rRNA datasets from MG-RAST metagenomic database. 26% of datasets contain Listeria sensu stricto sequences, and Lm is the most prevalent species, most abundant in soil and host-associated environments, including 5% of human stools. Lm is also detected in 10% of human stool samples from an independent cohort of 900 healthy asymptomatic donors. A specific microbiota signature is associated with Lm faecal carriage, both in humans and experimentally inoculated mice, in which it precedes Lm faecal carriage. These results indicate that Lm faecal carriage is common and depends on the gut microbiota, and suggest that Lm faecal carriage is a crucial yet overlooked consequence of its virulence.

Human commensal gut Proteobacteria withstand type VI secretion attacks through immunity protein-independent mechanisms.

While the major virulence factors for Vibrio cholerae, the cause of the devastating diarrheal disease cholera, have been extensively studied, the initial intestinal colonization of the bacterium is not well understood because non-human adult animals are refractory to its colonization. Recent studies suggest the involvement of an interbacterial killing device known as the type VI secretion system (T6SS). Here, we tested the T6SS-dependent interaction of V. cholerae with a selection of human gut commensal isolates. We show that the pathogen efficiently depleted representative genera of the Proteobacteria in vitro, while members of the Enterobacter cloacae complex and several Klebsiella species remained unaffected. We demonstrate that this resistance against T6SS assaults was mediated by the production of superior T6SS machinery or a barrier exerted by group I capsules. Collectively, our data provide new insights into immunity protein-independent T6SS resistance employed by the human microbiota and colonization resistance in general.

Macromolecular crowding links ribosomal protein gene dosage to growth rate in Vibrio cholerae.

BACKGROUND: In fast-growing bacteria, the genomic location of ribosomal protein (RP) genes is biased towards the replication origin (oriC). This trait allows optimizing their expression during exponential phase since oriC neighboring regions are in higher dose due to multifork replication. Relocation of s10-spc-α locus (S10), which codes for most of the RP, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction on its dosage, its expression, and bacterial growth rate. However, a mechanism linking S10 dosage to cell physiology has still not been determined. RESULTS: We hypothesized that S10 dosage perturbations impact protein synthesis capacity. Strikingly, we observed that in Vibrio cholerae, protein production capacity was independent of S10 position. Deep sequencing revealed that S10 relocation altered chromosomal replication dynamics and genome-wide transcription. Such changes increased as a function of oriC-S10 distance. Since RP constitutes a large proportion of cell mass, lower S10 dosage could lead to changes in macromolecular crowding, impacting cell physiology. Accordingly, cytoplasm fluidity was higher in mutants where S10 is most distant from oriC. In hyperosmotic conditions, when crowding differences are minimized, the growth rate and replication dynamics were highly alleviated in these strains. CONCLUSIONS: The genomic location of RP genes ensures its optimal dosage. However, besides of its essential function in translation, their genomic position sustains an optimal macromolecular crowding essential for maximizing growth. Hence, this could be another mechanism coordinating DNA replication to bacterial growth.

Community diversity and habitat structure shape the repertoire of extracellular proteins in bacteria.

We test the hypothesis that the frequency and cost of extracellular proteins produced by bacteria, which often depend on cooperative processes, vary with habitat structure and community diversity. The integration of the environmental distribution of bacteria (using 16S datasets) and their genomes shows that bacteria living in more structured habitats encode more extracellular proteins. In contrast, the effect of community diversity depends on protein function: it's positive for proteins implicated in antagonistic interactions and negative for those involved in nutrient acquisition. Extracellular proteins are costly and endure stronger selective pressure for low cost and for low diffusivity in less structured habitats and in more diverse communities. Finally, Bacteria found in multiple types of habitats, including host-associated generalists, encode more extracellular proteins than niche-restricted bacteria. These results show that ecological variables, notably habitat structure and community diversity, shape the evolution of the repertoires of genes encoding extracellular proteins and thus affect the ability of bacteria to manipulate their environment.

A simple, reproducible and cost-effective procedure to analyse gut phageome: from phage isolation to bioinformatic approach.

The microbiota of the human gut is a complex and rich community where bacteria and their viruses, the bacteriophages, are dominant. There are few studies on the phage community and no clear standard for isolating them, sequencing and analysing their genomes. Since this makes comparisons between studies difficult, we aimed at defining an easy, low-cost, and reproducible methodology. We analysed five different techniques to isolate phages from human adult faeces and developed an approach to analyse their genomes in order to quantify contamination and classify phage contigs in terms of taxonomy and lifestyle. We chose the polyethylene glycol concentration method to isolate phages because of its simplicity, low cost, reproducibility, and of the high number and diversity of phage sequences that we obtained. We also tested the reproducibility of this method with multiple displacement amplification (MDA) and showed that MDA severely decreases the phage genetic diversity of the samples and the reproducibility of the method. Lastly, we studied the influence of sequencing depth on the analysis of phage diversity and observed the beginning of a plateau for phage contigs at 20,000,000 reads. This work contributes to the development of methods for the isolation of phages in faeces and for their comparative analysis.

Challenges in Clinical Metaproteomics Highlighted by the Analysis of Acute Leukemia Patients with Gut Colonization by Multidrug-Resistant Enterobacteriaceae.

The microbiome has a strong impact on human health and disease and is, therefore, increasingly studied in a clinical context. Metaproteomics is also attracting considerable attention, and such data can be efficiently generated today owing to improvements in mass spectrometry-based proteomics. As we will discuss in this study, there are still major challenges notably in data analysis that need to be overcome. Here, we analyzed 212 fecal samples from 56 hospitalized acute leukemia patients with multidrug-resistant Enterobactericeae (MRE) gut colonization using metagenomics and metaproteomics. This is one of the largest clinical metaproteomic studies to date, and the first metaproteomic study addressing the gut microbiome in MRE colonized acute leukemia patients. Based on this substantial data set, we discuss major current limitations in clinical metaproteomic data analysis to provide guidance to researchers in the field. Notably, the results show that public metagenome databases are incomplete and that sample-specific metagenomes improve results. Furthermore, biological variation is tremendous which challenges clinical study designs and argues that longitudinal measurements of individual patients are a valuable future addition to the analysis of patient cohorts.

Gut colonization by a novel Clostridium species is associated with the onset of epizootic rabbit enteropathy.

Epizootic rabbit enteropathy (ERE) represents one of the most devastating diseases affecting rabbit farms. Previous studies showing transmissibility of disease symptoms through oral inoculation of intestinal contents from sick animals suggested a bacterial infectious origin for ERE. However, no etiological agent has been identified yet. On the other hand, ERE is associated with major changes in intestinal microbial communities, pinpointing dysbiosis as an alternative cause for the disease. To better understand the role of intestinal bacteria in ERE development, we have performed a prospective longitudinal study in which intestinal samples collected from the same animals before, during and after disease onset were analyzed using high-throughput sequencing. Changes in hundreds of bacterial groups were detected after the initiation of ERE. In contrast, before ERE onset, the microbiota from rabbits that developed ERE did not differ from those that remained healthy. Notably, an expansion of a single novel Clostridium species (Clostridium cuniculi) was detected the day of ERE onset. C. cuniculi encodes several putative toxins and it is phylogenetically related to the two well-characterized pathogens C. botulinum and C. perfringens. Our results are consistent with a bacterial infectious origin of ERE and discard dysbiosis as the initial trigger of the disease. Although experimental validation is required, results derived from sequencing analysis, propose a key role of C. cuniculi in ERE initiation.

Metagenomic assessment of the interplay between the environment and the genetic diversification of Acinetobacter.

Most bacteria have poorly characterized environmental reservoirs and unknown closely related species. This hampers the study of bacterial evolutionary ecology because both the environment and the genetic background of ancestral lineages are unknown. We combined metagenomics, comparative genomics and phylogenomics to overcome this limitation, to identify novel taxa and to propose environments where they can be isolated. We applied this method to characterize the ecological distribution of known and novel lineages of Acinetobacter spp. We observed two major environmental transitions at deep phylogenetic levels, splitting the genus into three ecologically differentiated clades. One of these has rapidly shifted towards host-association by acquiring genes involved in bacteria-eukaryote interactions. We show that environmental perturbations affect species distribution in predictable ways: bovines have very diverse communities of Acinetobacter, unless they were administered antibiotics, in which case they show highly uniform communities of Acinetobacter spp. that resemble those of humans. Our results uncover the diversity of bacterial lineages, overpassing the limitations of classical cultivation methods and highlight the role of the environment in shaping their evolution.

Abundance and co-occurrence of extracellular capsules increase environmental breadth: Implications for the emergence of pathogens.

Extracellular capsules constitute the outermost layer of many bacteria, are major virulence factors, and affect antimicrobial therapies. They have been used as epidemiological markers and recently became vaccination targets. Despite the efforts to biochemically serotype capsules in a few model pathogens, little is known of their taxonomic and environmental distribution. We developed, validated, and made available a computational tool, CapsuleFinder, to identify capsules in genomes. The analysis of over 2500 prokaryotic genomes, accessible in a database, revealed that ca. 50% of them-including Archaea-encode a capsule. The Wzx/Wzy-dependent capsular group was by far the most abundant. Surprisingly, a fifth of the genomes encode more than one capsule system-often from different groups-and their non-random co-occurrence suggests the existence of negative and positive epistatic interactions. To understand the role of multiple capsules, we queried more than 6700 metagenomes for the presence of species encoding capsules and showed that their distribution varied between environmental categories and, within the human microbiome, between body locations. Species encoding capsules, and especially those encoding multiple capsules, had larger environmental breadths than the other species. Accordingly, capsules were more frequent in environmental bacteria than in pathogens and, within the latter, they were more frequent among facultative pathogens. Nevertheless, capsules were frequent in clinical samples, and were usually associated with fast-growing bacteria with high infectious doses. Our results suggest that capsules increase the environmental range of bacteria and make them more resilient to environmental perturbations. Capsules might allow opportunistic pathogens to profit from empty ecological niches or environmental perturbations, such as those resulting from antibiotic therapy, to colonize the host. Capsule-associated virulence might thus be a by-product of environmental adaptation. Understanding the role of capsules in natural environments might enlighten their function in pathogenesis.

No Major Host Genetic Risk Factor Contributed to A(H1N1)2009 Influenza Severity.

While most patients affected by the influenza A(H1N1) pandemic experienced mild symptoms, a small fraction required hospitalization, often without concomitant factors that could explain such a severe course. We hypothesize that host genetic factors could contribute to aggravate the disease. To test this hypothesis, we compared the allele frequencies of 547,296 genome-wide single nucleotide polymorphisms (SNPs) between 49 severe and 107 mild confirmed influenza A cases, as well as against a general population sample of 549 individuals. When comparing severe vs. mild influenza A cases, only one SNP was close to the conventional p = 5×10-8. This SNP, rs28454025, sits in an intron of the GSK233 gene, which is involved in a neural development, but seems not to have any connections with immunological or inflammatory functions. Indirectly, a previous association reported with CD55 was replicated. Although sample sizes are low, we show that the statistical power in our design was sufficient to detect highly-penetrant, quasi-Mendelian genetic factors. Hence, and assuming that rs28454025 is likely to be a false positive, no major genetic factor was detected that could explain poor influenza A course.

Direct squencing from the minimal number of DNA molecules needed to fill a 454 picotiterplate.

The large amount of DNA needed to prepare a library in next generation sequencing protocols hinders direct sequencing of small DNA samples. This limitation is usually overcome by the enrichment of such samples with whole genome amplification (WGA), mostly by multiple displacement amplification (MDA) based on φ29 polymerase. However, this technique can be biased by the GC content of the sample and is prone to the development of chimeras as well as contamination during enrichment, which contributes to undesired noise during sequence data analysis, and also hampers the proper functional and/or taxonomic assignments. An alternative to MDA is direct DNA sequencing (DS), which represents the theoretical gold standard in genome sequencing. In this work, we explore the possibility of sequencing the genome of Escherichia coli fs 24 from the minimum number of DNA molecules required for pyrosequencing, according to the notion of one-bead-one-molecule. Using an optimized protocol for DS, we constructed a shotgun library containing the minimum number of DNA molecules needed to fill a selected region of a picotiterplate. We gathered most of the reference genome extension with uniform coverage. We compared the DS method with MDA applied to the same amount of starting DNA. As expected, MDA yielded a sparse and biased read distribution, with a very high amount of unassigned and unspecific DNA amplifications. The optimized DS protocol allows unbiased sequencing to be performed from samples with a very small amount of DNA.

Consistency of metagenomic assignment programs in simulated and real data.

BACKGROUND: Metagenomics is the genomic study of uncultured environmental samples, which has been greatly facilitated by the advent of shotgun-sequencing technologies. One of the main focuses of metagenomics is the discovery of previously uncultured microorganisms, which makes the assignment of sequences to a particular taxon a challenge and a crucial step. Recently, several methods have been developed to perform this task, based on different methodologies such as sequence composition or sequence similarity. The sequence composition methods have the ability to completely assign the whole dataset. However, their use in metagenomics and the study of their performance with real data is limited. In this work, we assess the consistency of three different methods (BLAST + Lowest Common Ancestor, Phymm, and Naïve Bayesian Classifier) in assigning real and simulated sequence reads. RESULTS: Both in real and in simulated data, BLAST + Lowest Common Ancestor (BLAST + LCA), Phymm, and Naïve Bayesian Classifier consistently assign a larger number of reads in higher taxonomic levels than in lower levels. However, discrepancies increase at lower taxonomic levels. In simulated data, consistent assignments between all three methods showed greater precision than assignments based on Phymm or Bayesian Classifier alone, since the BLAST + LCA algorithm performed best. In addition, assignment consistency in real data increased with sequence read length, in agreement with previously published simulation results. CONCLUSIONS: The use and combination of different approaches is advisable to assign metagenomic reads. Although the sensitivity could be reduced, the reliability can be increased by using the reads consistently assigned to the same taxa by, at least, two methods, and by training the programs using all available information.

A new method for extracting skin microbes allows metagenomic analysis of whole-deep skin.

In the last decade, an extensive effort has been made to characterize the human microbiota, due to its clinical and economic interests. However, a metagenomic approach to the skin microbiota is hampered by the high proportion of host DNA that is recovered. In contrast with the burgeoning field of gut metagenomics, skin metagenomics has been hindered by the absence of an efficient method to avoid sequencing the host DNA. We present here a method for recovering microbial DNA from skin samples, based on a combination of molecular techniques. We have applied this method to mouse skin, and have validated it by standard, quantitative PCR and amplicon sequencing of 16S rRNA. The taxonomic diversity recovered was not altered by this new method, as proved by comparing the phylogenetic structure revealed by 16S rRNA sequencing in untreated vs. treated samples. As proof of concept, we also present the first two mouse skin metagenomes, which allowed discovering new taxa (not only prokaryotes but also viruses and eukaryots) not reachable by 16S rRNA sequencing, as well as to characterize the skin microbiome functional landscape. Our method paves the way for the development of skin metagenomics, which will allow a much deeper knowledge of the skin microbiome and its relationship with the host, both in a healthy state and in relation to disease.

Staphylococcus prevails in the skin microbiota of long-term immunodeficient mice.

Host-commensal relationships in the skin are a complex system governed by variables related to the host, the bacteria and the environment. A disruption of this system may lead to new steady states, which, in turn, may lead to disease. We have studied one such disruption by characterizing the skin microbiota in healthy and immunodepressed (ID) mice. A detailed anatomopathological study failed to reveal any difference between the skin of healthy and ID mice. We sequenced the 16S rDNA V1-V2 gene region to saturation in 10 healthy and 10 ID 8 week-old mice, and found than all of the healthy and two of the ID mice had bacterial communities that were similar in composition to that of human skin, although, presumably because of the uniform raising conditions, less interindividual variation was found in mice. However, eight ID mice showed microbiota dominated by Staphylococcus epidermidis. Quantitative PCR amplification of 16S rDNA gene and of the Staphylococcus-specific TstaG region confirmed the previous results and indicated that the quantitative levels of Staphylococcus were similar in both groups while the total number of 16S copies was greater in the healthy mice. Thus, it is possible that, under long-term immunodeficiency, which removes the acquired but not the native immune system, S.epidermidis may inhibit the growth of other bacteria but does not cause a pathogenic state.

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics.

BACKGROUND: Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks. METHODOLOGY/PRINCIPALS FINDINGS: To investigate whether this pattern is present in ancient remains from a temperate environment, we have 454-FLX pyrosequenced different samples dated between 5,500 and 49,000 years ago: a bone from an extinct goat (Myotragus balearicus) that was treated with a depurinating agent (bleach), an Iberian lynx bone not subjected to any treatment, a human Neolithic sample from Barcelona (Spain), and a Neandertal sample from the El Sidrón site (Asturias, Spain). The efficiency of retrieval of endogenous sequences is below 1% in all cases. We have used the non-human samples to identify human sequences (0.35 and 1.4%, respectively), that we positively know are contaminants. CONCLUSIONS: We observed that bleach treatment appears to create a depurination-associated fragmentation pattern in resulting contaminant sequences that is indistinguishable from previously described endogenous sequences. Furthermore, the nucleotide composition pattern observed in 5' and 3' ends of contaminant sequences is much more complex than the flat pattern previously described in some Neandertal contaminants. Although much research on samples with known contaminant histories is needed, our results suggest that endogenous and contaminant sequences cannot be distinguished by the fragmentation pattern alone.