Consistency across multi-omics layers in a drug-perturbed gut microbial community

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

A faecal microbiota signature with high specificity for pancreatic cancer

BACKGROUND

Recent evidence suggests a role for the microbiome in pancreatic ductal adenocarcinoma (PDAC) aetiology and progression.

OBJECTIVE

To explore the faecal and salivary microbiota as potential diagnostic biomarkers.

METHODS

We applied shotgun metagenomic and 16S rRNA amplicon sequencing to samples from a Spanish case-control study (n=136), including 57 cases, 50 controls, and 29 patients with chronic pancreatitis in the discovery phase, and from a German case-control study (n=76), in the validation phase.

RESULTS

Faecal metagenomic classifiers performed much better than saliva-based classifiers and identified patients with PDAC with an accuracy of up to 0.84 area under the receiver operating characteristic curve (AUROC) based on a set of 27 microbial species, with consistent accuracy across early and late disease stages. Performance further improved to up to 0.94 AUROC when we combined our microbiome-based predictions with serum levels of carbohydrate antigen (CA) 19-9, the only current non-invasive, Food and Drug Administration approved, low specificity PDAC diagnostic biomarker. Furthermore, a microbiota-based classification model confined to PDAC-enriched species was highly disease-specific when validated against 25 publicly available metagenomic study populations for various health conditions (n=5792). Both microbiome-based models had a high prediction accuracy on a German validation population (n=76). Several faecal PDAC marker species were detectable in pancreatic tumour and non-tumour tissue using 16S rRNA sequencing and fluorescence in situ hybridisation.

CONCLUSION

Taken together, our results indicate that non-invasive, robust and specific faecal microbiota-based screening for the early detection of PDAC is feasible.

Antibiotics-induced monodominance of a novel gut bacterial order

OBJECTIVE

The composition of the healthy human adult gut microbiome is relatively stable over prolonged periods, and representatives of the most highly abundant and prevalent species have been cultured and described. However, microbial abundances can change on perturbations, such as antibiotics intake, enabling the identification and characterisation of otherwise low abundant species.

DESIGN

Analysing gut microbial time-series data, we used shotgun metagenomics to create strain level taxonomic and functional profiles. Community dynamics were modelled postintervention with a focus on conditionally rare taxa and previously unknown bacteria.

RESULTS

In response to a commonly prescribed cephalosporin (ceftriaxone), we observe a strong compositional shift in one subject, in which a previously unknown species, ^UBorkfalki ceftriaxensis, was identified, blooming to 92% relative abundance. The genome assembly reveals that this species (1) belongs to a so far undescribed order of Firmicutes, (2) is ubiquitously present at low abundances in at least one third of adults, (3) is opportunistically growing, being ecologically similar to typical probiotic species and (4) is stably associated to healthy hosts as determined by single nucleotide variation analysis. It was the first coloniser after the antibiotic intervention that led to a long-lasting microbial community shift and likely permanent loss of nine commensals.

CONCLUSION

The bloom of ^UB. ceftriaxensis and a subsequent one of Parabacteroides distasonis demonstrate the existence of monodominance community states in the gut. Our study points to an undiscovered wealth of low abundant but common taxa in the human gut and calls for more highly resolved longitudinal studies, in particular on ecosystem perturbations.

Sample Preservation and Storage Significantly Impact Taxonomic and Functional Profiles in Metaproteomics Studies of the Human Gut Microbiome

With the technological advances of the last decade, it is now feasible to analyze microbiome samples, such as human stool specimens, using multi-omic techniques. Given the inherent sample complexity, there exists a need for sample methods which preserve as much information as possible about the biological system at the time of sampling. Here, we analyzed human stool samples preserved and stored using different methods, applying metagenomics as well as metaproteomics. Our results demonstrate that sample preservation and storage have a significant effect on the taxonomic composition of identified proteins. The overall identification rates, as well as the proportion of proteins from Actinobacteria were much higher when samples were flash frozen. Preservation in RNAlater overall led to fewer protein identifications and a considerable increase in the share of Bacteroidetes, as well as Proteobacteria. Additionally, a decrease in the share of metabolism-related proteins and an increase of the relative amount of proteins involved in the processing of genetic information was observed for RNAlater-stored samples. This suggests that great care should be taken in choosing methods for the preservation and storage of microbiome samples, as well as in comparing the results of analyses using different sampling and storage methods. Flash freezing and subsequent storage at −80 °C should be chosen wherever possible.

Metagenomic analysis of gut microbial communities from a Central Asian population

OBJECTIVE

Changes in the gut microbiota are increasingly recognised to be involved in many diseases. This ecosystem is known to be shaped by many factors, including climate, geography, host nutrition, lifestyle and medication. Thus, knowledge of varying populations with different habits is important for a better understanding of the microbiome.

DESIGN

We therefore conducted a metagenomic analysis of intestinal microbiota from Kazakh donors, recruiting 84 subjects, including male and female healthy subjects and metabolic syndrome (MetS) patients aged 25-75 years, from the Kazakh administrative centre, Astana. We characterise and describe these microbiomes, the first deep-sequencing cohort from Central Asia, in comparison with a global dataset (832 individuals from five countries on three continents), and explore correlations between microbiota, clinical and laboratory parameters as well as with nutritional data from Food Frequency Questionnaires.

RESULTS

We observe that Kazakh microbiomes are relatively different from both European and East Asian counterparts, though similar to other Central Asian microbiomes, with the most striking difference being significantly more samples falling within the Prevotella-rich enterotype, potentially reflecting regional diet and lifestyle. We show that this enterotype designation remains stable within an individual over time in 82% of cases. We further observe gut microbiome features that distinguish MetS patients from controls (eg, significantly reduced Firmicutes to Bacteroidetes ratio, Bifidobacteria and Subdoligranulum, alongside increased Prevotella), though these overlap little with previously published reports and thus may reflect idiosyncrasies of the present cohort.

CONCLUSION

Taken together, this exploratory study describes gut microbiome data from an understudied population, providing a starting point for further comparative work on biogeography and research on widespread diseases.

TRIAL REGISTRATION NUMBER

ISRCTN37346212; Post-results.

Towards standards for human fecal sample processing in metagenomic studies

Technical variation in metagenomic analysis must be minimized to confidently assess the contributions of microbiota to human health. Here we tested 21 representative DNA extraction protocols on the same fecal samples and quantified differences in observed microbial community composition. We compared them with differences due to library preparation and sample storage, which we contrasted with observed biological variation within the same specimen or within an individual over time. We found that DNA extraction had the largest effect on the outcome of metagenomic analysis. To rank DNA extraction protocols, we considered resulting DNA quantity and quality, and we ascertained biases in estimates of community diversity and the ratio between Gram-positive and Gram-negative bacteria. We recommend a standardized DNA extraction method for human fecal samples, for which transferability across labs was established and which was further benchmarked using a mock community of known composition. Its adoption will improve comparability of human gut microbiome studies and facilitate meta-analyses.

Genotypes of SLC22A4 and SLC22A5 regulatory loci are predictive of the response of chronic myeloid leukemia patients to imatinib treatment

Background

Through high-throughput next-generation sequencing of promoters of solute carrier and ATP-binding cassette genes, which encode drug transporters, we aimed to identify SNPs associated with the response to imatinib administered for first-line treatment of patients with chronic myeloid leukemia.

Methods

In silico analysis using publicly available databases was done to select the SLC and ABC genes and their promoters for the next-generation sequencing. SNPs associated with the imatinib response were identified using Fisher’s exact probability tests and subjected to the linkage disequilibrium analyses with regulatory loci of concerned genes. We analyzed cumulative achievement of major molecular response and probability of event free survival in relation to identified SNP genotypes in 129 CML patients and performed multivariate analysis for determination of genotypes as independent predictors of outcome. Gene expression analysis of eight cell lines naturally carrying different genotypes was performed to outline an impact of genotypes on the gene expression.

Results

We observed significant differences in the frequencies of the rs460089-GC and rs460089-GG (SLC22A4) genotypes among rs2631365-TC (SLC22A5) genotype carriers that were associated with optimal and non-optimal responses, respectively. Loci rs460089 and rs2631365 were in highly significant linkage disequilibrium with 12 regulatory loci in introns of SLC22A4 and SLC22A5 encoding imatinib transporters. Genotype association analysis with the response to imatinib indicated that rs460089-GC carriers had a significantly higher probability of achieving a stable major molecular response (BCR-ABL1 transcript level below or equal to 0.1% in the international scale). In contrast, the rs460089-GG represented a risk factor for imatinib failure, which was significantly higher in rs460089-GG_rs2631365-TC carriers.

Conclusions

This exploratory study depicted potentially important genetic markers predicting outcome of imatinib treatment, which may be helpful for tailoring therapy in clinical practice.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-017-0523-3) contains supplementary material, which is available to authorized users.

Colorectal Cancer and the Human Gut Microbiome: Reproducibility with Whole-Genome Shotgun Sequencing

Accumulating evidence indicates that the gut microbiota affects colorectal cancer development, but previous studies have varied in population, technical methods, and associations with cancer. Understanding these variations is needed for comparisons and for potential pooling across studies. Therefore, we performed whole-genome shotgun sequencing on fecal samples from 52 pre-treatment colorectal cancer cases and 52 matched controls from Washington, DC. We compared findings from a previously published 16S rRNA study to the metagenomics-derived taxonomy within the same population. In addition, metagenome-predicted genes, modules, and pathways in the Washington, DC cases and controls were compared to cases and controls recruited in France whose specimens were processed using the same platform. Associations between the presence of fecal Fusobacteria, Fusobacterium, and Porphyromonas with colorectal cancer detected by 16S rRNA were reproduced by metagenomics, whereas higher relative abundance of Clostridia in cancer cases based on 16S rRNA was merely borderline based on metagenomics. This demonstrated that within the same sample set, most, but not all taxonomic associations were seen with both methods. Considering significant cancer associations with the relative abundance of genes, modules, and pathways in a recently published French metagenomics dataset, statistically significant associations in the Washington, DC population were detected for four out of 10 genes, three out of nine modules, and seven out of 17 pathways. In total, colorectal cancer status in the Washington, DC study was associated with 39% of the metagenome-predicted genes, modules, and pathways identified in the French study. More within and between population comparisons are needed to identify sources of variation and disease associations that can be reproduced despite these variations. Future studies should have larger sample sizes or pool data across studies to have sufficient power to detect associations that are reproducible and significant after correction for multiple testing.

Durable coexistence of donor and recipient strains after fecal microbiota transplantation

Fecal microbiota transplantation (FMT) has shown efficacy in treating recurrent Clostridium difficile infection and is increasingly being applied to other gastrointestinal disorders, yet the fate of native and introduced microbial strains remains largely unknown. To quantify the extent of donor microbiota colonization, we monitored strain populations in fecal samples from a recent FMT study on metabolic syndrome patients using single-nucleotide variants in metagenomes. We found extensive coexistence of donor and recipient strains, persisting 3 months after treatment. Colonization success was greater for conspecific strains than for new species, the latter falling within fluctuation levels observed in healthy individuals over a similar time frame. Furthermore, same-donor recipients displayed varying degrees of microbiota transfer, indicating individual patterns of microbiome resistance and donor-recipient compatibilities.

Potential of fecal microbiota for early-stage detection of colorectal cancer

Several bacterial species have been implicated in the development of colorectal carcinoma (CRC), but CRC-associated changes of fecal microbiota and their potential for cancer screening remain to be explored. Here, we used metagenomic sequencing of fecal samples to identify taxonomic markers that distinguished CRC patients from tumor-free controls in a study population of 156 participants. Accuracy of metagenomic CRC detection was similar to the standard fecal occult blood test (FOBT) and when both approaches were combined, sensitivity improved > 45% relative to the FOBT, while maintaining its specificity. Accuracy of metagenomic CRC detection did not differ significantly between early- and late-stage cancer and could be validated in independent patient and control populations (N = 335) from different countries. CRC-associated changes in the fecal microbiome at least partially reflected microbial community composition at the tumor itself, indicating that observed gene pool differences may reveal tumor-related host-microbe interactions. Indeed, we deduced a metabolic shift from fiber degradation in controls to utilization of host carbohydrates and amino acids in CRC patients, accompanied by an increase of lipopolysaccharide metabolism.

Identification of person and quantification of human DNA recovered from mosquitoes (Culicidae)

Mosquitoes (Culicidae) are widespread insects and can be important in forensic context as a source of human DNA. In order to establish the quantity of human DNA in mosquitoes' gut after different post-feeding interval and for how long after taking a bloodmeal the human donor could be identified, 174 blood-engorged mosquitoes (subfamily Anophelinae and Culicinae) were captured, kept alive and sacrificed at 8h intervals. Human DNA was amplified using forensic PCR kits (Identifiler, MiniFiler, and Quantifiler). A full DNA profiles were obtained from all Culicinae mosquitoes (74/74) up to 48 h and profiling was successful up to 88 h after a bloodmeal. Duration of post-feeding interval had a significant negative effect on the possibility of obtaining a full profile (p<0.001), and logistic regression found that the probability of obtaining the full profile is 15.5% less for every 8h increase in post-feeding interval. Culicinae mosquitoes are a suitable source of human DNA for forensic STR kits more than three days after a bloodmeal. Human DNA recovered from mosquito can be used for matching purposes and could be useful in revealing spatial and temporal relation of events that took place at the crime scene. Therefore, mosquitoes at the crime scene, dead or alive, could be a valuable piece of forensic evidence.