Application of machine learning techniques for creating urban microbial fingerprints

Background: Research has found that human associated microbial communities play a role in homeostasis and the disruption of these communities may be important in an array of medical conditions. However outside of the human body many of these communities remain poorly studied. The Metagenomics and Metadesign of the Subways and Urban Biomes (MetaSUB) International Consortium is characterizing the microbiomes of urban environments with the aim to improve design of mass transit systems. As part of the CAMDA 2018 MetaSUB Forensics Challenge 311 city microbiome samples were provided to create urban microbial fingerprints, as well as a further 3 mystery datasets for validation. Results: MetaSUB samples were clustered using t-SNE in an unsupervised fashion to almost discrete groups, which upon inspection represented city of origin. Based on this clustering, geographically close metropolitan areas appear to display similar microbial profiles such as those of Auckland and Hamilton. Mystery unlabeled samples were provided part of the challenge. A random forest classifier built on the initial dataset of 311 samples was capable of correctly classifying 83.3% of the mystery samples to their city of origin. Random Forest analyses also identified features with the highest discriminatory power, ranking bacterial species such as Campylobacter jejuni and Staphylococcus argenteus as highly predictive of city of origin. The surface from which the sample was collected displayed little detectable impact on the microbial profiles in the data generated here. The proportion of reads classified per sample varied greatly and so de-novo assembly was applied to recover genomic fragments representing organisms not captured in reference databases. Conclusions: Current methods can differentiate urban microbiome profiles from each other with relative ease. De-novo assembly indicated that the MetaSUB metagenomic data contains adequate depth to recover metagenomic assembled genomes and that current databases are not sufficient to fully characterize urban microbiomes. Profiles found here indicate there may be a relationship between geographical distance between areas and the urban microbiome composition although this will need further research. The impact of these different profiles on public health is currently unknown but the MetaSUB consortium is uniquely suited to evaluate these and provide a roadmap for the inclusion of urban microbiome information for city planning and public health policy

Choice of assembly software has a critical impact on virome characterisation

peer-reviewedBackground The viral component of microbial communities plays a vital role in driving bacterial diversity, facilitating nutrient turnover and shaping community composition. Despite their importance, the vast majority of viral sequences are poorly annotated and share little or no homology to reference databases. As a result, investigation of the viral metagenome (virome) relies heavily on de novo assembly of short sequencing reads to recover compositional and functional information. Metagenomic assembly is particularly challenging for virome data, often resulting in fragmented assemblies and poor recovery of viral community members. Despite the essential role of assembly in virome analysis and difficulties posed by these data, current assembly comparisons have been limited to subsections of virome studies or bacterial datasets. Design This study presents the most comprehensive virome assembly comparison to date, featuring 16 metagenomic assembly approaches which have featured in human virome studies. Assemblers were assessed using four independent virome datasets, namely, simulated reads, two mock communities, viromes spiked with a known phage and human gut viromes. Results Assembly performance varied significantly across all test datasets, with SPAdes (meta) performing consistently well. Performance of MIRA and VICUNA varied, highlighting the importance of using a range of datasets when comparing assembly programs. It was also found that while some assemblers addressed the challenges of virome data better than others, all assemblers had limitations. Low read coverage and genomic repeats resulted in assemblies with poor genome recovery, high degrees of fragmentation and low-accuracy contigs across all assemblers. These limitations must be considered when setting thresholds for downstream analysis and when drawing conclusions from virome data

Tumour-associated and non-tumour-associated microbiota in colorectal cancer

Objective: A signature that unifies the colorectal cancer (CRC) microbiota across multiple studies has not been identified. In addition to methodological variance, heterogeneity may be caused by both microbial and host response differences, which was addressed in this study. Design: We prospectively studied the colonic microbiota and the expression of specific host response genes using faecal and mucosal samples (‘ON’ and ‘OFF’ the tumour, proximal and distal) from 59 patients undergoing surgery for CRC, 21 individuals with polyps and 56 healthy controls. Microbiota composition was determined by 16S rRNA amplicon sequencing; expression of host genes involved in CRC progression and immune response was quantified by real-time quantitative PCR. Results: The microbiota of patients with CRC differed from that of controls, but alterations were not restricted to the cancerous tissue. Differences between distal and proximal cancers were detected and faecal microbiota only partially reflected mucosal microbiota in CRC. Patients with CRC can be stratified based on higher level structures of mucosal-associated bacterial co-abundance groups (CAGs) that resemble the previously formulated concept of enterotypes. Of these, Bacteroidetes Cluster 1 and Firmicutes Cluster 1 were in decreased abundance in CRC mucosa, whereas Bacteroidetes Cluster 2, Firmicutes Cluster 2, Pathogen Cluster and Prevotella Cluster showed increased abundance in CRC mucosa. CRC-associated CAGs were differentially correlated with the expression of host immunoinflammatory response genes. Conclusions: CRC-associated microbiota profiles differ from those in healthy subjects and are linked with distinct mucosal gene-expression profiles. Compositional alterations in the microbiota are not restricted to cancerous tissue and differ between distal and proximal cancers

Social interaction-induced activation of RNA splicing in the amygdala of microbiome-deficient mice

Social behaviour is regulated by activity of host-associated microbiota across multiple species. However, the molecular mechanisms mediating this relationship remain elusive. We therefore determined the dynamic, stimulus-dependent transcriptional regulation of germ-free (GF) and GF mice colonised post weaning (exGF) in the amygdala, a brain region critically involved in regulating social interaction. In GF mice the dynamic response seen in controls was attenuated and replaced by a marked increase in expression of splicing factors and alternative exon usage in GF mice upon stimulation, which was even more pronounced in exGF mice. In conclusion, we demonstrate a molecular basis for how the host microbiome is crucial for a normal behavioural response during social interaction. Our data further suggest that social behaviour is correlated with the gene-expression response in the amygdala, established during neurodevelopment as a result of host-microbe interactions. Our findings may help toward understanding neurodevelopmental events leading to social behaviour dysregulation, such as those found in autism spectrum disorders (ASDs)

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

peer-reviewedEstablishing 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

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

Regulation of prefrontal cortex myelination by the microbiota

The prefrontal cortex (PFC) is a key region implicated in a range of neuropsychiatric disorders such as depression, schizophrenia and autism. In parallel, the role of the gut microbiota in contributing to these disorders is emerging. Germ-free (GF) animals, microbiota-deficient throughout life, have been instrumental in elucidating the role of the microbiota in many aspects of physiology, especially the role of the microbiota in anxiety-related behaviours, impaired social cognition and stress responsivity. Here we aim to further elucidate the mechanisms of the microbial influence by investigating changes in the homeostatic regulation of neuronal transcription of GF mice within the PFC using a genome-wide transcriptome profiling approach. Our results reveal a marked, concerted upregulation of genes linked to myelination and myelin plasticity. This coincided with upregulation of neural activity-induced pathways, potentially driving myelin plasticity. Subsequent investigation at the ultrastructural level demonstrated the presence of hypermyelinated axons within the PFC of GF mice. Notably, these changes in myelin and activity-related gene expression could be reversed by colonization with a conventional microbiota following weaning. In summary, we believe we demonstrate for the first time that the microbiome is necessary for appropriate and dynamic regulation of myelin-related genes with clear implications for cortical myelination at an ultrastructural level. The microbiota is therefore a potential therapeutic target for psychiatric disorders involving dynamic myelination in the PFC

Reproducible protocols for metagenomic analysis of human faecal phageomes

peer-reviewedAll sequence data used in the analyses were deposited in the Sequence read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra) under BioProject PRJNA407341. Sample IDs, meta data and corresponding accession numbers are summarised in Additional file 2: Table S2. All raw count tables, 16S taxonomic assignments, BLAST top hits for viral contigs and R code used for the analysis are available at (https://figshare.com/s/71163558b4f78e3e7ed6).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