Curtis Huttenhower

Metagenomic biomarker discovery and explanation

This study describes and validates a new method for metagenomic biomarker discovery by way of class comparison, tests of biological consistency and effect size estimation. This addresses the challenge of finding organisms, genes, or pathways that consistently explain the differences between two or more microbial communities, which is a central problem to the study of metagenomics. We extensively validate our method on several microbiomes and a convenient online interface for the method is provided at http://huttenhower.sph.harvard.edu/lefse/.

Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment

BackgroundThe inflammatory bowel diseases (IBD) Crohns 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 Crohns 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.

Genomic analysis identifies association of Fusobacterium with colorectal carcinoma

The tumor microenvironment of colorectal carcinoma is a complex community of genomically altered cancer cells, nonneoplastic cells, and a diverse collection of microorganisms. Each of these components may contribute to carcinogenesis; however, the role of the microbiota is the least well understood. We have characterized the composition of the microbiota in colorectal carcinoma using whole genome sequences from nine tumor/normal pairs. Fusobacterium sequences were enriched in carcinomas, confirmed by quantitative PCR and 16S rDNA sequence analysis of 95 carcinoma/normal DNA pairs, while the Bacteroidetes and Firmicutes phyla were depleted in tumors. Fusobacteria were also visualized within colorectal tumors using FISH. These findings reveal alterations in the colorectal cancer microbiota; however, the precise role of Fusobacteria in colorectal carcinoma pathogenesis requires further investigation.

Metagenomic microbial community profiling using unique clade-specific marker genes

Metagenomic shotgun sequencing data can identify microbes populating a microbial community and their proportions, but existing taxonomic profiling methods are inefficient for increasingly large data sets. We present an approach that uses clade-specific marker genes to unambiguously assign reads to microbial clades more accurately and >50× faster than current approaches. We validated our metagenomic phylogenetic analysis tool, MetaPhlAn, on terabases of short reads and provide the largest metagenomic profiling to date of the human gut. It can be accessed at http://huttenhower.sph.harvard.edu/metaphlan/.

Microbial Co-occurrence Relationships in the Human Microbiome

The healthy microbiota show remarkable variability within and among individuals. In addition to external exposures, ecological relationships (both oppositional and symbiotic) between microbial inhabitants are important contributors to this variation. It is thus of interest to assess what relationships might exist among microbes and determine their underlying reasons. The initial Human Microbiome Project (HMP) cohort, comprising 239 individuals and 18 different microbial habitats, provides an unprecedented resource to detect, catalog, and analyze such relationships. Here, we applied an ensemble method based on multiple similarity measures in combination with generalized boosted linear models (GBLMs) to taxonomic marker (16S rRNA gene) profiles of this cohort, resulting in a global network of 3,005 significant co-occurrence and co-exclusion relationships between 197 clades occurring throughout the human microbiome. This network revealed strong niche specialization, with most microbial associations occurring within body sites and a number of accompanying inter-body site relationships. Microbial communities within the oropharynx grouped into three distinct habitats, which themselves showed no direct influence on the composition of the gut microbiota. Conversely, niches such as the vagina demonstrated little to no decomposition into region-specific interactions. Diverse mechanisms underlay individual interactions, with some such as the co-exclusion of Porphyromonaceae family members and Streptococcus in the subgingival plaque supported by known biochemical dependencies. These differences varied among broad phylogenetic groups as well, with the Bacilli and Fusobacteria, for example, both enriched for exclusion of taxa from other clades. Comparing phylogenetic versus functional similarities among bacteria, we show that dominant commensal taxa (such as Prevotellaceae and Bacteroides in the gut) often compete, while potential pathogens (e.g. Treponema and Prevotella in the dental plaque) are more likely to co-occur in complementary niches. This approach thus serves to open new opportunities for future targeted mechanistic studies of the microbial ecology of the human microbiome.

Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis

Rheumatoid arthritis (RA) is a prevalent systemic autoimmune disease, caused by a combination of genetic and environmental factors. Animal models suggest a role for intestinal bacteria in supporting the systemic immune response required for joint inflammation. Here we performed 16S sequencing on 114 stool samples from rheumatoid arthritis patients and controls, and shotgun sequencing on a subset of 44 such samples. We identified the presence of Prevotella copri as strongly correlated with disease in new-onset untreated rheumatoid arthritis (NORA) patients. Increases in Prevotella abundance correlated with a reduction in Bacteroides and a loss of reportedly beneficial microbes in NORA subjects. We also identified unique Prevotella genes that correlated with disease. Further, colonization of mice revealed the ability of P. copri to dominate the intestinal microbiota and resulted in an increased sensitivity to chemically induced colitis. This work identifies a potential role for P. copri in the pathogenesis of RA. DOI: http://dx.doi.org/10.7554/eLife.01202.001

Analysis of phosphorylation sites on proteins from Saccharomyces cerevisiae by electron transfer dissociation (ETD) mass spectrometry

We present a strategy for the analysis of the yeast phosphoproteome that uses endo-Lys C as the proteolytic enzyme, immobilized metal affinity chromatography for phosphopeptide enrichment, a 90-min nanoflow-HPLC/electrospray-ionization MS/MS experiment for phosphopeptide fractionation and detection, gas phase ion/ion chemistry, electron transfer dissociation for peptide fragmentation, and the Open Mass Spectrometry Search Algorithm for phosphoprotein identification and assignment of phosphorylation sites. From a 30-μg (≈600 pmol) sample of total yeast protein, we identify 1,252 phosphorylation sites on 629 proteins. Identified phosphoproteins have expression levels that range from <50 to 1,200,000 copies per cell and are encoded by genes involved in a wide variety of cellular processes. We identify a consensus site that likely represents a motif for one or more uncharacterized kinases and show that yeast kinases, themselves, contain a disproportionately large number of phosphorylation sites. Detection of a pHis containing peptide from the yeast protein, Cdc10, suggests an unexpected role for histidine phosphorylation in septin biology. From diverse functional genomics data, we show that phosphoproteins have a higher number of interactions than an average protein and interact with each other more than with a random protein. They are also likely to be conserved across large evolutionary distances.

Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples

BackgroundTo understand the relationship between our bacterial microbiome and health, it is essential to define the microbiome in the absence of disease. The digestive tract includes diverse habitats and hosts the human bodys greatest bacterial density. We describe the bacterial community composition of ten digestive tract sites from more than 200 normal adults enrolled in the Human Microbiome Project, and metagenomically determined metabolic potentials of four representative sites.ResultsThe microbiota of these diverse habitats formed four groups based on similar community compositions: buccal mucosa, keratinized gingiva, hard palate; saliva, tongue, tonsils, throat; sub- and supra-gingival plaques; and stool. Phyla initially identified from environmental samples were detected throughout this population, primarily TM7, SR1, and Synergistetes. Genera with pathogenic members were well-represented among this disease-free cohort. Tooth-associated communities were distinct, but not entirely dissimilar, from other oral surfaces. The Porphyromonadaceae, Veillonellaceae and Lachnospiraceae families were common to all sites, but the distributions of their genera varied significantly. Most metabolic processes were distributed widely throughout the digestive tract microbiota, with variations in metagenomic abundance between body habitats. These included shifts in sugar transporter types between the supragingival plaque, other oral surfaces, and stool; hydrogen and hydrogen sulfide production were also differentially distributed.ConclusionsThe microbiomes of ten digestive tract sites separated into four types based on composition. A core set of metabolic pathways was present across these diverse digestive tract habitats. These data provide a critical baseline for future studies investigating local and systemic diseases affecting human health.

The Dynamics of the Human Infant Gut Microbiome in Development and in Progression toward Type 1 Diabetes

Colonization of the fetal and infant gut microbiome results in dynamic changes in diversity, which can impact disease susceptibility. To examine the relationship between human gut microbiome dynamics throughout infancy and type 1 diabetes (T1D), we examined a cohort of 33 infants genetically predisposed to T1D. Modeling trajectories of microbial abundances through infancy revealed a subset of microbial relationships shared across most subjects. Although strain composition of a given species was highly variable between individuals, it was stable within individuals throughout infancy. Metabolic composition and metabolic pathway abundance remained constant across time. A marked drop in alpha-diversity was observed in T1D progressors in the time window between seroconversion and T1D diagnosis, accompanied by spikes in inflammation-favoring organisms, gene functions, and serum and stool metabolites. This work identifies trends in the development of the human infant gut microbiome along with specific alterations that precede T1D onset and distinguish T1D progressors from nonprogressors.

Assessment of colorectal cancer molecular features along bowel subsites challenges the conception of distinct dichotomy of proximal versus distal colorectum

Objective Colorectal cancer is typically classified into proximal colon, distal colon and rectal cancer. Tumour genetic and epigenetic features differ by tumour location. Considering a possible role of bowel contents (including microbiome) in carcinogenesis, this study hypothesised that tumour molecular features might gradually change along bowel subsites, rather than change abruptly at splenic flexure. Design Utilising 1443 colorectal cancers in two US nationwide prospective cohort studies, the frequencies of molecular features (CpG island methylator phenotype (CIMP), microsatellite instability (MSI), LINE-1 methylation and BRAF, KRAS and PIK3CA mutations) were examined along bowel subsites (rectum, rectosigmoid junction, sigmoid, descending colon, splenic flexure, transverse colon, hepatic flexure, ascending colon and caecum). The linearity and non-linearity of molecular relations along subsites were statistically tested by multivariate logistic or linear regression analysis. Results The frequencies of CIMP-high, MSI-high and BRAF mutations gradually increased from the rectum (<2.3%) to ascending colon (36–40%), followed by falls in the caecum (12–22%). By linearity tests, these molecular relations were significantly linear from rectum to ascending colon (p<0.0001), and there was little evidence of non-linearity (p>0.09). Caecal cancers exhibited the highest frequency of KRAS mutations (52% vs 27–35% in other sites; p<0.0001). Conclusions The frequencies of CIMP-high, MSI-high and BRAF mutations in cancer increased gradually along colorectum subsites from the rectum to ascending colon. These novel data challenge the common conception of discrete molecular features of proximal versus distal colorectal cancers, and have a substantial impact on clinical, translational and epidemiology research, which has typically been performed with the dichotomous classification of proximal versus distal tumours.