Autochthonous faecal viral transfer (FVT) impacts the murine microbiome after antibiotic perturbation.

BACKGROUND: It has become increasingly accepted that establishing and maintaining a complex and diverse gut microbiota is fundamental to human health. There are growing efforts to identify means of modulating and influencing the microbiota, especially in individuals who have experienced a disruption in their native microbiota. Faecal microbiota transplantation (FMT) is one method that restores diversity to the microbiota of an individual by introducing microbes from a healthy donor. FMT introduces the total microbial load into the recipient, including the bacteria, archaea, yeasts, protists and viruses. In this study, we investigated whether an autochthonous faecal viral transfer (FVT), in the form of a sterile faecal filtrate, could impact the recovery of a bacteriome disrupted by antibiotic treatment. RESULTS: Following antibiotic disruption of the bacteriome, test mice received an FVT harvested prior to antibiotic treatment, while control mice received a heat- and nuclease-treated FVT. In both groups of mice, the perturbed microbiome reverted over time to one more similar to the pre-treatment one. However, the bacteriomes of mice that received an FVT, in which bacteriophages predominate, separated from those of the control mice as determined by principal co-ordinate analysis (PCoA). Moreover, analysis of the differentially abundant taxa indicated a closer resemblance to the pre-treatment bacteriome in the test mice that had received an FVT. Similarly, metagenomic sequencing of the virome confirmed that faecal bacteriophages of FVT and control mice differed over time in both abundance and diversity, with the phages constituting the FVT persisting in mice that received them. CONCLUSIONS: An autochthonous virome transfer reshaped the bacteriomes of mice post-antibiotic treatment such that they more closely resembled the pre-antibiotic microbiota profile compared to mice that received non-viable phages. Thus, FVT may have a role in addressing antibiotic-associated microbiota alterations and potentially prevent the establishment of post-antibiotic infection. Given that bacteriophages are biologically inert in the absence of their host bacteria, they could form a safe and effective alternative to whole microbiota transplants that could be delivered during/following perturbation of the gut flora.

Exploratory analysis of covariation of microbiota-derived vitamin K and cognition in older adults.

BACKGROUND: Vitamin K has multiple important physiological roles, including blood coagulation and beneficial effects on myelin integrity in the brain. Some intestinal microbes possess the genes to produce vitamin K in the form of menaquinone (MK). MK appears in higher concentration in tissues, such as the brain, particularly MK4, than the dietary form of phylloquinone (PK). Lower PK concentrations have been reported in patients with Alzheimer disease while higher serum PK concentrations have been positively associated with verbal episodic memory. Despite knowledge of the importance of vitamin K for various health parameters, few studies have measured MK concentration and biosynthesis by gut commensals. OBJECTIVE: The aim of the current study was to investigate the relation between genes involved in gut-microbiota derived MK, concentrations of MK isoforms, and cognitive function. METHODS: Shotgun metagenomic sequencing of the gut microbiome of 74 elderly individuals with different cognitive ability levels was performed. From this, gene counts for microbial MK biosynthesis were determined. Associations between clusters of individuals, grouped based on a similar presence and prevalence of MK biosynthesis genes, and cognitive ability were investigated. Fecal MK concentrations were quantified by HPLC to investigate correlations with subject clusters. RESULTS: Separation of subject groups defined by banded quantification of the genetic potential of their microbiome to biosynthesize MK was associated with significant differences in cognitive ability [assessed using the Mini-Mental State Examination (MMSE)]. Three MK isoforms were found to be positively associated with MMSE, along with the identification of key components of the MK pathway that drive this association. Although the causality and direction of these associations remain unknown, these findings justify further studies. CONCLUSIONS: This study provides evidence that although total concentrations of MK did not covary with cognition, certain MK isoforms synthesized by the gut microbiome, particularly the longer chains, are positively associated with cognition.

Viromes of one year old infants reveal the impact of birth mode on microbiome diversity.

Establishing a diverse gut microbiota after birth is being increasingly recognised as important for preventing illnesses later in life. It is well established that bacterial diversity rapidly increases post-partum; however, few studies have examined the infant gut virome/phageome during this developmental period. We performed a metagenomic analysis of 20 infant faecal viromes at one year of age to determine whether spontaneous vaginal delivery (SVD) or caesarean section (CS) influenced viral composition. We find that birth mode results in distinctly different viral communities, with SVD infants having greater viral and bacteriophage diversity. We demonstrate that CrAssphage is acquired early in life, both in this cohort and two others, although no difference in birth mode is detected. A previous study has shown that bacterial OTU's (operational taxonomic units) identified in the same infants could not discriminate between birth mode at 12 months of age. Therefore, our results indicate that vertical transmission of viral communities from mother to child may play a role in shaping the early life microbiome, and that birth mode should be considered when studying the early life gut virome.

Reproducible protocols for metagenomic analysis of human faecal phageomes.

BACKGROUND: Recent studies have demonstrated that the human gut is populated by complex, highly individual and stable communities of viruses, the majority of which are bacteriophages. While disease-specific alterations in the gut phageome have been observed in IBD, AIDS and acute malnutrition, the human gut phageome remains poorly characterised. One important obstacle in metagenomic studies of the human gut phageome is a high level of discrepancy between results obtained by different research groups. This is often due to the use of different protocols for enriching virus-like particles, nucleic acid purification and sequencing. The goal of the present study is to develop a relatively simple, reproducible and cost-efficient protocol for the extraction of viral nucleic acids from human faecal samples, suitable for high-throughput studies. We also analyse the effect of certain potential confounding factors, such as storage conditions, repeated freeze-thaw cycles, and operator bias on the resultant phageome profile. Additionally, spiking of faecal samples with an exogenous phage standard was employed to quantitatively analyse phageomes following metagenomic sequencing. Comparative analysis of phageome profiles to bacteriome profiles was also performed following 16S rRNA amplicon sequencing. RESULTS: Faecal phageome profiles exhibit an overall greater individual specificity when compared to bacteriome profiles. The phageome and bacteriome both exhibited moderate change when stored at + 4 °C or room temperature. Phageome profiles were less impacted by multiple freeze-thaw cycles than bacteriome profiles, but there was a greater chance for operator effect in phageome processing. The successful spiking of faecal samples with exogenous bacteriophage demonstrated large variations in the total viral load between individual samples. CONCLUSIONS: The faecal phageome sequencing protocol developed in this study provides a valuable additional view of the human gut microbiota that is complementary to 16S amplicon sequencing and/or metagenomic sequencing of total faecal DNA. The protocol was optimised for several confounding factors that are encountered while processing faecal samples, to reduce discrepancies observed within and between research groups studying the human gut phageome. Rapid storage, limited freeze-thaw cycling and spiking of faecal samples with an exogenous phage standard are recommended for optimum results.

A long and abundant non-coding RNA in Lactobacillus salivarius.

Lactobacillus salivarius, found in the intestinal microbiota of humans and animals, is studied as an example of the sub-dominant intestinal commensals that may impart benefits upon their host. Strains typically harbour at least one megaplasmid that encodes functions contributing to contingency metabolism and environmental adaptation. RNA sequencing (RNA-seq)transcriptomic analysis of L. salivarius strain UCC118 identified the presence of a novel unusually abundant long non-coding RNA (lncRNA) encoded by the megaplasmid, and which represented more than 75 % of the total RNA-seq reads after depletion of rRNA species. The expression level of this 520 nt lncRNA in L. salivarius UCC118 exceeded that of the 16S rRNA, it accumulated during growth, was very stable over time and was also expressed during intestinal transit in a mouse. This lncRNA sequence is specific to the L. salivarius species; however, among 45 L. salivarius genomes analysed, not all (only 34) harboured the sequence for the lncRNA. This lncRNA was produced in 27 tested L. salivarius strains, but at strain-specific expression levels. High-level lncRNA expression correlated with high megaplasmid copy number. Transcriptome analysis of a deletion mutant lacking this lncRNA identified altered expression levels of genes in a number of pathways, but a definitive function of this new lncRNA was not identified. This lncRNA presents distinctive and unique properties, and suggests potential basic and applied scientific developments of this phenomenon.

Genomics and metagenomics of trimethylamine-utilizing Archaea in the human gut microbiome.

The biological significance of Archaea in the human gut microbiota is largely unclear. We recently reported genomic and biochemical analyses of the Methanomassiliicoccales, a novel order of methanogenic Archaea dwelling in soil and the animal digestive tract. We now show that these Methanomassiliicoccales are present in published microbiome data sets from eight countries. They are represented by five Operational Taxonomic Units present in at least four cohorts and phylogenetically distributed into two clades. Genes for utilizing trimethylamine (TMA), a bacterial precursor to an atherosclerogenic human metabolite, were present in four of the six novel Methanomassiliicoccales genomes assembled from ELDERMET metagenomes. In addition to increased microbiota TMA production capacity in long-term residential care subjects, abundance of TMA-utilizing Methanomassiliicoccales correlated positively with bacterial gene count for TMA production and negatively with fecal TMA concentrations. The two large Methanomassiliicoccales clades have opposite correlations with host health status in the ELDERMET cohort and putative distinct genomic signatures for gut adaptation.

Expanding the biotechnology potential of lactobacilli through comparative genomics of 213 strains and associated genera.

Lactobacilli are a diverse group of species that occupy diverse nutrient-rich niches associated with humans, animals, plants and food. They are used widely in biotechnology and food preservation, and are being explored as therapeutics. Exploiting lactobacilli has been complicated by metabolic diversity, unclear species identity and uncertain relationships between them and other commercially important lactic acid bacteria. The capacity for biotransformations catalysed by lactobacilli is an untapped biotechnology resource. Here we report the genome sequences of 213 Lactobacillus strains and associated genera, and their encoded genetic catalogue for modifying carbohydrates and proteins. In addition, we describe broad and diverse presence of novel CRISPR-Cas immune systems in lactobacilli that may be exploited for genome editing. We rationalize the phylogenomic distribution of host interaction factors and bacteriocins that affect their natural and industrial environments, and mechanisms to withstand stress during technological processes. We present a robust phylogenomic framework of existing species and for classifying new species.

The role of China in the global spread of the current cholera pandemic.

Epidemics and pandemics of cholera, a severe diarrheal disease, have occurred since the early 19th century and waves of epidemic disease continue today. Cholera epidemics are caused by individual, genetically monomorphic lineages of Vibrio cholerae: the ongoing seventh pandemic, which has spread globally since 1961, is associated with lineage L2 of biotype El Tor. Previous genomic studies of the epidemiology of the seventh pandemic identified three successive sub-lineages within L2, designated waves 1 to 3, which spread globally from the Bay of Bengal on multiple occasions. However, these studies did not include samples from China, which also experienced multiple epidemics of cholera in recent decades. We sequenced the genomes of 71 strains isolated in China between 1961 and 2010, as well as eight from other sources, and compared them with 181 published genomes. The results indicated that outbreaks in China between 1960 and 1990 were associated with wave 1 whereas later outbreaks were associated with wave 2. However, the previously defined waves overlapped temporally, and are an inadequate representation of the shape of the global genealogy. We therefore suggest replacing them by a series of tightly delineated clades. Between 1960 and 1990 multiple such clades were imported into China, underwent further microevolution there and then spread to other countries. China was thus both a sink and source during the pandemic spread of V. cholerae, and needs to be included in reconstructions of the global patterns of spread of cholera.

Detection and genomic characterization of motility in Lactobacillus curvatus: confirmation of motility in a species outside the Lactobacillus salivarius clade.

Lactobacillus is the largest genus within the lactic acid bacteria (LAB), with almost 180 species currently identified. Motility has been reported for at least 13 Lactobacillus species, all belonging to the Lactobacillus salivarius clade. Motility in lactobacilli is poorly characterized. It probably confers competitive advantages, such as superior nutrient acquisition and niche colonization, but it could also play an important role in innate immune system activation through flagellin­Toll-like receptor 5 (TLR5) interaction. We now report strong evidence of motility in a species outside the L. salivarius clade, Lactobacillus curvatus (strain NRIC0822). The motility of L. curvatus NRIC 0822 was revealed by phase-contrast microscopy and soft-agar motility assays. Strain NRIC 0822 was motile at temperatures between 15 °C and 37 °C, with a range of different carbohydrates, and under varying atmospheric conditions. We sequenced the L. curvatus NRIC 0822 genome, which revealed that the motility genes are organized in a single operon and that the products are very similar (>98.5% amino acid similarity over >11,000 amino acids) to those encoded by the motility operon of Lactobacillus acidipiscis KCTC 13900 (shown for the first time to be motile also). Moreover, the presence of a large number of mobile genetic elements within and flanking the motility operon of L. curvatus suggests recent horizontal transfer between members of two distinct Lactobacillus clades: L. acidipiscis in the L. salivarius clade and L. curvatus inthe L. sakei clade. This study provides novel phenotypic, genetic, and phylogenetic insights into flagellum-mediated motility in lactobacilli.

Transient Darwinian selection in Salmonella enterica serovar Paratyphi A during 450 years of global spread of enteric fever.

Multiple epidemic diseases have been designated as emerging or reemerging because the numbers of clinical cases have increased. Emerging diseases are often suspected to be driven by increased virulence or fitness, possibly associated with the gain of novel genes or mutations. However, the time period over which humans have been afflicted by such diseases is only known for very few bacterial pathogens, and the evidence for recently increased virulence or fitness is scanty. Has Darwinian (diversifying) selection at the genomic level recently driven microevolution within bacterial pathogens of humans? Salmonella enterica serovar Paratyphi A is a major cause of enteric fever, with a microbiological history dating to 1898. We identified seven modern lineages among 149 genomes on the basis of 4,584 SNPs in the core genome and estimated that Paratyphi A originated 450 y ago. During that time period, the effective population size has undergone expansion, reduction, and recent expansion. Mutations, some of which inactivate genes, have occurred continuously over the history of Paratyphi A, as has the gain or loss of accessory genes. We also identified 273 mutations that were under Darwinian selection. However, most genetic changes are transient, continuously being removed by purifying selection, and the genome of Paratyphi A has not changed dramatically over centuries. We conclude that Darwinian selection is not responsible for increased frequency of enteric fever and suggest that environmental changes may be more important for the frequency of disease.

Neutral genomic microevolution of a recently emerged pathogen, Salmonella enterica serovar Agona.

Salmonella enterica serovar Agona has caused multiple food-borne outbreaks of gastroenteritis since it was first isolated in 1952. We analyzed the genomes of 73 isolates from global sources, comparing five distinct outbreaks with sporadic infections as well as food contamination and the environment. Agona consists of three lineages with minimal mutational diversity: only 846 single nucleotide polymorphisms (SNPs) have accumulated in the non-repetitive, core genome since Agona evolved in 1932 and subsequently underwent a major population expansion in the 1960s. Homologous recombination with other serovars of S. enterica imported 42 recombinational tracts (360 kb) in 5/143 nodes within the genealogy, which resulted in 3,164 additional SNPs. In contrast to this paucity of genetic diversity, Agona is highly diverse according to pulsed-field gel electrophoresis (PFGE), which is used to assign isolates to outbreaks. PFGE diversity reflects a highly dynamic accessory genome associated with the gain or loss (indels) of 51 bacteriophages, 10 plasmids, and 6 integrative conjugational elements (ICE/IMEs), but did not correlate uniquely with outbreaks. Unlike the core genome, indels occurred repeatedly in independent nodes (homoplasies), resulting in inaccurate PFGE genealogies. The accessory genome contained only few cargo genes relevant to infection, other than antibiotic resistance. Thus, most of the genetic diversity within this recently emerged pathogen reflects changes in the accessory genome, or is due to recombination, but these changes seemed to reflect neutral processes rather than Darwinian selection. Each outbreak was caused by an independent clade, without universal, outbreak-associated genomic features, and none of the variable genes in the pan-genome seemed to be associated with an ability to cause outbreaks.

A genomic portrait of the emergence, evolution, and global spread of a methicillin-resistant Staphylococcus aureus pandemic.

The widespread use of antibiotics in association with high-density clinical care has driven the emergence of drug-resistant bacteria that are adapted to thrive in hospitalized patients. Of particular concern are globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clones that cause outbreaks and epidemics associated with health care. The most rapidly spreading and tenacious health-care-associated clone in Europe currently is EMRSA-15, which was first detected in the UK in the early 1990s and subsequently spread throughout Europe and beyond. Using phylogenomic methods to analyze the genome sequences for 193 S. aureus isolates, we were able to show that the current pandemic population of EMRSA-15 descends from a health-care-associated MRSA epidemic that spread throughout England in the 1980s, which had itself previously emerged from a primarily community-associated methicillin-sensitive population. The emergence of fluoroquinolone resistance in this EMRSA-15 subclone in the English Midlands during the mid-1980s appears to have played a key role in triggering pandemic spread, and occurred shortly after the first clinical trials of this drug. Genome-based coalescence analysis estimated that the population of this subclone over the last 20 yr has grown four times faster than its progenitor. Using comparative genomic analysis we identified the molecular genetic basis of 99.8% of the antimicrobial resistance phenotypes of the isolates, highlighting the potential of pathogen genome sequencing as a diagnostic tool. We document the genetic changes associated with adaptation to the hospital environment and with increasing drug resistance over time, and how MRSA evolution likely has been influenced by country-specific drug use regimens.

Archaea appear to dominate the microbiome of Inflatella pellicula deep sea sponges.

Microbes associated with marine sponges play significant roles in host physiology. Remarkable levels of microbial diversity have been observed in sponges worldwide through both culture-dependent and culture-independent studies. Most studies have focused on the structure of the bacterial communities in sponges and have involved sponges sampled from shallow waters. Here, we used pyrosequencing of 16S rRNA genes to compare the bacterial and archaeal communities associated with two individuals of the marine sponge Inflatella pellicula from the deep-sea, sampled from a depth of 2,900 m, a depth which far exceeds any previous sequence-based report of sponge-associated microbial communities. Sponge-microbial communities were also compared to the microbial community in the surrounding seawater. Sponge-associated microbial communities were dominated by archaeal sequencing reads with a single archaeal OTU, comprising ~60% and ~72% of sequences, being observed from Inflatella pellicula. Archaeal sequencing reads were less abundant in seawater (~11% of sequences). Sponge-associated microbial communities were less diverse and less even than any other sponge-microbial community investigated to date with just 210 and 273 OTUs (97% sequence identity) identified in sponges, with 4 and 6 dominant OTUs comprising ~88% and ~89% of sequences, respectively. Members of the candidate phyla, SAR406, NC10 and ZB3 are reported here from sponges for the first time, increasing the number of bacterial phyla or candidate divisions associated with sponges to 43. A minor cohort from both sponge samples (~0.2% and ~0.3% of sequences) were not classified to phylum level. A single OTU, common to both sponge individuals, dominates these unclassified reads and shares sequence homology with a sponge associated clone which itself has no known close relative and may represent a novel taxon.

The tree of genomes: an empirical comparison of genome-phylogeny reconstruction methods.

BACKGROUND: In the past decade or more, the emphasis for reconstructing species phylogenies has moved from the analysis of a single gene to the analysis of multiple genes and even completed genomes. The simplest method of scaling up is to use familiar analysis methods on a larger scale and this is the most popular approach. However, duplications and losses of genes along with horizontal gene transfer (HGT) can lead to a situation where there is only an indirect relationship between gene and genome phylogenies. In this study we examine five widely-used approaches and their variants to see if indeed they are more-or-less saying the same thing. In particular, we focus on Conditioned Reconstruction as it is a method that is designed to work well even if HGT is present. RESULTS: We confirm a previous suggestion that this method has a systematic bias. We show that no two methods produce the same results and most current methods of inferring genome phylogenies produce results that are significantly different to other methods. CONCLUSION: We conclude that genome phylogenies need to be interpreted differently, depending on the method used to construct them.