BackgroundThe inflammatory bowel diseases (IBD) Crohn's disease and ulcerative colitis result from alterations in intestinal microbes and the immune system. However, the precise dysfunctions of microbial metabolism in the gastrointestinal microbiome during IBD remain unclear. We analyzed the microbiota of intestinal biopsies and stool samples from 231 IBD and healthy subjects by 16S gene pyrosequencing and followed up a subset using shotgun metagenomics. Gene and pathway composition were assessed, based on 16S data from phylogenetically-related reference genomes, and associated using sparse multivariate linear modeling with medications, environmental factors, and IBD status.ResultsFirmicutes and Enterobacteriaceae abundances were associated with disease status as expected, but also with treatment and subject characteristics. Microbial function, though, was more consistently perturbed than composition, with 12% of analyzed pathways changed compared with 2% of genera. We identified major shifts in oxidative stress pathways, as well as decreased carbohydrate metabolism and amino acid biosynthesis in favor of nutrient transport and uptake. The microbiome of ileal Crohn's disease was notable for increases in virulence and secretion pathways.ConclusionsThis inferred functional metagenomic information provides the first insights into community-wide microbial processes and pathways that underpin IBD pathogenesis.

To the Editor: MetaPhlAn (metagenomic phylogenetic analysis) 1 is a method for characterizing the taxonomic profiles of whole-metagenome shotgun (WMS) samples that has been used successfully in large-scale microbial community studies 2,3 . This work complements the original species-level profiling method with a system for eukaryotic and viral quantitation, strain-level identification and strain tracking. These and other extensions make the MetaPhlAn2 computational package (http://segatalab. cibio.unitn.it/tools/metaphlan2/ and Supplementary Software) an efficient tool for mining WMS samples.Our method infers the presence and read coverage of cladespecific markers to unequivocally detect the taxonomic clades present in a microbiome sample and estimate their relative abundance 1 . MetaPhlAn2 includes an expanded set of ~1 million markers (184 ± 45 for each bacterial species) from >7,500 species (Supplementary Tables 1-3), based on the approximately tenfold increase in the number of sequenced genomes in the past 2 years. Subspecies markers enable strain-level analyses, and quasi-markers improve accuracy and allow the detection of viruses and eukaryotic microbes (a full list of additions is provided in Supplementary Notes 1-3 and Supplementary Fig. 1).We validated MetaPhlAn2 using 24 synthetic metagenomes comprising 656 million reads and 1,295 species (Supplementary Note 4 and Supplementary Table 4). MetaPhlAn2 proved more accurate (average correlation: 0.95 ± 0.05) than mOTU 4 and Kraken 5 (0.80 ± 0.21 and 0.75 ± 0.22, respectively) ( Fig. 1a, Supplementary Figs. 2-9 and Supplementary Tables 5-11),with fewer false positives (an average of 10, compared with 22 and 23 for mOTU and Kraken, respectively) and false negatives (an average of 12, compared with 27 for the other two methods), even when including genomes that were absent from the reference database (Supplementary Note 4). With the adoption of the BowTie2 fast mapper and support for parallelism, MetaPhlAn2 is more than ten times faster than MetaPhlAn, and its speed is comparable to that of other tested approaches ( Supplementary Fig. 10).We applied MetaPhlAn2 to four elbow-skin samples that we sequenced from three subjects (Fig. 1b, Supplementary Note 5 and Supplementary Table 12). Our data showed that Propionibacterium acnes and Staphylococcus epidermidis dominated these sites, in agreement with expected genus-level results 6 , while providing species-level resolution. Together with these core species, we found Malassezia globosa in 93.65% of samples and confirmed it by coverage analysis (Supplementary Fig. 11). Although M. globosa is a known colonizer of the skin, its metagenomic characterization highlights the ability of MetaPhlAn2 to identify non-prokaryotic species. Phages (e.g., for Propionibacterium) and double-stranded DNA viruses of the Polyomavirus genus were also consistently detected. We subsequently profiled the whole set of 982 samples from other body sites from the Human Microbiome Project (HMP), including 219 samples sequenced after the initi...

This is a PDF file of a peer-reviewed paper that has been accepted for publication. Although unedited, the content has been subjected to preliminary formatting. Nature is providing this early version of the typeset paper as a service to our authors and readers. The text and figures will undergo copyediting and a proof review before the paper is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers apply.

It is challenging to associate features such as human health outcomes, diet, environmental conditions, or other metadata to microbial community measurements, due in part to their quantitative properties. Microbiome multi-omics are typically noisy, sparse (zero-inflated), high-dimensional, extremely non-normal, and often in the form of count or compositional measurements. Here we introduce an optimized combination of novel and established methodology to assess multivariable association of microbial community features with complex metadata in population-scale observational studies. Our approach, MaAsLin 2 (Microbiome Multivariable Associations with Linear Models), uses generalized linear and mixed models to accommodate a wide variety of modern epidemiological studies, including cross-sectional and longitudinal designs, as well as a variety of data types (e.g., counts and relative abundances) with or without covariates and repeated measurements. To construct this method, we conducted a large-scale evaluation of a broad range of scenarios under which straightforward identification of meta-omics associations can be challenging. These simulation studies reveal that MaAsLin 2’s linear model preserves statistical power in the presence of repeated measures and multiple covariates, while accounting for the nuances of meta-omics features and controlling false discovery. We also applied MaAsLin 2 to a microbial multi-omics dataset from the Integrative Human Microbiome (HMP2) project which, in addition to reproducing established results, revealed a unique, integrated landscape of inflammatory bowel diseases (IBD) across multiple time points and omics profiles.

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells/multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Using single cell transcriptome profiling of ~140,000 liver and ~74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the impact of tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, NK and ILC precursors in the yolk sac. We demonstrate a shift in fetal liver haematopoietic composition during gestation away from being erythroid-predominant, accompanied by a parallel change in HSC/MPP differentiation potential, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a valuable reference for harnessing the therapeutic potential of HSC/MPPs.

SUMMARY Mammals have extremely limited regenerative capabilities; however, axolotls are profoundly regenerative and can replace entire limbs. The mechanisms underlying limb regeneration remain poorly understood, partly because the enormous and incompletely sequenced genomes of axolotls have hindered the study of genes facilitating regeneration. We assembled and annotated a de novo transcriptome using RNA-sequencing profiles for a broad spectrum of tissues that is estimated to have near-complete sequence information for 88% of axolotl genes. We devised expression analyses that identified the axolotl orthologs of cirbp and kazald1 as highly expressed and enriched in blastemas. Using morpholino anti-sense oligonucleotides, we find evidence that cirbp plays a cytoprotective role during limb regeneration while manipulation of kazald1 expression disrupts regeneration. Our transcriptome and annotation resources greatly complement previous transcriptomic studies and will be a valuable resource for future research in regenerative biology.

Type 2 innate lymphoid cells (ILC2s) both contribute to mucosal homeostasis and initiate pathologic inflammation in allergic asthma. However, the signals that direct ILC2s to promote homeostasis versus inflammation are unclear. To identify such molecular cues, we profiled mouse lung-resident ILCs using single-cell RNA sequencing at steady state and after in vivo stimulation with the alarmin cytokines IL-25 and IL-33. ILC2s were transcriptionally heterogeneous after activation, with subpopulations distinguished by expression of proliferative, homeostatic and effector genes. The neuropeptide receptor Nmur1 was preferentially expressed by ILC2s at steady state and after IL-25 stimulation. Neuromedin U (NMU), the ligand of NMUR1, activated ILC2s in vitro, and in vivo co-administration of NMU with IL-25 strongly amplified allergic inflammation. Loss of NMU–NMUR1 signalling reduced ILC2 frequency and effector function, and altered transcriptional programs following allergen challenge in vivo. Thus, NMUR1 signalling promotes inflammatory ILC2 responses, highlighting the importance of neuro-immune crosstalk in allergic inflammation at mucosal surfaces.

The increased availability of genomic and metagenomic data poses challenges at multiple analysis levels, including visualization of very large-scale microbial and microbial community data paired with rich metadata. We developed GraPhlAn (Graphical Phylogenetic Analysis), a computational tool that produces high-quality, compact visualizations of microbial genomes and metagenomes. This includes phylogenies spanning up to thousands of taxa, annotated with metadata ranging from microbial community abundances to microbial physiology or host and environmental phenotypes. GraPhlAn has been developed as an open-source command-driven tool in order to be easily integrated into complex, publication-quality bioinformatics pipelines. It can be executed either locally or through an online Galaxy web application. We present several examples including taxonomic and phylogenetic visualization of microbial communities, metabolic functions, and biomarker discovery that illustrate GraPhlAn’s potential for modern microbial and community genomics.