Linking the Human Gut Microbiome to Inflammatory Cytokine Production Capacity

Gut microbial dysbioses are linked to aberrant immune responses, which are often accompanied by abnormal production of inflammatory cytokines. As part of the Human Functional Genomics Project (HFGP), we investigate how differences in composition and function of gut microbial communities may contribute to inter-individual variation in cytokine responses to microbial stimulations in healthy humans. We observe microbiome-cytokine interaction patterns that are stimulus specific, cytokine specific, and cytokine and stimulus specific. Validation of two predicted host-microbial interactions reveal that TNFα and IFNγ production are associated with specific microbial metabolic pathways: palmitoleic acid metabolism and tryptophan degradation to tryptophol. Besides providing a resource of predicted microbially derived mediators that influence immune phenotypes in response to common microorganisms, these data can help to define principles for understanding disease susceptibility. The three HFGP studies presented in this issue lay the groundwork for further studies aimed at understanding the interplay between microbial, genetic, and environmental factors in the regulation of the immune response in humans. PAPERCLIP.

The effect of host genetics on the gut microbiome

The gut microbiome is affected by multiple factors, including genetics. In this study, we assessed the influence of host genetics on microbial species, pathways and gene ontology categories, on the basis of metagenomic sequencing in 1,514 subjects. In a genome-wide analysis, we identified associations of 9 loci with microbial taxonomies and 33 loci with microbial pathways and gene ontology terms at P < 5 × 10^-8. Additionally, in a targeted analysis of regions involved in complex diseases, innate and adaptive immunity, or food preferences, 32 loci were identified at the suggestive level of P < 5 × 10^-6. Most of our reported associations are new, including genome-wide significance for the C-type lectin molecules CLEC4F-CD207 at 2p13.3 and CLEC4A-FAM90A1 at 12p13. We also identified association of a functional LCT SNP with the Bifidobacterium genus (P = 3.45 × 10^-8) and provide evidence of a gene-diet interaction in the regulation of Bifidobacterium abundance. Our results demonstrate the importance of understanding host-microbe interactions to gain better insight into human health.

Lytic bacteriophage PM16 specific for Proteus mirabilis: a novel member of the genus Phikmvvirus

Lytic Proteus phage PM16, isolated from human faeces, is a novel virus that is specific for Proteus mirabilis cells. Bacteriophage PM16 is characterized by high stability, a short latency period, large burst size and the occurrence of low phage resistance. Phage PM16 was classified as a member of the genus Phikmvvirus on the basis of genome organization, gene synteny, and protein sequences similarities. Within the genus Phikmvvirus, phage PM16 is grouped with Vibrio phage VP93, Pantoea phage LIMElight, Acinetobacter phage Petty, Enterobacter phage phiKDA1, and KP34-like bacteriophages. An investigation of the phage-cell interaction demonstrated that phage PM16 attached to the cell surface, not to the bacterial flagella. The study of P. mirabilis mutant cells obtained during the phage-resistant bacterial cell assay that were resistant to phage PM16 re-infection revealed a non-swarming phenotype, changes in membrane characteristics, and the absence of flagella. Presumably, the resistance of non-swarming P. mirabilis cells to phage PM16 re-infection is determined by changes in membrane macromolecular composition and is associated with the absence of flagella and a non-swarming phenotype.

Population-level analysis of gut microbiome variation

Fecal microbiome variation in the average, healthy population has remained under-investigated. Here, we analyzed two independent, extensively phenotyped cohorts: the Belgian Flemish Gut Flora Project (FGFP; discovery cohort; N = 1106) and the Dutch LifeLines-DEEP study (LLDeep; replication; N = 1135). Integration with global data sets (N combined = 3948) revealed a 14-genera core microbiota, but the 664 identified genera still underexplore total gut diversity. Sixty-nine clinical and questionnaire-based covariates were found associated to microbiota compositional variation with a 92% replication rate. Stool consistency showed the largest effect size, whereas medication explained largest total variance and interacted with other covariate-microbiota associations. Early-life events such as birth mode were not reflected in adult microbiota composition. Finally, we found that proposed disease marker genera associated to host covariates, urging inclusion of the latter in study design.

Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity

Deep sequencing of the gut microbiomes of 1135 participants from a Dutch population-based cohort shows relations between the microbiome and 126 exogenous and intrinsic host factors, including 31 intrinsic factors, 12 diseases, 19 drug groups, 4 smoking categories, and 60 dietary factors. These factors collectively explain 18.7% of the variation seen in the interindividual distance of microbial composition. We could associate 110 factors to 125 species and observed that fecal chromogranin A (CgA), a protein secreted by enteroendocrine cells, was exclusively associated with 61 microbial species whose abundance collectively accounted for 53% of microbial composition. Low CgA concentrations were seen in individuals with a more diverse microbiome. These results are an important step toward a better understanding of environment-diet-microbe-host interactions.

Comparative Metagenomic Profiling of Symbiotic Bacterial Communities Associated with Ixodes persulcatus, Ixodes pavlovskyi and Dermacentor reticulatus Ticks

Ixodes persulcatus, Ixodes pavlovskyi, and Dermacentor reticulatus ticks inhabiting Western Siberia are responsible for the transmission of a number of etiological agents that cause human and animal tick-borne diseases. Because these ticks are abundant in the suburbs of large cities, agricultural areas, and popular tourist sites and frequently attack people and livestock, data regarding the microbiomes of these organisms are required. Using metagenomic 16S profiling, we evaluate bacterial communities associated with I. persulcatus, I. pavlovskyi, and D. reticulatus ticks collected from the Novosibirsk region of Russia. A total of 1214 ticks were used for this study. DNA extracted from the ticks was pooled according to tick species and sex. Sequencing of the V3-V5 domains of 16S rRNA genes was performed using the Illumina Miseq platform. The following bacterial genera were prevalent in the examined communities: Acinetobacter (all three tick species), Rickettsia (I. persulcatus and D. reticulatus) and Francisella (D. reticulatus). B. burgdorferi sensu lato and B. miyamotoi sequences were detected in I. persulcatus and I. pavlovskyi but not in D. reticulatus ticks. The pooled samples of all tick species studied contained bacteria from the Anaplasmataceae family, although their occurrence was low. DNA from A. phagocytophilum and Candidatus Neoehrlichia mikurensis was first observed in I. pavlovskyi ticks. Significant inter-species differences in the number of bacterial taxa as well as intra-species diversity related to tick sex were observed. The bacterial communities associated with the I. pavlovskyi ticks displayed a higher biodiversity compared with those of the I. persulcatus and D. reticulatus ticks. Bacterial community structure was also diverse across the studied tick species, as shown by permutational analysis of variance using the Bray-Curtis dissimilarity metric (p = 0.002). Between-sex variation was confirmed by PERMANOVA testing in I. persulcatus (p = 0.042) and I. pavlovskyi (p = 0.042) ticks. Our study indicated that 16S metagenomic profiling could be used for rapid assessment of the occurrence of medically important bacteria in tick populations inhabiting different natural biotopes and therefore the epidemic danger of studied foci.

MicroRNA drop in the bloodstream and microRNA boost in the tumour caused by treatment with ribonuclease A leads to an attenuation of tumour malignancy

Novel data showing an important role of microRNAs in mediating tumour progression opened a new field of possible molecular targets for cytotoxic ribonucleases. Recently, antitumour and antimetastatic activities of pancreatic ribonuclease A were demonstrated and here genome-wide profiles of microRNAs in the tumour and blood of mice bearing Lewis lung carcinoma after treatment with RNase A were analysed by high-throughput Sequencing by Oligonucleotide Ligation and Detection (SOLiD™) sequencing technology. Sequencing data showed that RNase A therapy resulted in the boost of 116 microRNAs in tumour tissue and a significant drop of 137 microRNAs in the bloodstream that were confirmed by qPCR. The microRNA boost in the tumour was accompanied by the overexpression of microRNA processing genes: RNASEN (Drosha), xpo5, dicer1, and eif2c2 (Ago2). Ribonuclease activity of RNase A was shown to be crucial for the activation of both microRNA synthesis and expression of the microRNA processing genes. In the tumour tissue, RNase A caused the upregulation of both oncomirs and tumour-suppressor microRNAs, including microRNAs of the let-7 family, known to negatively regulate tumour progression. Our results suggest that the alteration of microRNA signature caused by RNase A treatment leads to the attenuation of tumour malignancy.