The inflammatory bowel diseases (IBD), which include Crohn’s disease (CD) and ulcerative colitis (UC), are multifactorial, chronic conditions of the gastrointestinal tract. While IBD has been associated with dramatic changes in the gut microbiota, changes in the gut metabolome -- the molecular interface between host and microbiota -- are less-well understood. To address this gap, we performed untargeted LC-MS metabolomic and shotgun metagenomic profiling of cross-sectional stool samples from discovery ( n =155) and validation ( n =65) cohorts of CD, UC, and non-IBD control subjects. Metabolomic and metagenomic profiles were broadly correlated with fecal calprotectin levels (a measure of gut inflammation). Across >8,000 measured metabolite features, we identified chemicals and chemical classes that were differentially abundant (DA) in IBD, including enrichments for sphingolipids and bile acids, and depletions for triacylglycerols and tetrapyrroles. While >50% of DA metabolite features were uncharacterized, many could be assigned putative roles through metabolomic “guilt-by-association” (covariation with known metabolites). DA species and functions from the metagenomic profiles reflected adaptation to oxidative stress in the IBD gut, and were individually consistent with previous findings. Integrating these data, however, we identified 122 robust associations between DA species and well-characterized DA metabolites, indicating possible mechanistic relationships that are perturbed in IBD. Finally, we found that metabolome- and metagenome-based classifiers of IBD status were highly accurate and, like the vast majority of individual trends, generalized well to the independent validation cohort. Our findings thus provide an improved understanding of perturbations of the microbiome-metabolome interface in IBD, including identification of many potential diagnostic and therapeutic targets.

Mathematical reaction-diffusion models have been suggested to describe formation of animal pigmentation patterns and distribution of epidermal appendages. However, the crucial signals and in vivo mechanisms are still elusive. Here we identify WNT and its inhibitor DKK as primary determinants of murine hair follicle spacing, using a combined experimental and computational modeling approach. Transgenic DKK overexpression reduces overall appendage density. Moderate suppression of endogenous WNT signaling forces follicles to form clusters during an otherwise normal morphogenetic program. These results confirm predictions of a WNT/DKK-specific mathematical model and provide in vivo corroboration of the reaction-diffusion mechanism for epidermal appendage formation.

Gene regulatory, signal transduction and metabolic networks are major areas of interest in the newly emerging field of systems biology. In living cells, stochastic dynamics play an important role; however, the kinetic parameters of biochemical reactions necessary for modelling these processes are often not accessible directly through experiments. The problem of estimating stochastic reaction constants from molecule count data measured, with error, at discrete time points is considered. For modelling the system, a hidden Markov process is used, where the hidden states are the true molecule counts, and the transitions between those states correspond to reaction events following collisions of molecules. Two different algorithms are proposed for estimating the unknown model parameters. The first is an approximate maximum likelihood method that gives good estimates of the reaction parameters in systems with few possible reactions in each sampling interval. The second algorithm, treating the data as exact measurements, approximates the number of reactions in each sampling interval by solving a simple linear equation. Maximising the likelihood based on these approximations can provide good results, even in complex reaction systems.

No abstract