Evaluation of an Adapted Semi-Automated DNA Extraction for Human Salivary Shotgun Metagenomics

Recent attention has highlighted the importance of oral microbiota in human health and disease, e.g., in Parkinson's disease, notably using shotgun metagenomics. One key aspect for efficient shotgun metagenomic analysis relies on optimal microbial sampling and DNA extraction, generally implementing commercial solutions developed to improve sample collection and preservation, and provide high DNA quality and quantity for downstream analysis. As metagenomic studies are today performed on a large number of samples, the next evolution to increase study throughput is with DNA extraction automation. In this study, we proposed a semi-automated DNA extraction protocol for human salivary samples collected with a commercial kit, and compared the outcomes with the DNA extraction recommended by the manufacturer. While similar DNA yields were observed between the protocols, our semi-automated DNA protocol generated significantly higher DNA fragment sizes. Moreover, we showed that the oral microbiome composition was equivalent between DNA extraction methods, even at the species level. This study demonstrates that our semi-automated protocol is suitable for shotgun metagenomic analysis, while allowing for improved sample treatment logistics with reduced technical variability and without compromising the structure of the oral microbiome.

Evidence of a causal and modifiable relationship between kidney function and circulating trimethylamine N-oxide

The host-microbiota co-metabolite trimethylamine N-oxide (TMAO) is linked to increased cardiovascular risk but how its circulating levels are regulated remains unclear. We applied "explainable" machine learning, univariate, multivariate and mediation analyses of fasting plasma TMAO concentration and a multitude of phenotypes in 1,741 adult Europeans of the MetaCardis study. Here we show that next to age, kidney function is the primary variable predicting circulating TMAO, with microbiota composition and diet playing minor, albeit significant, roles. Mediation analysis suggests a causal relationship between TMAO and kidney function that we corroborate in preclinical models where TMAO exposure increases kidney scarring. Consistent with our findings, patients receiving glucose-lowering drugs with reno-protective properties have significantly lower circulating TMAO when compared to propensity-score matched control individuals. Our analyses uncover a bidirectional relationship between kidney function and TMAO that can potentially be modified by reno-protective anti-diabetic drugs and suggest a clinically actionable intervention for decreasing TMAO-associated excess cardiovascular risk.

The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice

Anorexia nervosa (AN) is an eating disorder with a high mortality. About 95% of cases are women and it has a population prevalence of about 1%, but evidence-based treatment is lacking. The pathogenesis of AN probably involves genetics and various environmental factors, and an altered gut microbiota has been observed in individuals with AN using amplicon sequencing and relatively small cohorts. Here we investigated whether a disrupted gut microbiota contributes to AN pathogenesis. Shotgun metagenomics and metabolomics were performed on faecal and serum samples, respectively, from a cohort of 77 females with AN and 70 healthy females. Multiple bacterial taxa (for example, Clostridium species) were altered in AN and correlated with estimates of eating behaviour and mental health. The gut virome was also altered in AN including a reduction in viral-bacterial interactions. Bacterial functional modules associated with the degradation of neurotransmitters were enriched in AN and various structural variants in bacteria were linked to metabolic features of AN. Serum metabolomics revealed an increase in metabolites associated with reduced food intake (for example, indole-3-propionic acid). Causal inference analyses implied that serum bacterial metabolites are potentially mediating the impact of an altered gut microbiota on AN behaviour. Further, we performed faecal microbiota transplantation from AN cases to germ-free mice under energy-restricted feeding to mirror AN eating behaviour. We found that the reduced weight gain and induced hypothalamic and adipose tissue gene expression were related to aberrant energy metabolism and eating behaviour. Our 'omics' and mechanistic studies imply that a disruptive gut microbiome may contribute to AN pathogenesis.

In Vitro Modelling of Oral Microbial Invasion in the Human Colon

Recent advances in the human microbiome characterization have revealed significant oral microbial detection in stools of dysbiotic patients. However, little is known about the potential interactions of these invasive oral microorganisms with commensal intestinal microbiota and the host. In this proof-of-concept study, we proposed a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Oral invasion of the intestinal microbiota was simulated by injection of enriched saliva in the in vitro colon model inoculated with a fecal sample from the same healthy adult donor. The mucosal compartment of M-ARCOL was able to retain the highest species richness levels over time, while species richness levels decreased in the luminal compartment. This study also showed that oral microorganisms preferably colonized the mucosal microenvironment, suggesting potential oral-to-intestinal mucosal competitions. This new model of oral-to-gut invasion can provide useful mechanistic insights into the role of oral microbiome in various disease processes. IMPORTANCE Here, we propose a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Our study revealed the importance of integrating the mucus compartment, which retained higher microbial richness during fermentation, showed the preference of oral microbial invaders for the mucosal resources, and indicated potential oral-to-intestinal mucosal competitions. It also underlined promising opportunities to further understand mechanisms of oral invasion into the human gut microbiome, define microbe-microbe and mucus-microbe interactions in a compartmentalized fashion, and help to better characterize the potential of oral microbial invasion and their persistence in the gut.

Metagenomic and metabolomic remodeling in nonagenarians and centenarians and its association with genetic and socioeconomic factors

A better understanding of the biological and environmental variables that contribute to exceptional longevity has the potential to inform the treatment of geriatric diseases and help achieve healthy aging. Here, we compared the gut microbiome and blood metabolome of extremely long-lived individuals (94-105 years old) to that of their children (50-79 years old) in 116 Han Chinese families. We found extensive metagenomic and metabolomic remodeling in advanced age and observed a generational divergence in the correlations with socioeconomic factors. An analysis of quantitative trait loci revealed that genetic associations with metagenomic and metabolomic features were largely generation-specific, but we also found 131 plasma metabolic quantitative trait loci associations that were cross-generational with the genetic variants concentrated in six loci. These included associations between FADS1/2 and arachidonate, PTPA and succinylcarnitine and FLVCR1 and choline. Our characterization of the extensive metagenomic and metabolomic remodeling that occurs in people reaching extreme ages may offer new targets for aging-related interventions.

Metabolomic and microbiome profiling reveals personalized risk factors for coronary artery disease

Complex diseases, such as coronary artery disease (CAD), are often multifactorial, caused by multiple underlying pathological mechanisms. Here, to study the multifactorial nature of CAD, we performed comprehensive clinical and multi-omic profiling, including serum metabolomics and gut microbiome data, for 199 patients with acute coronary syndrome (ACS) recruited from two major Israeli hospitals, and validated these results in a geographically distinct cohort. ACS patients had distinct serum metabolome and gut microbial signatures as compared with control individuals, and were depleted in a previously unknown bacterial species of the Clostridiaceae family. This bacterial species was associated with levels of multiple circulating metabolites in control individuals, several of which have previously been linked to an increased risk of CAD. Metabolic deviations in ACS patients were found to be person specific with respect to their potential genetic or environmental origin, and to correlate with clinical parameters and cardiovascular outcomes. Moreover, metabolic aberrations in ACS patients linked to microbiome and diet were also observed to a lesser extent in control individuals with metabolic impairment, suggesting the involvement of these aberrations in earlier dysmetabolic phases preceding clinically overt CAD. Finally, a metabolomics-based model of body mass index (BMI) trained on the non-ACS cohort predicted higher-than-actual BMI when applied to ACS patients, and the excess BMI predictions independently correlated with both diabetes mellitus (DM) and CAD severity, as defined by the number of vessels involved. These results highlight the utility of the serum metabolome in understanding the basis of risk-factor heterogeneity in CAD.

Conjugated C-6 hydroxylated bile acids in serum relate to human metabolic health and gut Clostridia species

Knowledge about in vivo effects of human circulating C-6 hydroxylated bile acids (BAs), also called muricholic acids, is sparse. It is unsettled if the gut microbiome might contribute to their biosynthesis. Here, we measured a range of serum BAs and related them to markers of human metabolic health and the gut microbiome. We examined 283 non-obese and obese Danish adults from the MetaHit study. Fasting concentrations of serum BAs were quantified using ultra-performance liquid chromatography-tandem mass-spectrometry. The gut microbiome was characterized with shotgun metagenomic sequencing and genome-scale metabolic modeling. We find that tauro- and glycohyocholic acid correlated inversely with body mass index (P = 4.1e-03, P = 1.9e-05, respectively), waist circumference (P = 0.017, P = 1.1e-04, respectively), body fat percentage (P = 2.5e-03, P = 2.3e-06, respectively), insulin resistance (P = 0.051, P = 4.6e-4, respectively), fasting concentrations of triglycerides (P = 0.06, P = 9.2e-4, respectively) and leptin (P = 0.067, P = 9.2e-4). Tauro- and glycohyocholic acids, and tauro-a-muricholic acid were directly linked with a distinct gut microbial community primarily composed of Clostridia species (P = 0.037, P = 0.013, P = 0.027, respectively). We conclude that serum conjugated C-6-hydroxylated BAs associate with measures of human metabolic health and gut communities of Clostridia species. The findings merit preclinical interventions and human feasibility studies to explore the therapeutic potential of these BAs in obesity and type 2 diabetes.

Imidazole propionate is increased in diabetes and associated with dietary patterns and altered microbial ecology

Microbiota-host-diet interactions contribute to the development of metabolic diseases. Imidazole propionate is a novel microbially produced metabolite from histidine, which impairs glucose metabolism. Here, we show that subjects with prediabetes and diabetes in the MetaCardis cohort from three European countries have elevated serum imidazole propionate levels. Furthermore, imidazole propionate levels were increased in subjects with low bacterial gene richness and Bacteroides 2 enterotype, which have previously been associated with obesity. The Bacteroides 2 enterotype was also associated with increased abundance of the genes involved in imidazole propionate biosynthesis from dietary histidine. Since patients and controls did not differ in their histidine dietary intake, the elevated levels of imidazole propionate in type 2 diabetes likely reflects altered microbial metabolism of histidine, rather than histidine intake per se. Thus the microbiota may contribute to type 2 diabetes by generating imidazole propionate that can modulate host inflammation and metabolism.

Gut microbial metabolism of nutrients contributes to metabolic diseases, and the histidine metabolite imidazole propionate (ImP) is produced by type 2 diabetes (T2D) associated microbiome. Here the authors report that circulating ImP levels are increased in subjects with prediabetes or T2D in three European populations, and this increase associates with altered gut microbiota rather than dietary histidine.

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.

Impact of laparoscopic Roux-en-Y gastric bypass and sleeve gastrectomy on gut microbiota: a metagenomic comparative analysis

BACKGROUND

Bariatric surgery is an effective therapeutic procedure for morbidly obese patients. The 2 most common interventions are sleeve gastrectomy (SG) and laparoscopic Roux-en-Y gastric bypass (LRYGB).

OBJECTIVES

The aim of this study was to compare microbiome long-term microbiome after SG and LRYGB surgery in obese patients.

SETTING

University Hospital, France; University Hospital, United States; and University Hospital, Switzerland.

METHODS

Eighty-nine and 108 patients who underwent SG and LRYGB, respectively, were recruited. Stools were collected before and 6 months after surgery. Microbial DNA was analyzed with shotgun metagenomic sequencing (SOLiD 5500 xl Wildfire). MSPminer, a novel innovative tool to characterize new in silico biological entities, was used to identify 715 Metagenomic Species Pan-genome. One hundred forty-eight functional modules were analyzed using GOmixer and KEGG database.

RESULTS

Both interventions resulted in a similar increase of Shannon's diversity index and gene richness of gut microbiota, in parallel with weight loss, but the changes of microbial composition were different. LRYGB led to higher relative abundance of aero-tolerant bacteria, such as Escherichia coli and buccal species, such as Streptococcus and Veillonella spp. In contrast, anaerobes, such as Clostridium, were more abundant after SG, suggesting better conservation of anaerobic conditions in the gut. Enrichment of Akkermansia muciniphila was also observed after both surgeries. Function-level changes included higher potential for bacterial use of supplements, such as vitamin B12, B1, and iron upon LRYGB.

CONCLUSION

Microbiota changes after bariatric surgery depend on the nature of the intervention. LRYGB induces greater taxonomic and functional changes in gut microbiota than SG. Possible long-term health consequences of these alterations remain to be established.

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.

Microbial ecology perturbation in human IgA deficiency

Paradoxically, loss of immunoglobulin A (IgA), one of the most abundant antibodies, does not irrevocably lead to severe infections in humans but rather is associated with relatively mild respiratory infections, atopy, and autoimmunity. IgA might therefore also play covert roles, not uniquely associated with control of pathogens. We show that human IgA deficiency is not associated with massive quantitative perturbations of gut microbial ecology. Metagenomic analysis highlights an expected pathobiont expansion but a less expected depletion in some typically beneficial symbionts. Gut colonization by species usually present in the oropharynx is also reminiscent of spatial microbiota disorganization. IgM only partially rescues IgA deficiency because not all typical IgA targets are efficiently bound by IgM in the intestinal lumen. Together, IgA appears to play a nonredundant role at the forefront of the immune/microbial interface, away from the intestinal barrier, ranging from pathobiont control and regulation of systemic inflammation to preservation of commensal diversity and community networks.

Prediction of the intestinal resistome by a three-dimensional structure-based method

The intestinal microbiota is considered to be a major reservoir of antibiotic resistance determinants (ARDs) that could potentially be transferred to bacterial pathogens via mobile genetic elements. Yet, this assumption is poorly supported by empirical evidence due to the distant homologies between known ARDs (mostly from culturable bacteria) and ARDs from the intestinal microbiota. Consequently, an accurate census of intestinal ARDs (that is, the intestinal resistome) has not yet been fully determined. For this purpose, we developed and validated an annotation method (called pairwise comparative modelling) on the basis of a three-dimensional structure (homology comparative modelling), leading to the prediction of 6,095 ARDs in a catalogue of 3.9 million proteins from the human intestinal microbiota. We found that the majority of predicted ARDs (pdARDs) were distantly related to known ARDs (mean amino acid identity 29.8%) and found little evidence supporting their transfer between species. According to the composition of their resistome, we were able to cluster subjects from the MetaHIT cohort (n = 663) into six resistotypes that were connected to the previously described enterotypes. Finally, we found that the relative abundance of pdARDs was positively associated with gene richness, but not when subjects were exposed to antibiotics. Altogether, our results indicate that the majority of intestinal microbiota ARDs can be considered intrinsic to the dominant commensal microbiota and that these genes are rarely shared with bacterial pathogens.

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.

Transcriptional interactions suggest niche segregation among microorganisms in the human gut

The human gastrointestinal (GI) tract is the habitat for hundreds of microbial species, of which many cannot be cultivated readily, presumably because of the dependencies between species¹. Studies of microbial co-occurrence in the gut have indicated community substructures that may reflect functional and metabolic interactions between cohabiting species^2,3. To move beyond species co-occurrence networks, we systematically identified transcriptional interactions between pairs of coexisting gut microbes using metagenomics and microarray-based metatranscriptomics data from 233 stool samples from Europeans. In 102 significantly interacting species pairs, the transcriptional changes led to a reduced expression of orthologous functions between the coexisting species. Specific species-species transcriptional interactions were enriched for functions important for H₂ and CO₂ homeostasis, butyrate biosynthesis, ATP-binding cassette (ABC) transporters, flagella assembly and bacterial chemotaxis, as well as for the metabolism of carbohydrates, amino acids and cofactors. The analysis gives the first insight into the microbial community-wide transcriptional interactions, and suggests that the regulation of gene expression plays an important role in species adaptation to coexistence and that niche segregation takes place at the transcriptional level.

Tongue Coating and the Salivary Microbial Communities Vary in Children with Halitosis

Halitosis is a common symptom mainly caused by microbial activities in the oral cavity. Here, we used 16S rRNA gene pyrosequencing and metagenomic sequencing to examine oral microbial compositions and their functional variations in children with halitosis. We found that the tongue coating of subjects with halitosis had greater bacterial richness than those of healthy subjects. The relative abundance and prevalence of Leptotrichia wadei and Peptostreptococcus stomatis were higher in tongue coating samples from children with halitosis than those from children without halitosis; Prevotella shahii had higher relative abundance and prevalence in saliva samples from children with halitosis. We present the first comprehensive evaluation of the co-occurrence networks of saliva and tongue coating communities under healthy and halitosis conditions, and investigated patterns of significant differences between these communities. Moreover, we observed that bacterial genes associated with responses to infectious diseases and terpenoid and polyketide metabolism were enriched in subjects with halitosis, but not in healthy subjects. Hydrogen sulphide (H₂S)-related metabolic pathways suggested that there was higher microbial production and less usage of H₂S in subjects with halitosis. Thus, the mechanism of halitosis was implied for the first time via metagenomic sequencing.

Influence of H7N9 virus infection and associated treatment on human gut microbiota

Between March and June, 2013, forty H7N9 patients were hospitalized in our hospital. Next-generation sequencing technologies have been used to sequence the fecal DNA samples of the patient, the within sample diversity analysis, enterotyping, functional gene and metagenomic species analysis have been carried on both the patients and healthy controls. The influence of associated treatment in H7N9 infected patients is dramatic and was firstly revealed in species level due to deep sequencing technology. We found that most of the MetaGenomic Species (MGS) enriched in the control samples were Roseburia inulinivorans DSM 16841, butyrate producing bacterium SS3/4 and most of MGS enriched in the H7N9 patients were Clostridium sp. 7 2 43FAA and Enterococcus faecium. It was concluded that H7N9 viral infection and antibiotic administration have a significant effect on the microbiota community with decreased diversity and overgrowth of the bacteria such as Escherichia coli and Enterococcus faecium. Enterotype analysis showed that the communities were unstable. Treatment including antivirals, probiotics and antibiotics helps to improve the microbiota diversity and the abundance of beneficial bacteria in the gut.

A catalog of the mouse gut metagenome

We established a catalog of the mouse gut metagenome comprising ∼2.6 million nonredundant genes by sequencing DNA from fecal samples of 184 mice. To secure high microbiome diversity, we used mouse strains of diverse genetic backgrounds, from different providers, kept in different housing laboratories and fed either a low-fat or high-fat diet. Similar to the human gut microbiome, >99% of the cataloged genes are bacterial. We identified 541 metagenomic species and defined a core set of 26 metagenomic species found in 95% of the mice. The mouse gut microbiome is functionally similar to its human counterpart, with 95.2% of its Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologous groups in common. However, only 4.0% of the mouse gut microbial genes were shared (95% identity, 90% coverage) with those of the human gut microbiome. This catalog provides a useful reference for future studies.

Specific gut microbiota features and metabolic markers in postmenopausal women with obesity

BACKGROUND

Gut microbial gene richness and specific bacterial species are associated with metabolic risk markers in humans, but the impact of host physiology and dietary habits on the link between the gut microbiota and metabolic markers remain unclear. The objective of this study was to identify gut metagenomic markers associated with estimates of insulin resistance, lipid metabolism and inflammation in obesity, and to explore whether the associations between metagenomic and metabolic markers persisted after adjustment for body fat, age and habitual dietary intake.

METHODS

Faecal DNA from 53 women with obesity was analysed through quantitative metagenomic sequencing and analysis, and a systematic search was performed for bacterial genes associated with estimates of insulin resistance, inflammation and lipid metabolism. Subsequently, the correlations between metagenomic species and metabolic markers were tested by linear regression models, with and without covariate adjustment.

RESULTS

One hundred and fourteen metagenomic species correlated with metabolic markers (P<0.001) including Akkermansia muciniphila, Bilophila wadsworthia, Bifidobacterium longum and Faecalibacterium prausnitzii, but also species not previously associated with metabolic markers including Bacteroides faecis and Dorea longicatena. The majority of the identified correlations between bacterial species and metabolic markers persisted after adjustment for differences in body fat, age and dietary macronutrient composition; however, the negative correlation with insulin resistance observed for B. longum and F. prausnitzii appeared to be modified by the intake of dietary fibre and fat, respectively.

CONCLUSIONS

This study shows that several gut bacterial species are linked to metabolic risk markers in obesity, also after adjustment for potential confounders, such as long-term diet composition. The study supports the use of gut metagenomic markers for metabolic disease prediction and warrants further investigation of causality.

Effect of genome and environment on metabolic and inflammatory profiles

Twin and family studies have established the contribution of genetic factors to variation in metabolic, hematologic and immunological parameters. The majority of these studies analyzed single or combined traits into pre-defined syndromes. In the present study, we explore an alternative multivariate approach in which a broad range of metabolic, hematologic, and immunological traits are analyzed simultaneously to determine the resemblance of monozygotic (MZ) twin pairs, twin-spouse pairs and unrelated, non-cohabiting individuals. A total of 517 participants from the Netherlands Twin Register, including 210 MZ twin pairs and 64 twin-spouse pairs, took part in the study. Data were collected on body composition, blood pressure, heart rate, and multiple biomarkers assessed in fasting blood samples, including lipid levels, glucose, insulin, liver enzymes, hematological measurements and cytokine levels. For all 51 measured traits, pair-wise Pearson correlations, correcting for family relatedness, were calculated across all the individuals in the cohort. Hierarchical clustering techniques were applied to group the measured traits into sub-clusters based on similarity. Sub-clusters were observed among metabolic traits and among inflammatory markers. We defined a phenotypic profile as the collection of all the traits measured for a given individual. Average within-pair similarity of phenotypic profiles was determined for the groups of MZ twin pairs, spouse pairs and pairs of unrelated individuals. The average similarity across the full phenotypic profile was higher for MZ twin pairs than for spouse pairs, and lowest for pairs of unrelated individuals. Cohabiting MZ twins were more similar in their phenotypic profile compared to MZ twins who no longer lived together. The correspondence in the phenotypic profile is therefore determined to a large degree by familial, mostly genetic, factors, while household factors contribute to a lesser degree to profile similarity.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1406-7) contains supplementary material, which is available to authorized users.

Bacterial protein signals are associated with Crohn's disease

OBJECTIVE

No Crohn's disease (CD) molecular maker has advanced to clinical use, and independent lines of evidence support a central role of the gut microbial community in CD. Here we explore the feasibility of extracting bacterial protein signals relevant to CD, by interrogating myriads of intestinal bacterial proteomes from a small number of patients and healthy controls.

DESIGN

We first developed and validated a workflow-including extraction of microbial communities, two-dimensional difference gel electrophoresis (2D-DIGE), and LC-MS/MS-to discover protein signals from CD-associated gut microbial communities. Then we used selected reaction monitoring (SRM) to confirm a set of candidates. In parallel, we used 16S rRNA gene sequencing for an integrated analysis of gut ecosystem structure and functions.

RESULTS

Our 2D-DIGE-based discovery approach revealed an imbalance of intestinal bacterial functions in CD. Many proteins, largely derived from Bacteroides species, were over-represented, while under-represented proteins were mostly from Firmicutes and some Prevotella members. Most overabundant proteins could be confirmed using SRM. They correspond to functions allowing opportunistic pathogens to colonise the mucus layers, breach the host barriers and invade the mucosae, which could still be aggravated by decreased host-derived pancreatic zymogen granule membrane protein GP2 in CD patients. Moreover, although the abundance of most protein groups reflected that of related bacterial populations, we found a specific independent regulation of bacteria-derived cell envelope proteins.

CONCLUSIONS

This study provides the first evidence that quantifiable bacterial protein signals are associated with CD, which can have a profound impact on future molecular diagnosis.

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.

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/.

Building the repertoire of dispensable chromosome regions in Bacillus subtilis entails major refinement of cognate large-scale metabolic model

The nonessential regions in bacterial chromosomes are ill-defined due to incomplete functional information. Here, we establish a comprehensive repertoire of the genome regions that are dispensable for growth of Bacillus subtilis in a variety of media conditions. In complex medium, we attempted deletion of 157 individual regions ranging in size from 2 to 159 kb. A total of 146 deletions were successful in complex medium, whereas the remaining regions were subdivided to identify new essential genes (4) and coessential gene sets (7). Overall, our repertoire covers ∼76% of the genome. We screened for viability of mutant strains in rich defined medium and glucose minimal media. Experimental observations were compared with predictions by the iBsu1103 model, revealing discrepancies that led to numerous model changes, including the large-scale application of model reconciliation techniques. We ultimately produced the iBsu1103V2 model and generated predictions of metabolites that could restore the growth of unviable strains. These predictions were experimentally tested and demonstrated to be correct for 27 strains, validating the refinements made to the model. The iBsu1103V2 model has improved considerably at predicting loss of viability, and many insights gained from the model revisions have been integrated into the Model SEED to improve reconstruction of other microbial models.

Enterotypes of the human gut microbiome

Our knowledge of species and functional composition of the human gut microbiome is rapidly increasing, but it is still based on very few cohorts and little is known about variation across the world. By combining 22 newly sequenced faecal metagenomes of individuals from four countries with previously published data sets, here we identify three robust clusters (referred to as enterotypes hereafter) that are not nation or continent specific. We also confirmed the enterotypes in two published, larger cohorts, indicating that intestinal microbiota variation is generally stratified, not continuous. This indicates further the existence of a limited number of well-balanced host-microbial symbiotic states that might respond differently to diet and drug intake. The enterotypes are mostly driven by species composition, but abundant molecular functions are not necessarily provided by abundant species, highlighting the importance of a functional analysis to understand microbial communities. Although individual host properties such as body mass index, age, or gender cannot explain the observed enterotypes, data-driven marker genes or functional modules can be identified for each of these host properties. For example, twelve genes significantly correlate with age and three functional modules with the body mass index, hinting at a diagnostic potential of microbial markers.

[Metagenomics of the intestinal microbiota: potential applications]

A major challenge in the human metagenomics field is to identify associations of the bacterial genes and human phenotypes and act to modulate microbial populations in order to improve human health and wellbeing. MetaHIT project addresses this ambitious challenge by developing and integrating a number of necessary approaches within the context of the gut microbiome. Among the first results is the establishment of a broad catalog of the human gut microbial genes, which was achieved by an original application of the new generation sequencing technology. The catalog contains 3.3 million non-redundant genes, 150-fold more than the human genome equivalent and includes a large majority of the gut metagenomic sequences determined across three continents, Europe, America and Asia. Its content corresponds to some 1000 bacterial species, which likely represent a large fraction of species associated with humankind intestinal tract. The catalog enables development of the gene profiling approaches aiming to detect associations of bacterial genes and phenotypes. These should lead to the speedy development of diagnostic and prognostic tools and open avenues to reasoned approaches to the modulation of the individual's microbiota in order to optimize health and well-being.

Functional metagenomics: a high throughput screening method to decipher microbiota-driven NF-κB modulation in the human gut

BACKGROUND/AIM

The human intestinal microbiota plays an important role in modulation of mucosal immune responses. To study interactions between intestinal epithelial cells (IECs) and commensal bacteria, a functional metagenomic approach was developed. One interest of metagenomics is to provide access to genomes of uncultured microbes. We aimed at identifying bacterial genes involved in regulation of NF-κB signaling in IECs. A high throughput cell-based screening assay allowing rapid detection of NF-κB modulation in IECs was established using the reporter-gene strategy to screen metagenomic libraries issued from the human intestinal microbiota.

METHODS

A plasmid containing the secreted alkaline phosphatase (SEAP) gene under the control of NF-κB binding elements was stably transfected in HT-29 cells. The reporter clone HT-29/kb-seap-25 was selected and characterized. Then, a first screening of a metagenomic library from Crohn's disease patients was performed to identify NF-κB modulating clones. Furthermore, genes potentially involved in the effect of one stimulatory metagenomic clone were determined by sequence analysis associated to mutagenesis by transposition.

RESULTS

The two proinflammatory cytokines, TNF-α and IL-1β, were able to activate the reporter system, translating the activation of the NF-κB signaling pathway and NF-κB inhibitors, BAY 11-7082, caffeic acid phenethyl ester and MG132 were efficient. A screening of 2640 metagenomic clones led to the identification of 171 modulating clones. Among them, one stimulatory metagenomic clone, 52B7, was further characterized. Sequence analysis revealed that its metagenomic DNA insert might belong to a new Bacteroides strain and we identified 2 loci encoding an ABC transport system and a putative lipoprotein potentially involved in 52B7 effect on NF-κB.

CONCLUSIONS

We have established a robust high throughput screening assay for metagenomic libraries derived from the human intestinal microbiota to study bacteria-driven NF-κB regulation. This opens a strategic path toward the identification of bacterial strains and molecular patterns presenting a potential therapeutic interest.

A human gut microbial gene catalogue established by metagenomic sequencing

To understand the impact of gut microbes on human health and well-being it is crucial to assess their genetic potential. Here we describe the Illumina-based metagenomic sequencing, assembly and characterization of 3.3 million non-redundant microbial genes, derived from 576.7 gigabases of sequence, from faecal samples of 124 European individuals. The gene set, approximately 150 times larger than the human gene complement, contains an overwhelming majority of the prevalent (more frequent) microbial genes of the cohort and probably includes a large proportion of the prevalent human intestinal microbial genes. The genes are largely shared among individuals of the cohort. Over 99% of the genes are bacterial, indicating that the entire cohort harbours between 1,000 and 1,150 prevalent bacterial species and each individual at least 160 such species, which are also largely shared. We define and describe the minimal gut metagenome and the minimal gut bacterial genome in terms of functions present in all individuals and most bacteria, respectively.

Activation of mRNA translation by phage protein and low temperature: the case of Lactococcus lactis abortive infection system AbiD1

Background

Abortive infection (Abi) mechanisms comprise numerous strategies developed by bacteria to avoid being killed by bacteriophage (phage). Escherichia coli Abis are considered as mediators of programmed cell death, which is induced by infecting phage. Abis were also proposed to be stress response elements, but no environmental activation signals have yet been identified. Abis are widespread in Lactococcus lactis, but regulation of their expression remains an open question. We previously showed that development of AbiD1 abortive infection against phage bIL66 depends on orf1, which is expressed in mid-infection. However, molecular basis for this activation remains unclear.

Results

In non-infected AbiD1+ cells, specific abiD1 mRNA is unstable and present in low amounts. It does not increase during abortive infection of sensitive phage. Protein synthesis directed by the abiD1 translation initiation region is also inefficient. The presence of the phage orf1 gene, but not its mutant AbiD1^R allele, strongly increases abiD1 translation efficiency. Interestingly, cell growth at low temperature also activates translation of abiD1 mRNA and consequently the AbiD1 phenotype, and occurs independently of phage infection. There is no synergism between the two abiD1 inducers. Purified Orf1 protein binds mRNAs containing a secondary structure motif, identified within the translation initiation regions of abiD1, the mid-infection phage bIL66 M-operon, and the L. lactis osmC gene.

Conclusion

Expression of the abiD1 gene and consequently AbiD1 phenotype is specifically translationally activated by the phage Orf1 protein. The loss of ability to activate translation of abiD1 mRNA determines the molecular basis for phage resistance to AbiD1. We show for the first time that temperature downshift also activates abortive infection by activation of abiD1 mRNA translation.

Anticipating chromosomal replication fork arrest: SSB targets repair DNA helicases to active forks

In bacteria, several salvage responses to DNA replication arrest culminate in reassembly of the replisome on inactivated forks to resume replication. The PriA DNA helicase is a prominent trigger of this replication restart process, preceded in many cases by a repair and/or remodeling of the arrested fork, which can be performed by many specific proteins. The mechanisms that target these rescue effectors to damaged forks in the cell are unknown. We report that the single-stranded DNA binding (SSB) protein is the key factor that links PriA to active chromosomal replication forks in vivo. This targeting mechanism determines the efficiency by which PriA reaches its specific DNA-binding site in vitro and directs replication restart in vivo. The RecG and RecQ DNA helicases, which are involved in intricate replication reactivation pathways, also associate with the chromosomal replication forks by similarly interacting with SSB. These results identify SSB as a platform for linking a 'repair toolbox' with active replication forks, providing a first line of rescue responses to accidental arrest.

Extensive horizontal transfer of core genome genes between two Lactobacillus species found in the gastrointestinal tract

BACKGROUND

While genes that are conserved between related bacterial species are usually thought to have evolved along with the species, phylogenetic trees reconstructed for individual genes may contradict this picture and indicate horizontal gene transfer. Individual trees are often not resolved with high confidence, however, and in that case alternative trees are generally not considered as contradicting the species tree, although not confirming it either. Here we conduct an in-depth analysis of 401 protein phylogenetic trees inferred with varying levels of confidence for three lactobacilli from the acidophilus complex. At present the relationship between these bacteria, isolated from environments as diverse as the gastrointestinal tract (Lactobacillus acidophilus and Lactobacillus johnsonii) and yogurt (Lactobacillus delbrueckii ssp. bulgaricus), is ambiguous due to contradictory phenotypical and 16S rRNA based classifications.

RESULTS

Among the 401 phylogenetic trees, those that could be reconstructed with high confidence support the 16S-rRNA tree or one alternative topology in an astonishing 3:2 ratio, while the third possible topology is practically absent. Lowering the confidence threshold for trees to be taken into consideration does not significantly affect this ratio, and therefore suggests that gene transfer may have affected as much as 40% of the core genome genes. Gene function bias suggests that the 16S rRNA phylogeny of the acidophilus complex, which indicates that L. acidophilus and L. delbrueckii ssp. bulgaricus are the closest related of these three species, is correct. A novel approach of comparison of interspecies protein divergence data employed in this study allowed to determine that gene transfer most likely took place between the lineages of the two species found in the gastrointestinal tract.

CONCLUSION

This case-study reports an unprecedented level of phylogenetic incongruence, presumably resulting from extensive horizontal gene transfer. The data give a first indication of the large extent of gene transfer that may take place in the gastrointestinal tract and its accumulated effect. For future studies, our results should encourage a careful weighing of data on phylogenetic tree topology, confidence and distribution to conclude on the absence or presence and extent of horizontal gene transfer.

KSY1, a lactococcal phage with a T7-like transcription

The virulent lactococcal phage KSY1 possesses a large elongated capsid (223 nm long, 45 nm wide) and a short tail (32 nm). This phage of the Podoviridae group (C3 morphotype) has a linear 79,232-bp double-stranded DNA genome, which encodes 131 putative proteins and 3 tRNAs. This is the first description of the genome of a phage of this morphotype. KSY1 possesses a T7-like transcription system, including an RNA polymerase and a series of specific promoters, showing sequence homology to other known T7-like RNA polymerase promoters. Late stages of KSY1 multiplication are resistant to rifampicin. Otherwise, KSY1 shares limited similarity with other Podoviridae phages. Fourteen KSY1 structural proteins were identified by SDS-PAGE analysis. Among these proteins, those forming the distal tail structure and likely involved in host recognition are encoded by a 5-kb genomic region of KSY1. This region consists of a mosaic of DNA segments highly homologous to DNA of other lactococcal phages, suggesting an horizontal gene transfer.

Genetic evidence for a link between glycolysis and DNA replication

Background

A challenging goal in biology is to understand how the principal cellular functions are integrated so that cells achieve viability and optimal fitness in a wide range of nutritional conditions.

Methodology/Principal Findings

We report here a tight link between glycolysis and DNA synthesis. The link, discovered during an analysis of suppressors of thermosensitive replication mutants in bacterium Bacillus subtilis, is very strong as some metabolic alterations fully restore viability to replication mutants in which a lethal arrest of DNA synthesis otherwise occurs at a high, restrictive, temperature. Full restoration of viability by such alterations was limited to cells with mutations in three elongation factors (the lagging strand DnaE polymerase, the primase and the helicase) out of a large set of thermosensitive mutants affected in most of the replication proteins. Restoration of viability resulted, at least in part, from maintenance of replication protein activity at high temperature. Physiological studies suggested that this restoration depended on the activity of the three-carbon part of the glycolysis/gluconeogenesis pathway and occurred in both glycolytic and gluconeogenic regimens. Restoration took place abruptly over a narrow range of expression of genes in the three-carbon part of glycolysis. However, the absolute value of this range varied greatly with the allele in question. Finally, restoration of cell viability did not appear to be the result of a decrease in growth rate or an induction of major stress responses.

Conclusions/Significance

Our findings provide the first evidence for a genetic system that connects DNA chain elongation to glycolysis. Its role may be to modulate some aspect of DNA synthesis in response to the energy provided by the environment and the underlying mechanism is discussed. It is proposed that related systems are ubiquitous.

Extending the Bacillus cereus group genomics to putative food-borne pathogens of different toxicity

The Bacillus cereus group represents sporulating soil bacteria containing pathogenic strains which may cause diarrheic or emetic food poisoning outbreaks. Multiple locus sequence typing revealed a presence in natural samples of these bacteria of about 30 clonal complexes. Application of genomic methods to this group was however biased due to the major interest for representatives closely related to Bacillus anthracis. Albeit the most important food-borne pathogens were not yet defined, existing data indicate that they are scattered all over the phylogenetic tree. The preliminary analysis of the sequences of three genomes discussed in this paper narrows down the gaps in our knowledge of the B. cereus group. The strain NVH391-98 is a rare but particularly severe food-borne pathogen. Sequencing revealed that the strain should be a representative of a novel bacterial species, for which the name Bacillus cytotoxis or Bacillus cytotoxicus is proposed. This strain has a reduced genome size compared to other B. cereus group strains. Genome analysis revealed absence of sigma B factor and the presence of genes encoding diarrheic Nhe toxin, not detected earlier. The strain B. cereus F837/76 represents a clonal complex close to that of B. anthracis. Including F837/76, three such B. cereus strains had been sequenced. Alignment of genomes suggests that B. anthracis is their common ancestor. Since such strains often emerge from clinical cases, they merit a special attention. The third strain, KBAB4, is a typical facultative psychrophile generally found in soil. Phylogenic studies show that in nature it is the most active group in terms of gene exchange. Genomic sequence revealed high presence of extra-chromosomal genetic material (about 530kb) that may account for this phenomenon. Genes coding Nhe-like toxin were found on a big plasmid in this strain. This may indicate a potential mechanism of toxicity spread from the psychrophile strain community. The results of this genomic work and ecological compartments of different strains incite to consider a necessity of creating prophylactic vaccines against bacteria closely related to NVH391-98 and F837/76. Presumably developing of such vaccines can be based on the properties of non-pathogenic strains such as KBAB4 or ATCC14579 reported here or earlier. By comparing the protein coding genes of strains being sequenced in this project to others we estimate the shared proteome, or core genome, in the B. cereus group to be 3000+/-200 genes and the total proteome, or pan-genome, to be 20-25,000 genes.