High-sensitivity pattern discovery in large, paired multiomic datasets.

MOTIVATION: Modern biological screens yield enormous numbers of measurements, and identifying and interpreting statistically significant associations among features are essential. In experiments featuring multiple high-dimensional datasets collected from the same set of samples, it is useful to identify groups of associated features between the datasets in a way that provides high statistical power and false discovery rate (FDR) control. RESULTS: Here, we present a novel hierarchical framework, HAllA (Hierarchical All-against-All association testing), for structured association discovery between paired high-dimensional datasets. HAllA efficiently integrates hierarchical hypothesis testing with FDR correction to reveal significant linear and non-linear block-wise relationships among continuous and/or categorical data. We optimized and evaluated HAllA using heterogeneous synthetic datasets of known association structure, where HAllA outperformed all-against-all and other block-testing approaches across a range of common similarity measures. We then applied HAllA to a series of real-world multiomics datasets, revealing new associations between gene expression and host immune activity, the microbiome and host transcriptome, metabolomic profiling and human health phenotypes. AVAILABILITY AND IMPLEMENTATION: An open-source implementation of HAllA is freely available at http://huttenhower.sph.harvard.edu/halla along with documentation, demo datasets and a user group. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration.

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Linking Strain Engraftment in Fecal Microbiota Transplantation With Maintenance of Remission in Crohn's Disease.

BACKGROUND & AIMS: Crohn's disease (CD) is a chronic gastrointestinal disease resulting from the dysfunctional interplay between genetic susceptibility, the immune system, and commensal intestinal microbiota. Emerging evidence suggests that treatment by suppression of the immune response and replacement of the microbiota through fecal microbiota transplantation (FMT) is a promising approach for the treatment of CD. METHODS: We obtained stool metagenomes from CD patients in remission and assessed gut microbiome composition before and after FMT at the species and strain levels. Longitudinal follow-up evaluation allowed us to identify the gain, loss, and strain replacement of specific species and link these events to the maintenance of remission in CD. RESULTS: We found that FMT had a significant long-term effect on patient microbial compositions, although this was primarily driven by the engraftment of donor species, which remained at low abundance. Thirty-eight percent of FMT-driven changes were strain replacements, emphasizing the importance of detailed profiling methods, such as metagenomics. Several instances of long-term coexistence between donor and patient strains were also observed. Engraftment of some Actinobacteria, and engraftment or loss of Proteobacteria, were related to better disease outcomes in CD patients who received FMT, and transmission of Bacteroidetes was deleterious. CONCLUSIONS: Our results suggest clades that may be beneficial to transmit/eliminate through FMT, and provide criteria that may help identify personalized FMT donors to more effectively maintain remission in CD patients. The framework established here creates a foundation for future studies centered around the application of FMT and defined microbial communities as a therapeutic approach for treating CD.

Structure of the Mucosal and Stool Microbiome in Lynch Syndrome.

The gut microbiota has been associated with colorectal cancer (CRC), but causal alterations preceding CRC have not been elucidated. To prospectively assess microbiome changes prior to colorectal neoplasia, we investigated samples from 100 Lynch syndrome patients using 16S rRNA gene sequencing of colon biopsies, coupled with metagenomic and metatranscriptomic sequencing of feces. Colectomy and CRC history represented the largest effects on microbiome profiles. A subset of Clostridiaceae were depleted in stool corresponding with baseline adenomas, while Desulfovibrio was enriched both in stool and in mucosal biopsies. A classifier leveraging stool metatranscriptomes resulted in modest power to predict interval development of preneoplastic colonic adenoma. Predictive transcripts corresponded with a shift in flagellin contributors and oxidative metabolic microenvironment, potentially factors in local CRC pathogenesis. This suggests that the effectiveness of prospective microbiome monitoring for adenomas may be limited but supports the potential causality of these consistent, early microbial changes in colonic neoplasia.

Global chemical effects of the microbiome include new bile-acid conjugations.

A mosaic of cross-phylum chemical interactions occurs between all metazoans and their microbiomes. A number of molecular families that are known to be produced by the microbiome have a marked effect on the balance between health and disease1-9. Considering the diversity of the human microbiome (which numbers over 40,000 operational taxonomic units10), the effect of the microbiome on the chemistry of an entire animal remains underexplored. Here we use mass spectrometry informatics and data visualization approaches11-13 to provide an assessment of the effects of the microbiome on the chemistry of an entire mammal by comparing metabolomics data from germ-free and specific-pathogen-free mice. We found that the microbiota affects the chemistry of all organs. This included the amino acid conjugations of host bile acids that were used to produce phenylalanocholic acid, tyrosocholic acid and leucocholic acid, which have not previously been characterized despite extensive research on bile-acid chemistry14. These bile-acid conjugates were also found in humans, and were enriched in patients with inflammatory bowel disease or cystic fibrosis. These compounds agonized the farnesoid X receptor in vitro, and mice gavaged with the compounds showed reduced expression of bile-acid synthesis genes in vivo. Further studies are required to confirm whether these compounds have a physiological role in the host, and whether they contribute to gut diseases that are associated with microbiome dysbiosis.

Association Between Sulfur-Metabolizing Bacterial Communities in Stool and Risk of Distal Colorectal Cancer in Men.

BACKGROUND & AIMS: Sulfur-metabolizing microbes, which convert dietary sources of sulfur into genotoxic hydrogen sulfide (H2S), have been associated with development of colorectal cancer (CRC). We identified a dietary pattern associated with sulfur-metabolizing bacteria in stool and then investigated its association with risk of incident CRC using data from a large prospective study of men. METHODS: We collected data from 51,529 men enrolled in the Health Professionals Follow-up Study since 1986 to determine the association between sulfur-metabolizing bacteria in stool and risk of CRC over 26 years of follow-up. First, in a subcohort of 307 healthy men, we profiled serial stool metagenomes and metatranscriptomes and assessed diet using semiquantitative food frequency questionnaires to identify food groups associated with 43 bacterial species involved in sulfur metabolism. We used these data to develop a sulfur microbial dietary score. We then used Cox proportional hazards modeling to evaluate adherence to this pattern among eligible individuals (n = 48,246) from 1986 through 2012 with risk for incident CRC. RESULTS: Foods associated with higher sulfur microbial diet scores included increased consumption of processed meats and low-calorie drinks and lower consumption of vegetables and legumes. Increased sulfur microbial diet scores were associated with risk of distal colon and rectal cancers, after adjusting for other risk factors (multivariable relative risk, highest vs lowest quartile, 1.43; 95% confidence interval 1.14-1.81; P-trend = .002). In contrast, sulfur microbial diet scores were not associated with risk of proximal colon cancer (multivariable relative risk 0.86; 95% CI 0.65-1.14; P-trend = .31). CONCLUSIONS: In an analysis of participants in the Health Professionals Follow-up Study, we found that long-term adherence to a dietary pattern associated with sulfur-metabolizing bacteria in stool was associated with an increased risk of distal CRC. Further studies are needed to determine how sulfur-metabolizing bacteria might contribute to CRC pathogenesis.

Growth effects of N-acylethanolamines on gut bacteria reflect altered bacterial abundances in inflammatory bowel disease.

Inflammatory bowel diseases (IBD) are associated with alterations in gut microbial abundances and lumenal metabolite concentrations, but the effects of specific metabolites on the gut microbiota in health and disease remain largely unknown. Here, we analysed the influences of metabolites that are differentially abundant in IBD on the growth and physiology of gut bacteria that are also differentially abundant in IBD. We found that N-acylethanolamines (NAEs), a class of endogenously produced signalling lipids elevated in the stool of IBD patients and a T-cell transfer model of colitis, stimulated growth of species over-represented in IBD and inhibited that of species depleted in IBD in vitro. Using metagenomic sequencing, we recapitulated the effects of NAEs in complex microbial communities ex vivo, with Proteobacteria blooming and Bacteroidetes declining in the presence of NAEs. Metatranscriptomic analysis of the same communities identified components of the respiratory chain as important for the metabolism of NAEs, and this was verified using a mutant deficient for respiratory complex I. In this study, we identified NAEs as a class of metabolites that are elevated in IBD and have the potential to shift gut microbiota towards an IBD-like composition.

Obese Individuals with and without Type 2 Diabetes Show Different Gut Microbial Functional Capacity and Composition.

Obesity and type 2 diabetes (T2D) are metabolic disorders that are linked to microbiome alterations. However, their co-occurrence poses challenges in disentangling microbial features unique to each condition. We analyzed gut microbiomes of lean non-diabetic (n = 633), obese non-diabetic (n = 494), and obese individuals with T2D (n = 153) from German population and metabolic disease cohorts. Microbial taxonomic and functional profiles were analyzed along with medical histories, serum metabolomics, biometrics, and dietary data. Obesity was associated with alterations in microbiome composition, individual taxa, and functions with notable changes in Akkermansia, Faecalibacterium, Oscillibacter, and Alistipes, as well as in serum metabolites that correlated with gut microbial patterns. However, microbiome associations were modest for T2D, with nominal increases in Escherichia/Shigella. Medications, including antihypertensives and antidiabetics, along with dietary supplements including iron, were significantly associated with microbiome variation. These results differentiate microbial components of these interrelated metabolic diseases and identify dietary and medication exposures to consider in future studies.

Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases.

Inflammatory bowel diseases, which include Crohn's disease and ulcerative colitis, affect several million individuals worldwide. Crohn's disease and ulcerative colitis are complex diseases that are heterogeneous at the clinical, immunological, molecular, genetic, and microbial levels. Individual contributing factors have been the focus of extensive research. As part of the Integrative Human Microbiome Project (HMP2 or iHMP), we followed 132 subjects for one year each to generate integrated longitudinal molecular profiles of host and microbial activity during disease (up to 24 time points each; in total 2,965 stool, biopsy, and blood specimens). Here we present the results, which provide a comprehensive view of functional dysbiosis in the gut microbiome during inflammatory bowel disease activity. We demonstrate a characteristic increase in facultative anaerobes at the expense of obligate anaerobes, as well as molecular disruptions in microbial transcription (for example, among clostridia), metabolite pools (acylcarnitines, bile acids, and short-chain fatty acids), and levels of antibodies in host serum. Periods of disease activity were also marked by increases in temporal variability, with characteristic taxonomic, functional, and biochemical shifts. Finally, integrative analysis identified microbial, biochemical, and host factors central to this dysregulation. The study's infrastructure resources, results, and data, which are available through the Inflammatory Bowel Disease Multi'omics Database ( http://ibdmdb.org ), provide the most comprehensive description to date of host and microbial activities in inflammatory bowel diseases.

Chromosome and plasmid-borne PLacO3O1 promoters differ in sensitivity to critically low temperatures.

Temperature shifts trigger genome-wide changes in Escherichia coli's gene expression. We studied if chromosome integration impacts on a gene's sensitivity to these shifts, by comparing the single-RNA production kinetics of a PLacO3O1 promoter, when chromosomally-integrated and when single-copy plasmid-borne. At suboptimal temperatures their induction range, fold change, and response to decreasing temperatures are similar. At critically low temperatures, the chromosome-integrated promoter becomes weaker and noisier. Dissection of its initiation kinetics reveals longer lasting states preceding open complex formation, suggesting enhanced supercoiling buildup. Measurements with Gyrase and Topoisomerase I inhibitors suggest hindrance to escape supercoiling buildup at low temperatures. Consistently, similar phenomena occur in energy-depleted cells by DNP at 30 °C. Transient, critically-low temperatures have no long-term consequences, as raising temperature quickly restores transcription rates. We conclude that the chromosomally-integrated PLacO3O1 has higher sensitivity to low temperatures, due to longer-lasting super-coiled states. A lesser active, chromosome-integrated native lac is shown to be insensitive to Gyrase overexpression, even at critically low temperatures, indicating that the rate of escaping positive supercoiling buildup is temperature and transcription rate dependent. A genome-wide analysis supports this, since cold-shock genes exhibit atypical supercoiling-sensitivities. This phenomenon might partially explain the temperature-sensitivity of some transcriptional programs of E. coli.

Metatranscriptome of human faecal microbial communities in a cohort of adult men.

The gut microbiome is intimately related to human health, but it is not yet known which functional activities are driven by specific microorganisms' ecological configurations or transcription. We report a large-scale investigation of 372 human faecal metatranscriptomes and 929 metagenomes from a subset of 308 men in the Health Professionals Follow-Up Study. We identified a metatranscriptomic 'core' universally transcribed over time and across participants, often by different microorganisms. In contrast to the housekeeping functions enriched in this core, a 'variable' metatranscriptome included specialized pathways that were differentially expressed both across participants and among microorganisms. Finally, longitudinal metagenomic profiles allowed ecological interaction network reconstruction, which remained stable over the six-month timespan, as did strain tracking within and between participants. These results provide an initial characterization of human faecal microbial ecology into core, subject-specific, microorganism-specific and temporally variable transcription, and they differentiate metagenomically versus metatranscriptomically informative aspects of the human faecal microbiome.

Stability of the human faecal microbiome in a cohort of adult men.

Characterizing the stability of the gut microbiome is important to exploit it as a therapeutic target and diagnostic biomarker. We metagenomically and metatranscriptomically sequenced the faecal microbiomes of 308 participants in the Health Professionals Follow-Up Study. Participants provided four stool samples-one pair collected 24-72 h apart and a second pair ~6 months later. Within-person taxonomic and functional variation was consistently lower than between-person variation over time. In contrast, metatranscriptomic profiles were comparably variable within and between subjects due to higher within-subject longitudinal variation. Metagenomic instability accounted for ~74% of corresponding metatranscriptomic instability. The rest was probably attributable to sources such as regulation. Among the pathways that were differentially regulated, most were consistently over- or under-transcribed at each time point. Together, these results suggest that a single measurement of the faecal microbiome can provide long-term information regarding organismal composition and functional potential, but repeated or short-term measures may be necessary for dynamic features identified by metatranscriptomics.

Dynamics of metatranscription in the inflammatory bowel disease gut microbiome.

Inflammatory bowel disease (IBD) is a group of chronic diseases of the digestive tract that affects millions of people worldwide. Genetic, environmental and microbial factors have been implicated in the onset and exacerbation of IBD. However, the mechanisms associating gut microbial dysbioses and aberrant immune responses remain largely unknown. The integrative Human Microbiome Project seeks to close these gaps by examining the dynamics of microbiome functionality in disease by profiling the gut microbiomes of >100 individuals sampled over a 1-year period. Here, we present the first results based on 78 paired faecal metagenomes and metatranscriptomes, and 222 additional metagenomes from 59 patients with Crohn's disease, 34 with ulcerative colitis and 24 non-IBD control patients. We demonstrate several cases in which measures of microbial gene expression in the inflamed gut can be informative relative to metagenomic profiles of functional potential. First, although many microbial organisms exhibited concordant DNA and RNA abundances, we also detected species-specific biases in transcriptional activity, revealing predominant transcription of pathways by individual microorganisms per host (for example, by Faecalibacterium prausnitzii). Thus, a loss of these organisms in disease may have more far-reaching consequences than suggested by their genomic abundances. Furthermore, we identified organisms that were metagenomically abundant but inactive or dormant in the gut with little or no expression (for example, Dialister invisus). Last, certain disease-specific microbial characteristics were more pronounced or only detectable at the transcript level, such as pathways that were predominantly expressed by different organisms in patients with IBD (for example, Bacteroides vulgatus and Alistipes putredinis). This provides potential insights into gut microbial pathway transcription that can vary over time, inducing phenotypical changes that are complementary to those linked to metagenomic abundances. The study's results highlight the strength of analysing both the activity and the presence of gut microorganisms to provide insight into the role of the microbiome in IBD.

Erratum: Strains, functions and dynamics in the expanded Human Microbiome Project.

This corrects the article DOI: 10.1038/nature23889.

Strains, functions and dynamics in the expanded Human Microbiome Project.

The characterization of baseline microbial and functional diversity in the human microbiome has enabled studies of microbiome-related disease, diversity, biogeography, and molecular function. The National Institutes of Health Human Microbiome Project has provided one of the broadest such characterizations so far. Here we introduce a second wave of data from the study, comprising 1,631 new metagenomes (2,355 total) targeting diverse body sites with multiple time points in 265 individuals. We applied updated profiling and assembly methods to provide new characterizations of microbiome personalization. Strain identification revealed subspecies clades specific to body sites; it also quantified species with phylogenetic diversity under-represented in isolate genomes. Body-wide functional profiling classified pathways into universal, human-enriched, and body site-enriched subsets. Finally, temporal analysis decomposed microbial variation into rapidly variable, moderately variable, and stable subsets. This study furthers our knowledge of baseline human microbial diversity and enables an understanding of personalized microbiome function and dynamics.

sgnesR: An R package for simulating gene expression data from an underlying real gene network structure considering delay parameters.

BACKGROUND: sgnesR (Stochastic Gene Network Expression Simulator in R) is an R package that provides an interface to simulate gene expression data from a given gene network using the stochastic simulation algorithm (SSA). The package allows various options for delay parameters and can easily included in reactions for promoter delay, RNA delay and Protein delay. A user can tune these parameters to model various types of reactions within a cell. As examples, we present two network models to generate expression profiles. We also demonstrated the inference of networks and the evaluation of association measure of edge and non-edge components from the generated expression profiles. RESULTS: The purpose of sgnesR is to enable an easy to use and a quick implementation for generating realistic gene expression data from biologically relevant networks that can be user selected. CONCLUSIONS: sgnesR is freely available for academic use. The R package has been tested for R 3.2.0 under Linux, Windows and Mac OS X.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements.

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the Plac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the PtetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

The healthy human microbiome.

Humans are virtually identical in their genetic makeup, yet the small differences in our DNA give rise to tremendous phenotypic diversity across the human population. By contrast, the metagenome of the human microbiome-the total DNA content of microbes inhabiting our bodies-is quite a bit more variable, with only a third of its constituent genes found in a majority of healthy individuals. Understanding this variability in the "healthy microbiome" has thus been a major challenge in microbiome research, dating back at least to the 1960s, continuing through the Human Microbiome Project and beyond. Cataloguing the necessary and sufficient sets of microbiome features that support health, and the normal ranges of these features in healthy populations, is an essential first step to identifying and correcting microbial configurations that are implicated in disease. Toward this goal, several population-scale studies have documented the ranges and diversity of both taxonomic compositions and functional potentials normally observed in the microbiomes of healthy populations, along with possible driving factors such as geography, diet, and lifestyle. Here, we review several definitions of a 'healthy microbiome' that have emerged, the current understanding of the ranges of healthy microbial diversity, and gaps such as the characterization of molecular function and the development of ecological therapies to be addressed in the future.

Dissecting the stochastic transcription initiation process in live Escherichia coli.

We investigate the hypothesis that, in Escherichia coli, while the concentration of RNA polymerases differs in different growth conditions, the fraction of RNA polymerases free for transcription remains approximately constant within a certain range of these conditions. After establishing this, we apply a standard model-fitting procedure to fully characterize the in vivo kinetics of the rate-limiting steps in transcription initiation of the Plac/ara-1 promoter from distributions of intervals between transcription events in cells with different RNA polymerase concentrations. We find that, under full induction, the closed complex lasts ∼788 s while subsequent steps last ∼193 s, on average. We then establish that the closed complex formation usually occurs multiple times prior to each successful initiation event. Furthermore, the promoter intermittently switches to an inactive state that, on average, lasts ∼87 s. This is shown to arise from the intermittent repression of the promoter by LacI. The methods employed here should be of use to resolve the rate-limiting steps governing the in vivo dynamics of initiation of prokaryotic promoters, similar to established steady-state assays to resolve the in vitro dynamics.

Increased cytoplasm viscosity hampers aggregate polar segregation in Escherichia coli.

In Escherichia coli, under optimal conditions, protein aggregates associated with cellular aging are excluded from midcell by the nucleoid. We study the functionality of this process under sub-optimal temperatures from population and time lapse images of individual cells and aggregates and nucleoids within. We show that, as temperature decreases, aggregates become homogeneously distributed and uncorrelated with nucleoid size and location. We present evidence that this is due to increased cytoplasm viscosity, which weakens the anisotropy in aggregate displacements at the nucleoid borders that is responsible for their preference for polar localisation. Next, we show that in plasmolysed cells, which have increased cytoplasm viscosity, aggregates are also not preferentially located at the poles. Finally, we show that the inability of cells with increased viscosity to exclude aggregates from midcell results in enhanced aggregate concentration in between the nucleoids in cells close to dividing. This weakens the asymmetries in aggregate numbers between sister cells of subsequent generations required for rejuvenating cell lineages. We conclude that the process of exclusion of protein aggregates from midcell is not immune to stress conditions affecting the cytoplasm viscosity. The findings contribute to our understanding of E. coli's internal organisation and functioning, and its fragility to stressful conditions.