Human milk variation is shaped by maternal genetics and impacts the infant gut microbiome.

Human milk is a complex mix of nutritional and bioactive components that provide complete nutrition for the infant. However, we lack a systematic knowledge of the factors shaping milk composition and how milk variation influences infant health. Here, we used multi-omic profiling to characterize interactions between maternal genetics, milk gene expression, milk composition, and the infant fecal microbiome in 242 exclusively breastfeeding mother-infant pairs. We identified 487 genetic loci associated with milk gene expression unique to the lactating mammary gland, including loci that impacted breast cancer risk and human milk oligosaccharide concentration. Integrative analyses uncovered connections between milk gene expression and infant gut microbiome, including an association between the expression of inflammation-related genes with IL-6 concentration in milk and the abundance of Bifidobacteria in the infant gut. Our results show how an improved understanding of the genetics and genomics of human milk connects lactation biology with maternal and infant health.

Identification of shared and disease-specific host gene-microbiome associations across human diseases using multi-omic integration.

While gut microbiome and host gene regulation independently contribute to gastrointestinal disorders, it is unclear how the two may interact to influence host pathophysiology. Here we developed a machine learning-based framework to jointly analyse paired host transcriptomic (n = 208) and gut microbiome (n = 208) profiles from colonic mucosal samples of patients with colorectal cancer, inflammatory bowel disease and irritable bowel syndrome. We identified associations between gut microbes and host genes that depict shared as well as disease-specific patterns. We found that a common set of host genes and pathways implicated in gastrointestinal inflammation, gut barrier protection and energy metabolism are associated with disease-specific gut microbes. Additionally, we also found that mucosal gut microbes that have been implicated in all three diseases, such as Streptococcus, are associated with different host pathways in each disease, suggesting that similar microbes can affect host pathophysiology in a disease-specific manner through regulation of different host genes. Our framework can be applied to other diseases for the identification of host gene-microbiome associations that may influence disease outcomes.

Interspecies variation in hominid gut microbiota controls host gene regulation.

The gut microbiome exhibits extreme compositional variation between hominid hosts. However, it is unclear how this variation impacts host physiology across species and whether this effect can be mediated through microbial regulation of host gene expression in interacting epithelial cells. Here, we characterize the transcriptional response of human colonic epithelial cells in vitro to live microbial communities extracted from humans, chimpanzees, gorillas, and orangutans. We find that most host genes exhibit a conserved response, whereby they respond similarly to the four hominid microbiomes. However, hundreds of host genes exhibit a divergent response, whereby they respond only to microbiomes from specific host species. Such genes are associated with intestinal diseases in humans, including inflammatory bowel disease and Crohn's disease. Last, we find that inflammation-associated microbial species regulate the expression of host genes previously associated with inflammatory bowel disease, suggesting health-related consequences for species-specific host-microbiome interactions across hominids.

Longitudinal Multi-omics Reveals Subset-Specific Mechanisms Underlying Irritable Bowel Syndrome.

The gut microbiome has been implicated in multiple human chronic gastrointestinal (GI) disorders. Determining its mechanistic role in disease has been difficult due to apparent disconnects between animal and human studies and lack of an integrated multi-omics view of disease-specific physiological changes. We integrated longitudinal multi-omics data from the gut microbiome, metabolome, host epigenome, and transcriptome in the context of irritable bowel syndrome (IBS) host physiology. We identified IBS subtype-specific and symptom-related variation in microbial composition and function. A subset of identified changes in microbial metabolites correspond to host physiological mechanisms that are relevant to IBS. By integrating multiple data layers, we identified purine metabolism as a novel host-microbial metabolic pathway in IBS with translational potential. Our study highlights the importance of longitudinal sampling and integrating complementary multi-omics data to identify functional mechanisms that can serve as therapeutic targets in a comprehensive treatment strategy for chronic GI diseases. VIDEO ABSTRACT.

Interactions between the gut microbiome and host gene regulation in cystic fibrosis.

BACKGROUND: Cystic fibrosis is the most common autosomal recessive genetic disease in Caucasians. It is caused by mutations in the CFTR gene, leading to poor hydration of mucus and impairment of the respiratory, digestive, and reproductive organ functions. Advancements in medical care have led to markedly increased longevity of patients with cystic fibrosis, but new complications have emerged, such as early onset of colorectal cancer. Although the pathogenesis of colorectal cancer in cystic fibrosis remains unclear, altered host-microbe interactions might play a critical role. To investigate this, we characterized changes in the microbiome and host gene expression in the colonic mucosa of cystic fibrosis patients relative to healthy controls, and identified host gene-microbiome interactions in the colon of cystic fibrosis patients. METHODS: We performed RNA-seq on colonic mucosa samples from cystic fibrosis patients and healthy controls to determine differentially expressed host genes. We also performed 16S rRNA sequencing to characterize the colonic mucosal microbiome and identify gut microbes that are differentially abundant between patients and healthy controls. Lastly, we modeled associations between relative abundances of specific bacterial taxa in the gut mucosa and host gene expression. RESULTS: We find that 1543 genes, including CFTR, show differential expression in the colon of cystic fibrosis patients compared to healthy controls. These genes are enriched with functions related to gastrointestinal and colorectal cancer, such as metastasis of colorectal cancer, tumor suppression, p53, and mTOR signaling pathways. In addition, patients with cystic fibrosis show decreased gut microbial diversity, decreased abundance of butyrate producing bacteria, such as Ruminococcaceae and Butyricimonas, and increased abundance of other taxa, such as Actinobacteria and Clostridium. An integrative analysis identified colorectal cancer-related genes, including LCN2 and DUOX2, for which gene expression is correlated with the abundance of colorectal cancer-associated bacteria, such as Ruminococcaceae and Veillonella. CONCLUSIONS: In addition to characterizing host gene expression and mucosal microbiome in cystic fibrosis patients, our study explored the potential role of host-microbe interactions in the etiology of colorectal cancer in cystic fibrosis. Our results provide biomarkers that may potentially serve as targets for stratifying risk of colorectal cancer in patients with cystic fibrosis.

Integrating tumor genomics into studies of the microbiome in colorectal cancer.

Although the gut microbiome has been linked to colorectal cancer (CRC) development, associations of microbial taxa with CRC status are often inconsistent across studies. We have recently shown that tumor genomics, a factor that is rarely incorporated in analyses of the CRC microbiome, has a strong effect on the composition of the microbiota. Here, we discuss these results in the wider context of studies characterizing interaction between host genetics and the microbiome, and describe the implications of our findings for understanding the role of the microbiome in CRC.

Distinct microbes, metabolites, and ecologies define the microbiome in deficient and proficient mismatch repair colorectal cancers.

BACKGROUND: Links between colorectal cancer (CRC) and the gut microbiome have been established, but the specific microbial species and their role in carcinogenesis remain an active area of inquiry. Our understanding would be enhanced by better accounting for tumor subtype, microbial community interactions, metabolism, and ecology. METHODS: We collected paired colon tumor and normal-adjacent tissue and mucosa samples from 83 individuals who underwent partial or total colectomies for CRC. Mismatch repair (MMR) status was determined in each tumor sample and classified as either deficient MMR (dMMR) or proficient MMR (pMMR) tumor subtypes. Samples underwent 16S rRNA gene sequencing and a subset of samples from 50 individuals were submitted for targeted metabolomic analysis to quantify amino acids and short-chain fatty acids. A PERMANOVA was used to identify the biological variables that explained variance within the microbial communities. dMMR and pMMR microbial communities were then analyzed separately using a generalized linear mixed effects model that accounted for MMR status, sample location, intra-subject variability, and read depth. Genome-scale metabolic models were then used to generate microbial interaction networks for dMMR and pMMR microbial communities. We assessed global network properties as well as the metabolic influence of each microbe within the dMMR and pMMR networks. RESULTS: We demonstrate distinct roles for microbes in dMMR and pMMR CRC. Bacteroides fragilis and sulfidogenic Fusobacterium nucleatum were significantly enriched in dMMR CRC, but not pMMR CRC. These findings were further supported by metabolic modeling and metabolomics indicating suppression of B. fragilis in pMMR CRC and increased production of amino acid proxies for hydrogen sulfide in dMMR CRC. CONCLUSIONS: Integrating tumor biology and microbial ecology highlighted distinct microbial, metabolic, and ecological properties unique to dMMR and pMMR CRC. This approach could critically improve our ability to define, predict, prevent, and treat colorectal cancers.

Transposon mutagenesis screen in mice identifies TM9SF2 as a novel colorectal cancer oncogene.

New therapeutic targets for advanced colorectal cancer (CRC) are critically needed. Our laboratory recently performed an insertional mutagenesis screen in mice to identify novel CRC driver genes and, thus, potential drug targets. Here, we define Transmembrane 9 Superfamily 2 (TM9SF2) as a novel CRC oncogene. TM9SF2 is an understudied protein, belonging to a well conserved protein family characterized by their nine putative transmembrane domains. Based on our transposon screen we found that TM9SF2 is a candidate progression driver in digestive tract tumors. Analysis of The Cancer Genome Atlas (TCGA) data revealed that approximately 35% of CRC patients have elevated levels of TM9SF2 mRNA, data we validated using an independent set of CRC samples. RNAi silencing of TM9SF2 reduced CRC cell growth in an anchorage-independent manner, a hallmark of cancer. Furthermore, CRISPR/Cas9 knockout of TM9SF2 substantially diminished CRC tumor fitness in vitro and in vivo. Transcriptome analysis of TM9SF2 knockout cells revealed a potential role for TM9SF2 in cell cycle progression, oxidative phosphorylation, and ceramide signaling. Lastly, we report that increased TM9SF2 expression correlates with disease stage and low TM9SF2 expression correlate with a more favorable relapse-free survival. Taken together, this study provides evidence that TM9SF2 is a novel CRC oncogene.

Colorectal cancer mutational profiles correlate with defined microbial communities in the tumor microenvironment.

Variation in the gut microbiome has been linked to colorectal cancer (CRC), as well as to host genetic variation. However, we do not know whether, in addition to baseline host genetics, somatic mutational profiles in CRC tumors interact with the surrounding tumor microbiome, and if so, whether these changes can be used to understand microbe-host interactions with potential functional biological relevance. Here, we characterized the association between CRC microbial communities and tumor mutations using microbiome profiling and whole-exome sequencing in 44 pairs of tumors and matched normal tissues. We found statistically significant associations between loss-of-function mutations in tumor genes and shifts in the abundances of specific sets of bacterial taxa, suggestive of potential functional interaction. This correlation allows us to statistically predict interactions between loss-of-function tumor mutations in cancer-related genes and pathways, including MAPK and Wnt signaling, solely based on the composition of the microbiome. In conclusion, our study shows that CRC microbiomes are correlated with tumor mutational profiles, pointing towards possible mechanisms of molecular interaction.

Interaction between Host MicroRNAs and the Gut Microbiota in Colorectal Cancer.

Although variation in gut microbiome composition has been linked with colorectal cancer (CRC), the factors that mediate the interactions between CRC tumors and the microbiome are poorly understood. MicroRNAs (miRNAs) are known to regulate CRC progression and are associated with patient survival outcomes. In addition, recent studies suggested that host miRNAs can also regulate bacterial growth and influence the composition of the gut microbiome. Here, we investigated the association between miRNA expression and microbiome composition in human CRC tumor and normal tissues. We identified 76 miRNAs as differentially expressed (DE) in tissue from CRC tumors and normal tissue, including the known oncogenic miRNAs miR-182, miR-503, and mir-17~92 cluster. These DE miRNAs were correlated with the relative abundances of several bacterial taxa, including Firmicutes, Bacteroidetes, and Proteobacteria. Bacteria correlated with DE miRNAs were enriched with distinct predicted metabolic categories. Additionally, we found that miRNAs that correlated with CRC-associated bacteria are predicted to regulate targets that are relevant for host-microbiome interactions and highlight a possible role for miRNA-driven glycan production in the recruitment of pathogenic microbial taxa. Our work characterized a global relationship between microbial community composition and miRNA expression in human CRC tissues. IMPORTANCE Recent studies have found an association between colorectal cancer (CRC) and the gut microbiota. One potential mechanism by which the microbiota can influence host physiology is through affecting gene expression in host cells. MicroRNAs (miRNAs) are small noncoding RNA molecules that can regulate gene expression and have important roles in cancer development. Here, we investigated the link between the gut microbiota and the expression of miRNA in CRC. We found that dozens of miRNAs are differentially regulated in CRC tumors and adjacent normal colon and that these miRNAs are correlated with the abundance of microbes in the tumor microenvironment. Moreover, we found that microbes that have been previously associated with CRC are correlated with miRNAs that regulate genes related to interactions with microbes. Notably, these miRNAs likely regulate glycan production, which is important for the recruitment of pathogenic microbial taxa to the tumor. This work provides a first systems-level map of the association between microbes and host miRNAs in the context of CRC and provides targets for further experimental validation and potential interventions.

Genetic and transcriptional analysis of human host response to healthy gut microbiota.

Many studies have demonstrated the importance of the gut microbiota in healthy and disease states. However, establishing the causality of host-microbiota interactions in humans is still challenging. Here, we describe a novel experimental system to define the transcriptional response induced by the microbiota in human cells and to shed light on the molecular mechanisms underlying host-gut microbiota interactions. In primary human colonic epithelial cells, we identified over 6,000 genes that change expression at various time points following co-culturing with the gut microbiota of a healthy individual. Among the differentially expressed genes we found a 1.8-fold enrichment of genes associated with diseases that have been previously linked to the microbiome, such as obesity and colorectal cancer. In addition, our experimental system allowed us to identify 87 host SNPs that show allele-specific expression in 69 genes. Furthermore, for 12 SNPs in 12 different genes, allele-specific expression is conditional on the exposure to the microbiota. Of these 12 genes, eight have been associated with diseases linked to the gut microbiota, specifically colorectal cancer, obesity and type 2 diabetes. Our study demonstrates a scalable approach to study host-gut microbiota interactions and can be used to identify putative mechanisms for the interplay between host genetics and microbiota in health and disease.

Virulence genes are a signature of the microbiome in the colorectal tumor microenvironment.

BACKGROUND: The human gut microbiome is associated with the development of colon cancer, and recent studies have found changes in the microbiome in cancer patients compared to healthy controls. Studying the microbial communities in the tumor microenvironment may shed light on the role of host-bacteria interactions in colorectal cancer. Here, we highlight the major shifts in the colorectal tumor microbiome relative to that of matched normal colon tissue from the same individual, allowing us to survey the microbial communities in the tumor microenvironment and providing intrinsic control for environmental and host genetic effects on the microbiome. METHODS: We sequenced the microbiome in 44 primary tumor and 44 patient-matched normal colon tissue samples to determine differentially abundant microbial taxa These data were also used to functionally characterize the microbiome of the cancer and normal sample pairs and identify functional pathways enriched in the tumor-associated microbiota. RESULTS: We find that tumors harbor distinct microbial communities compared to nearby healthy tissue. Our results show increased microbial diversity in the tumor microenvironment, with changes in the abundances of commensal and pathogenic bacterial taxa, including Fusobacterium and Providencia. While Fusobacterium has previously been implicated in colorectal cancer, Providencia is a novel tumor-associated agent which has not been identified in previous studies. Additionally, we identified a clear, significant enrichment of predicted virulence-associated genes in the colorectal cancer microenvironment, likely dependent upon the genomes of Fusobacterium and Providencia. CONCLUSIONS: This work identifies bacterial taxa significantly correlated with colorectal cancer, including a novel finding of an elevated abundance of Providencia in the tumor microenvironment. We also describe the predicted metabolic pathways and enzymes differentially present in the tumor-associated microbiome, and show an enrichment of virulence-associated bacterial genes in the tumor microenvironment. This predicted virulence enrichment supports the hypothesis that the microbiome plays an active role in colorectal cancer development and/or progression. Our results provide a starting point for future prognostic and therapeutic research with the potential to improve patient outcomes.

Functional comparison of innate immune signaling pathways in primates.

Humans respond differently than other primates to a large number of infections. Differences in susceptibility to infectious agents between humans and other primates are probably due to inter-species differences in immune response to infection. Consistent with that notion, genes involved in immunity-related processes are strongly enriched among recent targets of positive selection in primates, suggesting that immune responses evolve rapidly, yet providing only indirect evidence for possible inter-species functional differences. To directly compare immune responses among primates, we stimulated primary monocytes from humans, chimpanzees, and rhesus macaques with lipopolysaccharide (LPS) and studied the ensuing time-course regulatory responses. We find that, while the universal Toll-like receptor response is mostly conserved across primates, the regulatory response associated with viral infections is often lineage-specific, probably reflecting rapid host-virus mutual adaptation cycles. Additionally, human-specific immune responses are enriched for genes involved in apoptosis, as well as for genes associated with cancer and with susceptibility to infectious diseases or immune-related disorders. Finally, we find that chimpanzee-specific immune signaling pathways are enriched for HIV-interacting genes. Put together, our observations lend strong support to the notion that lineage-specific immune responses may help explain known inter-species differences in susceptibility to infectious diseases.

A combination of genomic approaches reveals the role of FOXO1a in regulating an oxidative stress response pathway.

BACKGROUND: While many of the phenotypic differences between human and chimpanzee may result from changes in gene regulation, only a handful of functionally important regulatory differences are currently known. As a first step towards identifying transcriptional pathways that have been remodeled in the human lineage, we focused on a transcription factor, FOXO1a, which we had previously found to be up-regulated in the human liver compared to that of three other primate species. We concentrated on this gene because of its known role in the regulation of metabolism and in longevity. METHODOLOGY: Using a combination of expression profiling following siRNA knockdown and chromatin immunoprecipitation in a human liver cell line, we identified eight novel direct transcriptional targets of FOXO1a. This set includes the gene for thioredoxin-interacting protein (TXNIP), the expression of which is directly repressed by FOXO1a. The thioredoxin-interacting protein is known to inhibit the reducing activity of thioredoxin (TRX), thereby hindering the cellular response to oxidative stress and affecting life span. CONCLUSIONS: Our results provide an explanation for the repeated observations that differences in the regulation of FOXO transcription factors affect longevity. Moreover, we found that TXNIP is down-regulated in human compared to chimpanzee, consistent with the up-regulation of its direct repressor FOXO1a in humans, and with differences in longevity between the two species.