Combinatorial, additive and dose-dependent drug-microbiome associations.

During the transition from a healthy state to cardiometabolic disease, patients become heavily medicated, which leads to an increasingly aberrant gut microbiome and serum metabolome, and complicates biomarker discovery1-5. Here, through integrated multi-omics analyses of 2,173 European residents from the MetaCardis cohort, we show that the explanatory power of drugs for the variability in both host and gut microbiome features exceeds that of disease. We quantify inferred effects of single medications, their combinations as well as additive effects, and show that the latter shift the metabolome and microbiome towards a healthier state, exemplified in synergistic reduction in serum atherogenic lipoproteins by statins combined with aspirin, or enrichment of intestinal Roseburia by diuretic agents combined with beta-blockers. Several antibiotics exhibit a quantitative relationship between the number of courses prescribed and progression towards a microbiome state that is associated with the severity of cardiometabolic disease. We also report a relationship between cardiometabolic drug dosage, improvement in clinical markers and microbiome composition, supporting direct drug effects. Taken together, our computational framework and resulting resources enable the disentanglement of the effects of drugs and disease on host and microbiome features in multimedicated individuals. Furthermore, the robust signatures identified using our framework provide new hypotheses for drug-host-microbiome interactions in cardiometabolic disease.

Human gut microbes impact host serum metabolome and insulin sensitivity.

Insulin resistance is a forerunner state of ischaemic cardiovascular disease and type 2 diabetes. Here we show how the human gut microbiome impacts the serum metabolome and associates with insulin resistance in 277 non-diabetic Danish individuals. The serum metabolome of insulin-resistant individuals is characterized by increased levels of branched-chain amino acids (BCAAs), which correlate with a gut microbiome that has an enriched biosynthetic potential for BCAAs and is deprived of genes encoding bacterial inward transporters for these amino acids. Prevotella copri and Bacteroides vulgatus are identified as the main species driving the association between biosynthesis of BCAAs and insulin resistance, and in mice we demonstrate that P. copri can induce insulin resistance, aggravate glucose intolerance and augment circulating levels of BCAAs. Our findings suggest that microbial targets may have the potential to diminish insulin resistance and reduce the incidence of common metabolic and cardiovascular disorders.

Effects of Low FODMAP Diet on Symptoms, Fecal Microbiome, and Markers of Inflammation in Patients With Quiescent Inflammatory Bowel Disease in a Randomized Trial.

BACKGROUND & AIMS: There is limited evidence that a diet low in fermentable oligosaccharides, disaccharides, monosaccharides, and polyols (FODMAPs) reduces gut symptoms in quiescent inflammatory bowel disease (IBD). We performed a randomized, controlled trial to investigate the effects of a low FODMAP diet on persistent gut symptoms, the intestinal microbiome, and circulating markers of inflammation in patients with quiescent IBD. METHODS: We performed a single-blind trial of 52 patients with quiescent Crohn's disease or ulcerative colitis and persistent gut symptoms at 2 large gastroenterology clinics in the United Kingdom. Patients were randomly assigned to groups that followed a diet low in FODMAPs (n = 27) or a control diet (n = 25), with dietary advice, for 4 weeks. Gut symptoms and health-related quality of life were measured using validated questionnaires. Stool and blood samples were collected at baseline and end of trial. We assessed fecal microbiome composition and function using shotgun metagenomic sequencing and phenotypes of T cells in blood using flow cytometry. RESULTS: A higher proportion of patients reported adequate relief of gut symptoms following the low FODMAP diet (14/27, 52%) than the control diet (4/25, 16%, P=.007). Patients had a greater reduction in irritable bowel syndrome severity scores following the low FODMAP diet (mean reduction of 67; standard error, 78) than the control diet (mean reduction of 34; standard error, 50), although this difference was not statistically significant (P = .075). Following the low FODMAP diet, patients had higher health-related quality of life scores (81.9 ± 1.2) than patients on the control diet (78.3 ± 1.2, P = .042). A targeted analysis revealed that in stool samples collected at the end of the study period, patients on the low FODMAP diet had significantly lower abundance of Bifidobacterium adolescentis, Bifidobacterium longum, and Faecalibacterium prausnitzii than patients on control diet. However, microbiome diversity and markers of inflammation did not differ significantly between groups. CONCLUSIONS: In a trial of the low FODMAP diet vs a control diet in patients with quiescent IBD, we found no significant difference after 4 weeks in change in irritable bowel syndrome severity scores, but significant improvements in specific symptom scores and numbers reporting adequate symptom relief. The low FODMAP diet reduced fecal abundance of microbes believed to regulate the immune response, compared with the control diet, but had no significant effect on markers of inflammation. We conclude that a 4-week diet low in FODMAPs is safe and effective for managing persistent gut symptoms in patients with quiescent IBD. www.isrctn.com no.: ISRCTN17061468.

Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota.

In recent years, several associations between common chronic human disorders and altered gut microbiome composition and function have been reported. In most of these reports, treatment regimens were not controlled for and conclusions could thus be confounded by the effects of various drugs on the microbiota, which may obscure microbial causes, protective factors or diagnostically relevant signals. Our study addresses disease and drug signatures in the human gut microbiome of type 2 diabetes mellitus (T2D). Two previous quantitative gut metagenomics studies of T2D patients that were unstratified for treatment yielded divergent conclusions regarding its associated gut microbial dysbiosis. Here we show, using 784 available human gut metagenomes, how antidiabetic medication confounds these results, and analyse in detail the effects of the most widely used antidiabetic drug metformin. We provide support for microbial mediation of the therapeutic effects of metformin through short-chain fatty acid production, as well as for potential microbiota-mediated mechanisms behind known intestinal adverse effects in the form of a relative increase in abundance of Escherichia species. Controlling for metformin treatment, we report a unified signature of gut microbiome shifts in T2D with a depletion of butyrate-producing taxa. These in turn cause functional microbiome shifts, in part alleviated by metformin-induced changes. Overall, the present study emphasizes the need to disentangle gut microbiota signatures of specific human diseases from those of medication.

Commensal Hafnia alvei strain reduces food intake and fat mass in obese mice-a new potential probiotic for appetite and body weight management.

BACKGROUND/OBJECTIVES: Based on the recent identification of E.coli heat shock protein ClpB as a mimetic of the anorexigenic α-melanocyte stimulating hormone (α-MSH), the objective of this study was to preclinically validate Hafnia alvei, a ClpB-producing commensal bacterium as a potential probiotic for appetite and body weight management in overweight and obesity. METHODS: The involvement of enterobacterial ClpB in the putative anti-obesity effects was studied using ClpB-deficient E.coli. A food-grade H. alvei HA4597 strain synthetizing the ClpB protein with an α-MSH-like motif was selected as a candidate probiotic to be tested in ob/ob and high-fat diet (HFD)-fed obese and overweight mice. The relevance of the enterobacterial ClpB gene to human obesity was studied by in silico analysis of fecal metagenomes of 569 healthy individuals from the "MetaHIT" database. RESULTS: Chronic per os administration of native but not ClpB-deficient E.coli strain reduced body weight gain (p < 0.05) and daily meal frequency (p < 0.001) in ob/ob mice. Oral gavage of H.alvei for 18 and 46 days in ob/ob and HFD-fed obese mice, respectively, was well tolerated, reduced body weight gain and fat mass in both obesity models (p < 0.05) and decreased food intake in hyperphagic ob/ob mice (p < 0.001). Elevated fat tissue levels of phosphorylated hormone-sensitive lipase were detected in H.alvei -treated ob/ob mice (p < 0.01). Enterobacterial ClpB gene richness was lower in obese vs. non-obese humans (p < 0.0001) and correlated negatively with BMI in genera of Enterobacter, Klebsiella and Hafnia. CONCLUSIONS: H.alvei HA4597 strain reduces food intake, body weight and fat mass gain in hyperphagic and obese mice. These data combined with low enterobacterial ClpB gene abundance in the microbiota of obese humans provide the rationale for using H.alvei as a probiotic for appetite and body weight management in overweight and obesity.

MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data.

MOTIVATION: Analysis toolkits for shotgun metagenomic data achieve strain-level characterization of complex microbial communities by capturing intra-species gene content variation. Yet, these tools are hampered by the extent of reference genomes that are far from covering all microbial variability, as many species are still not sequenced or have only few strains available. Binning co-abundant genes obtained from de novo assembly is a powerful reference-free technique to discover and reconstitute gene repertoire of microbial species. While current methods accurately identify species core parts, they miss many accessory genes or split them into small gene groups that remain unassociated to core clusters. RESULTS: We introduce MSPminer, a computationally efficient software tool that reconstitutes Metagenomic Species Pan-genomes (MSPs) by binning co-abundant genes across metagenomic samples. MSPminer relies on a new robust measure of proportionality coupled with an empirical classifier to group and distinguish not only species core genes but accessory genes also. Applied to a large scale metagenomic dataset, MSPminer successfully delineates in a few hours the gene repertoires of 1661 microbial species with similar specificity and higher sensitivity than existing tools. The taxonomic annotation of MSPs reveals microorganisms hitherto unknown and brings coherence in the nomenclature of the species of the human gut microbiota. The provided MSPs can be readily used for taxonomic profiling and biomarkers discovery in human gut metagenomic samples. In addition, MSPminer can be applied on gene count tables from other ecosystems to perform similar analyses. AVAILABILITY AND IMPLEMENTATION: The binary is freely available for non-commercial users at www.enterome.com/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alterations of the human gut microbiome in liver cirrhosis.

Liver cirrhosis occurs as a consequence of many chronic liver diseases that are prevalent worldwide. Here we characterize the gut microbiome in liver cirrhosis by comparing 98 patients and 83 healthy control individuals. We build a reference gene set for the cohort containing 2.69 million genes, 36.1% of which are novel. Quantitative metagenomics reveals 75,245 genes that differ in abundance between the patients and healthy individuals (false discovery rate < 0.0001) and can be grouped into 66 clusters representing cognate bacterial species; 28 are enriched in patients and 38 in control individuals. Most (54%) of the patient-enriched, taxonomically assigned species are of buccal origin, suggesting an invasion of the gut from the mouth in liver cirrhosis. Biomarkers specific to liver cirrhosis at gene and function levels are revealed by a comparison with those for type 2 diabetes and inflammatory bowel disease. On the basis of only 15 biomarkers, a highly accurate patient discrimination index is created and validated on an independent cohort. Thus microbiota-targeted biomarkers may be a powerful tool for diagnosis of different diseases.

Correction: Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa.

The fortieth author's name was listed incorrectly. The correct presentation is A Keski-Rahkonen.

Widespread white matter microstructural differences in schizophrenia across 4322 individuals: results from the ENIGMA Schizophrenia DTI Working Group.

The regional distribution of white matter (WM) abnormalities in schizophrenia remains poorly understood, and reported disease effects on the brain vary widely between studies. In an effort to identify commonalities across studies, we perform what we believe is the first ever large-scale coordinated study of WM microstructural differences in schizophrenia. Our analysis consisted of 2359 healthy controls and 1963 schizophrenia patients from 29 independent international studies; we harmonized the processing and statistical analyses of diffusion tensor imaging (DTI) data across sites and meta-analyzed effects across studies. Significant reductions in fractional anisotropy (FA) in schizophrenia patients were widespread, and detected in 20 of 25 regions of interest within a WM skeleton representing all major WM fasciculi. Effect sizes varied by region, peaking at (d=0.42) for the entire WM skeleton, driven more by peripheral areas as opposed to the core WM where regions of interest were defined. The anterior corona radiata (d=0.40) and corpus callosum (d=0.39), specifically its body (d=0.39) and genu (d=0.37), showed greatest effects. Significant decreases, to lesser degrees, were observed in almost all regions analyzed. Larger effect sizes were observed for FA than diffusivity measures; significantly higher mean and radial diffusivity was observed for schizophrenia patients compared with controls. No significant effects of age at onset of schizophrenia or medication dosage were detected. As the largest coordinated analysis of WM differences in a psychiatric disorder to date, the present study provides a robust profile of widespread WM abnormalities in schizophrenia patients worldwide. Interactive three-dimensional visualization of the results is available at www.enigma-viewer.org.

Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa.

Anorexia nervosa (AN) is a complex neuropsychiatric disorder presenting with dangerously low body weight, and a deep and persistent fear of gaining weight. To date, only one genome-wide significant locus associated with AN has been identified. We performed an exome-chip based genome-wide association studies (GWAS) in 2158 cases from nine populations of European origin and 15 485 ancestrally matched controls. Unlike previous studies, this GWAS also probed association in low-frequency and rare variants. Sixteen independent variants were taken forward for in silico and de novo replication (11 common and 5 rare). No findings reached genome-wide significance. Two notable common variants were identified: rs10791286, an intronic variant in OPCML (P=9.89 × 10-6), and rs7700147, an intergenic variant (P=2.93 × 10-5). No low-frequency variant associations were identified at genome-wide significance, although the study was well-powered to detect low-frequency variants with large effect sizes, suggesting that there may be no AN loci in this genomic search space with large effect sizes.

Comparison of epidemiological, clinical and microbiological characteristics of bloodstream infection in children with solid tumours and haematological malignancies.

Bloodstream infection (BSI) is a serious complication in immunocompromised hosts. This study compares epidemiological, clinical and microbiological characteristics of BSI among children with haematological malignancies (HM) and solid tumours (ST). The study was conducted from October 2012 through to November 2015 at a referral hospital for cancer care and included the first BSI episode detected in 210 patients aged 18 years or less. BSI cases were prospectively detected by daily laboratory-based surveillance. The Centers for Disease Control and Prevention definitions for primary or secondary BSI were used. A higher proportion of use of corticosteroids (P = 0.02), chemotherapy (P = 0.01) and antibiotics (P = 0.05) before the BSI diagnosis; as well as of neutropenia (P < 0.001) and mucositis (P < 0.001) at the time of BSI diagnosis was observed in patients with HM than with ST. Previous surgical procedures (P = 0.03), mechanical ventilation (P = 0.01) and bed confinement (P < 0.001) were more frequent among children with ST. The frequency of use of temporary (P = 0.01) and implanted vascular lines (P < 0.01) was significantly higher in children with ST than with HM while the tunnelled line (P = 0.01) use was more frequent in children with HM as compared to ST. Most (n = 181) BSI cases were primary BSI. BSI associated with a tunnelled catheter was more frequent in children with HM (P < 0.01), whereas BSI associated with an implanted (P < 0.01) or temporary central line (P < 0.02) was more common in patients with ST. BSI associated with mucosal barrier injury was more frequent (P = 0.01) in children with HM. Indication for intensive care was more frequent in children (P = 0.05) with ST. Mortality ratio was similar in children with ST and HM, and length of hospital stay after BSI was higher in patients with HM than with ST (median of 19 vs. 13 days; P = 0.02). Infection caused by Gram-negative bacteria (P = 0.04) and polymicrobial infections (P = 0.05) due to Gram-positive cocci plus fungus was more common in patients with HM. These findings suggest that the characteristics of BSI acquisition and mortality can be cancer-specific.

Richness of human gut microbiome correlates with metabolic markers.

We are facing a global metabolic health crisis provoked by an obesity epidemic. Here we report the human gut microbial composition in a population sample of 123 non-obese and 169 obese Danish individuals. We find two groups of individuals that differ by the number of gut microbial genes and thus gut bacterial richness. They contain known and previously unknown bacterial species at different proportions; individuals with a low bacterial richness (23% of the population) are characterized by more marked overall adiposity, insulin resistance and dyslipidaemia and a more pronounced inflammatory phenotype when compared with high bacterial richness individuals. The obese individuals among the lower bacterial richness group also gain more weight over time. Only a few bacterial species are sufficient to distinguish between individuals with high and low bacterial richness, and even between lean and obese participants. Our classifications based on variation in the gut microbiome identify subsets of individuals in the general white adult population who may be at increased risk of progressing to adiposity-associated co-morbidities.

Subliminal and supraliminal processing of reward-related stimuli in anorexia nervosa.

BACKGROUND: Previous studies have highlighted the role of the brain reward and cognitive control systems in the etiology of anorexia nervosa (AN). In an attempt to disentangle the relative contribution of these systems to the disorder, we used functional magnetic resonance imaging (fMRI) to investigate hemodynamic responses to reward-related stimuli presented both subliminally and supraliminally in acutely underweight AN patients and age-matched healthy controls (HC). METHODS: fMRI data were collected from a total of 35 AN patients and 35 HC, while they passively viewed subliminally and supraliminally presented streams of food, positive social, and neutral stimuli. Activation patterns of the group × stimulation condition × stimulus type interaction were interrogated to investigate potential group differences in processing different stimulus types under the two stimulation conditions. Moreover, changes in functional connectivity were investigated using generalized psychophysiological interaction analysis. RESULTS: AN patients showed a generally increased response to supraliminally presented stimuli in the inferior frontal junction (IFJ), but no alterations within the reward system. Increased activation during supraliminal stimulation with food stimuli was observed in the AN group in visual regions including superior occipital gyrus and the fusiform gyrus/parahippocampal gyrus. No group difference was found with respect to the subliminal stimulation condition and functional connectivity. CONCLUSION: Increased IFJ activation in AN during supraliminal stimulation may indicate hyperactive cognitive control, which resonates with clinical presentation of excessive self-control in AN patients. Increased activation to food stimuli in visual regions may be interpreted in light of an attentional food bias in AN.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index.

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10-5, Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10-06/Pfemales: 3.45 × 10-07/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium.

The profile of brain structural abnormalities in schizophrenia is still not fully understood, despite decades of research using brain scans. To validate a prospective meta-analysis approach to analyzing multicenter neuroimaging data, we analyzed brain MRI scans from 2028 schizophrenia patients and 2540 healthy controls, assessed with standardized methods at 15 centers worldwide. We identified subcortical brain volumes that differentiated patients from controls, and ranked them according to their effect sizes. Compared with healthy controls, patients with schizophrenia had smaller hippocampus (Cohen's d=-0.46), amygdala (d=-0.31), thalamus (d=-0.31), accumbens (d=-0.25) and intracranial volumes (d=-0.12), as well as larger pallidum (d=0.21) and lateral ventricle volumes (d=0.37). Putamen and pallidum volume augmentations were positively associated with duration of illness and hippocampal deficits scaled with the proportion of unmedicated patients. Worldwide cooperative analyses of brain imaging data support a profile of subcortical abnormalities in schizophrenia, which is consistent with that based on traditional meta-analytic approaches. This first ENIGMA Schizophrenia Working Group study validates that collaborative data analyses can readily be used across brain phenotypes and disorders and encourages analysis and data sharing efforts to further our understanding of severe mental illness.

Metagenomic species profiling using universal phylogenetic marker genes.

To quantify known and unknown microorganisms at species-level resolution using shotgun sequencing data, we developed a method that establishes metagenomic operational taxonomic units (mOTUs) based on single-copy phylogenetic marker genes. Applied to 252 human fecal samples, the method revealed that on average 43% of the species abundance and 58% of the richness cannot be captured by current reference genome-based methods. An implementation of the method is available at http://www.bork.embl.de/software/mOTU/.

Prepregnancy Weight in Women with Type I Diabetes Mellitus: Effect on Pregnancy Outcomes.

Objective This study aims to evaluate the association between prepregnancy body mass index (BMI) and adverse pregnancy outcomes in women with type 1 diabetes mellitus (DM). Methods This is a secondary analysis of a cohort of 426 pregnancies in women with type 1 DM recruited before 20 weeks gestation. Women were categorized according to prepregnancy BMI: low BMI (< 20 kg/m2), normal BMI (20 to < 25 kg/m2), and high BMI (≥ 25 kg/m2). The outcomes of interest were: spontaneous abortion (delivery < 20 weeks gestation); preeclampsia; emergent delivery for maternal indications (hypertension or placental abruption); and preterm delivery (< 37 weeks gestation). Analyses included proportional hazards and multiple logistic regression models with covariates: age, age at diagnosis of type 1 DM, previous spontaneous abortion, microvascular disease (nephropathy or retinopathy), and glycohemoglobin A1 concentrations. Results Low BMI was associated with preterm delivery. High BMI was associated with emergent delivery for maternal indications. Glycemic control as measured by glycohemoglobin A1 was associated with increased risk of spontaneous abortion, attenuating the association with low prepregnancy weight. Conclusion Prepregnancy BMI is a risk factor to be considered when caring for women with type 1 DM, in particular for preterm delivery (low BMI) and emergent delivery for maternal indications (high BMI).

Quality control of microbiota metagenomics by k-mer analysis.

BACKGROUND: The biological and clinical consequences of the tight interactions between host and microbiota are rapidly being unraveled by next generation sequencing technologies and sophisticated bioinformatics, also referred to as microbiota metagenomics. The recent success of metagenomics has created a demand to rapidly apply the technology to large case-control cohort studies and to studies of microbiota from various habitats, including habitats relatively poor in microbes. It is therefore of foremost importance to enable a robust and rapid quality assessment of metagenomic data from samples that challenge present technological limits (sample numbers and size). Here we demonstrate that the distribution of overlapping k-mers of metagenome sequence data predicts sequence quality as defined by gene distribution and efficiency of sequence mapping to a reference gene catalogue. RESULTS: We used serial dilutions of gut microbiota metagenomic datasets to generate well-defined high to low quality metagenomes. We also analyzed a collection of 52 microbiota-derived metagenomes. We demonstrate that k-mer distributions of metagenomic sequence data identify sequence contaminations, such as sequences derived from "empty" ligation products. Of note, k-mer distributions were also able to predict the frequency of sequences mapping to a reference gene catalogue not only for the well-defined serial dilution datasets, but also for 52 human gut microbiota derived metagenomic datasets. CONCLUSIONS: We propose that k-mer analysis of raw metagenome sequence reads should be implemented as a first quality assessment prior to more extensive bioinformatics analysis, such as sequence filtering and gene mapping. With the rising demand for metagenomic analysis of microbiota it is crucial to provide tools for rapid and efficient decision making. This will eventually lead to a faster turn-around time, improved analytical quality including sample quality metrics and a significant cost reduction. Finally, improved quality assessment will have a major impact on the robustness of biological and clinical conclusions drawn from metagenomic studies.