Gut-associated functions are favored during microbiome assembly across a major part of C. elegans life.

The microbiome expresses a variety of functions that influence host biology. The range of functions depends on the microbiome's composition, which can change during the host's lifetime due to neutral assembly processes, host-mediated selection, and environmental conditions. To date, the exact dynamics of microbiome assembly, the underlying determinants, and the effects on host-associated functions remain poorly understood. Here, we used the nematode Caenorhabditis elegans and a defined community of fully sequenced, naturally associated bacteria to study microbiome dynamics and functions across a major part of the worm's lifetime of hosts under controlled experimental conditions. Bacterial community composition initially shows strongly declining levels of stochasticity, which increases during later time points, suggesting selective effects in younger animals as opposed to more random processes in older animals. The adult microbiome is enriched in genera Ochrobactrum and Enterobacter compared to the direct substrate and a host-free control environment. Using pathway analysis, metabolic, and ecological modeling, we further find that the lifetime assembly dynamics increase competitive strategies and gut-associated functions in the host-associated microbiome, indicating that the colonizing bacteria benefit the worm. Overall, our study introduces a framework for studying microbiome assembly dynamics based on stochastic, ecological, and metabolic models, yielding new insights into the processes that determine host-associated microbiome composition and function. IMPORTANCE: The microbiome plays a crucial role in host biology. Its functions depend on the microbiome composition that can change during a host's lifetime. To date, the dynamics of microbiome assembly and the resulting functions still need to be better understood. This study introduces a new approach to characterize the functional consequences of microbiome assembly by modeling both the relevance of stochastic processes and metabolic characteristics of microbial community changes. The approach was applied to experimental time-series data obtained for the microbiome of the nematode Caenorhabditis elegans across the major part of its lifetime. Stochastic processes played a minor role, whereas beneficial bacteria as well as gut-associated functions enriched in hosts. This indicates that the host might actively shape the composition of its microbiome. Overall, this study provides a framework for studying microbiome assembly dynamics and yields new insights into C. elegans microbiome functions.

Forensically relevant anatomical brain regions cannot be sub-differentiated by RNA expression analysis.

The contextualization of biological traces generated by severe head injuries can be beneficial for criminal investigations. Here we aimed to identify and validate mRNA candidates for a robust sub-differentiation of forensically and traumatologically relevant brain regions. To this purpose, massively parallel sequencing of whole transcriptomes in sample material taken from four different areas of the cerebral cortex (frontal, temporal, parietal, occipital lobe) was performed, followed by bioinformatical data analysis, classification, and biostatistical candidate selection. Candidates were evaluated by Multiplex-RT-PCR and capillary electrophoresis. Only a weak relative upregulation and solely for candidates expressed in the parietal lobe was observed. Two candidates with upregulation in the cerebellar region (PVALB and CDR2L) were chosen for further investigation; however, PVALB could not reliably and repeatedly be detected in any lobe whereas CDR2L was detectable in all lobes. Consequently, we suggest that differences in mRNA expression between four regions of the cerebral cortex are too small and less pronounced to be useful for and applicable in forensic RNA analysis. We conclude that sub-differentiation of these brain regions via RNA expression analysis is generally not feasible within a forensic scope.

Dupilumab but not cyclosporine treatment shifts the microbiome toward a healthy skin flora in patients with moderate-to-severe atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) patients display an altered skin microbiome which may not only be an indicator but also a driver of inflammation. We aimed to investigate associations among AD patients' skin microbiome, clinical data, and response to systemic therapy in patients of the TREATgermany registry. METHODS: Skin swabs of 157 patients were profiled with 16S rRNA gene amplicon sequencing before and after 3 months of treatment with dupilumab or cyclosporine. For comparison, 16s microbiome data from 258 population-based healthy controls were used. Disease severity was assessed using established instruments such as the Eczema Area and Severity Index (EASI). RESULTS: We confirmed the previously shown correlation of Staphylococcus aureus abundance and bacterial alpha diversity with AD severity as measured by EASI. Therapy with Dupilumab shifted the bacterial community toward the pattern seen in healthy controls. The relative abundance of Staphylococci and in particular S. aureus significantly decreased on both lesional and non-lesional skin, whereas the abundance of Staphylococcus hominis increased. These changes were largely independent from the degree of clinical improvement and were not observed for cyclosporine. CONCLUSIONS: Systemic treatment with dupilumab but not cyclosporine tends to restore a healthy skin microbiome largely independent of the clinical response indicating potential effects of IL-4RA blockade on the microbiome.

Sphingolipids Are Depleted in Alcohol-Related Liver Fibrosis.

BACKGROUND & AIMS: Alcohol disturbs hepatic lipid synthesis and transport, but the role of lipid dysfunction in alcohol-related liver disease (ALD) is unclear. In this biopsy-controlled, prospective, observational study, we characterized the liver and plasma lipidomes in patients with early ALD. METHODS: We performed mass spectrometry-based lipidomics of paired liver and plasma samples from 315 patients with ALD and of plasma from 51 matched healthy controls. We associated lipid levels with histologic fibrosis, inflammation, and steatosis with correction for multiple testing and adjustment for confounders. We further investigated sphingolipid regulation by means of quantitative real-time polymerase chain reaction sequencing of microRNA, prediction of liver-related events, and tested causality with Mendelian randomization. RESULTS: We detected 198 lipids in the liver and 236 lipids in the circulation from 18 lipid classes. Most sphingolipids (sphingomyelins and ceramides) and phosphocholines were co-down-regulated in both liver and plasma, where lower abundance correlated with higher fibrosis stage. Sphingomyelins showed the most pronounced negative correlation to fibrosis, mirrored by negative correlations in both liver and plasma with hepatic inflammation. Reduced sphingomyelins predicted future liver-related events. This seemed to be characteristic of "pure ALD," as sphingomyelin levels were higher in patients with concomitant metabolic syndrome and ALD/nonalcoholic fatty liver disease overlap. Mendelian randomization in FinnGen and UK Biobanks indicated ALD as the cause of low sphingomyelins, and alcohol use disorder did not correlate with genetic susceptibility to low sphingomyelin levels. CONCLUSIONS: Alcohol-related liver fibrosis is characterized by selective and progressive lipid depletion in liver and blood, particularly sphingomyelins, which also associates with progression to liver-related events.

Host genetic factors related to innate immunity, environmental sensing and cellular functions are associated with human skin microbiota.

Despite the increasing knowledge about factors shaping the human microbiome, the host genetic factors that modulate the skin-microbiome interactions are still largely understudied. This contrasts with recent efforts to characterize host genes that influence the gut microbiota. Here, we investigated the effect of genetics on skin microbiota across three different skin microenvironments through meta-analyses of genome-wide association studies (GWAS) of two population-based German cohorts. We identified 23 genome-wide significant loci harboring 30 candidate genes involved in innate immune signaling, environmental sensing, cell differentiation, proliferation and fibroblast activity. However, no locus passed the strict threshold for study-wide significance (P < 6.3 × 10-10 for 80 features included in the analysis). Mendelian randomization (MR) analysis indicated the influence of staphylococci on eczema/dermatitis and suggested modulating effects of the microbiota on other skin diseases. Finally, transcriptional profiles of keratinocytes significantly changed after in vitro co-culturing with Staphylococcus epidermidis, chosen as a representative of skin commensals. Seven candidate genes from the GWAS were found overlapping with differential expression in the co-culturing experiments, warranting further research of the skin commensal and host genetic makeup interaction.

High-fat meals do not affect thrombin formation and fibrin clot lysis in individuals with obesity during intentional weight loss.

Repeated weight loss cycles are associated with increased cardiovascular morbidity. Meal-induced thrombin formation, measured as prothrombin fragment 1+2 (F1+2), is observed in individuals with overweight after weight loss, and postprandial effects can be one of the mechanisms underlying harmful effects during intentional weight loss. We hypothesize that consumption of high-fat meals during intentional weight loss triggers a prothrombotic state by increasing postprandial F1+2 or decreasing fibrin clot lysis in individuals with obesity, and that the response associates with the gut bacteria composition. A cross-over meal study was conducted in patients admitted to bariatric surgery during dietary weight loss (N = 20) and surgical weight loss (N = 16) (weight loss groups). High-fat (67 E%) and low-fat (16 E%) meals were served at 08:15 and 10:00 on 2 study days. Blood samples collected at 08:00 (fasting), 12:00, and 14:00 were analyzed for triglycerides, activated factor VII (FVIIa), F1+2, D-dimer, fibrinogen, tissue factor , and fibrin clot lysis. The proportion of Gram-negative bacteria and bacterial diversity were analyzed in fecal samples obtained less than 24 hours before the meal test. Triglyceride and FVIIa increased after high-fat meals in both weight loss groups, whereas D-dimer (dietary group) and F1+2 decreased and tissue factor and fibrin clot lysis did not change. There was a negative association between the proportion of Gram-negative bacteria and changes in FVIIa in the surgery group. Postprandial FVII activation after high-fat meals is not accompanied by increased F1+2, irrespective of the weight loss intervention, but might be associated with the proportion of Gram-negative gut bacteria.

Combinatorial, additive and dose-dependent drug-microbiome associations.

During the transition from a healthy state to cardiometabolic disease, patients become heavily medicated, which leads to an increasingly aberrant gut microbiome and serum metabolome, and complicates biomarker discovery1-5. Here, through integrated multi-omics analyses of 2,173 European residents from the MetaCardis cohort, we show that the explanatory power of drugs for the variability in both host and gut microbiome features exceeds that of disease. We quantify inferred effects of single medications, their combinations as well as additive effects, and show that the latter shift the metabolome and microbiome towards a healthier state, exemplified in synergistic reduction in serum atherogenic lipoproteins by statins combined with aspirin, or enrichment of intestinal Roseburia by diuretic agents combined with beta-blockers. Several antibiotics exhibit a quantitative relationship between the number of courses prescribed and progression towards a microbiome state that is associated with the severity of cardiometabolic disease. We also report a relationship between cardiometabolic drug dosage, improvement in clinical markers and microbiome composition, supporting direct drug effects. Taken together, our computational framework and resulting resources enable the disentanglement of the effects of drugs and disease on host and microbiome features in multimedicated individuals. Furthermore, the robust signatures identified using our framework provide new hypotheses for drug-host-microbiome interactions in cardiometabolic disease.

The power and potential of BIOMAP to elucidate host-microbiome interplay in skin inflammatory diseases.

The two most common chronic inflammatory skin diseases are atopic dermatitis (AD) and psoriasis. The underpinnings of the remarkable degree of clinical heterogeneity of AD and psoriasis are poorly understood and, as a consequence, disease onset and progression are unpredictable and the optimal type and time point for intervention are as yet unknown. The BIOMAP project is the first IMI (Innovative Medicines Initiative) project dedicated to investigating the causes and mechanisms of AD and psoriasis and to identify potential biomarkers responsible for the variation in disease outcome. The consortium includes 7 large pharmaceutical companies and 25 non-industry partners including academia. Since there is mounting evidence supporting an important role for microbial exposures and our microbiota as factors mediating immune polarization and AD and psoriasis pathogenesis, an entire work package is dedicated to the investigation of skin and gut microbiome linked to AD or psoriasis. The large collaborative BIOMAP project will enable the integration of patient cohorts, data and knowledge in unprecedented proportions. The project has a unique opportunity with a potential to bridge and fill the gaps between current problems and solutions. This review highlights the power and potential of the BIOMAP project in the investigation of microbe-host interplay in AD and psoriasis.

Short-term physical exercise impacts on the human holobiont obtained by a randomised intervention study.

BACKGROUND: Human well-being has been linked to the composition and functional capacity of the intestinal microbiota. As regular exercise is known to improve human health, it is not surprising that exercise was previously described to positively modulate the gut microbiota, too. However, most previous studies mainly focused on either elite athletes or animal models. Thus, we conducted a randomised intervention study that focused on the effects of different types of training (endurance and strength) in previously physically inactive, healthy adults in comparison to controls that did not perform regular exercise. Overall study duration was ten weeks including six weeks of intervention period. In addition to 16S rRNA gene amplicon sequencing of longitudinally sampled faecal material of participants (six time points), detailed body composition measurements and analysis of blood samples (at baseline and after the intervention) were performed to obtain overall physiological changes within the intervention period. Activity tracker devices (wrist-band wearables) provided activity status and sleeping patterns of participants as well as exercise intensity and heart measurements. RESULTS: Different biometric responses between endurance and strength activities were identified, such as a significant increase of lymphocytes and decrease of mean corpuscular haemoglobin concentration (MCHC) only within the strength intervention group. In the endurance group, we observed a significant reduction in hip circumference and an increase in physical working capacity (PWC). Though a large variation of microbiota changes were observed between individuals of the same group, we did not find specific collective alterations in the endurance nor the strength groups, arguing for microbiome variations specific to individuals, and therefore, were not captured in our analysis. CONCLUSIONS: We could show that different types of exercise have distinct but moderate effects on the overall physiology of humans and very distinct microbial changes in the gut. The observed overall changes during the intervention highlight the importance of physical activity on well-being. Future studies should investigate the effect of exercise on a longer timescale, investigate different training intensities and consider high-resolution shotgun metagenomics technology. TRIAL REGISTRATION: DRKS, DRKS00015873 . Registered 12 December 2018; Retrospectively registered.

Altered Gut Microbial Metabolism of Essential Nutrients in Primary Sclerosing Cholangitis.

BACKGROUND & AIMS: To influence host and disease phenotype, compositional microbiome changes, which have been demonstrated in patients with primary sclerosing cholangitis (PSC), must be accompanied by functional changes. We therefore aimed to characterize the genetic potential of the gut microbiome in patients with PSC compared with healthy controls (HCs) and patients with inflammatory bowel disease (IBD). METHODS: Fecal DNA from 2 cohorts (1 Norwegian and 1 German), in total comprising 136 patients with PSC (58% with IBD), 158 HCs, and 93 patients with IBD without PSC, were subjected to metagenomic shotgun sequencing, generating 17 billion paired-end sequences, which were processed using HUMAnN2 and MetaPhlAn2, and analyzed using generalized linear models and random effects meta-analyses. RESULTS: Patients with PSC had fewer microbial genes compared with HCs (P < .0001). Compared with HCs, patients with PSC showed enrichment and increased prevalence of Clostridium species and a depletion of, for example, Eubacterium spp and Ruminococcus obeum. Patients with PSC showed marked differences in the abundance of genes related to vitamin B6 synthesis and branched-chain amino acid synthesis (Qfdr < .05). Targeted metabolomics of plasma from an independent set of patients with PSC and controls found reduced concentrations of vitamin B6 and branched-chain amino acids in PSC (P < .0001), which strongly associated with reduced liver transplantation-free survival (log-rank P < .001). No taxonomic or functional differences were detected between patients with PSC with and without IBD. CONCLUSIONS: The gut microbiome in patients with PSC exhibits large functional differences compared with that in HCs, including microbial metabolism of essential nutrients. Alterations in related circulating metabolites associated with disease course, suggesting that microbial functions may be relevant for the disease process in PSC.

Genome-wide association study in 8,956 German individuals identifies influence of ABO histo-blood groups on gut microbiome.

The intestinal microbiome is implicated as an important modulating factor in multiple inflammatory1,2, neurologic3 and neoplastic diseases4. Recent genome-wide association studies yielded inconsistent, underpowered and rarely replicated results such that the role of human host genetics as a contributing factor to microbiome assembly and structure remains uncertain5-11. Nevertheless, twin studies clearly suggest host genetics as a driver of microbiome composition11. In a genome-wide association analysis of 8,956 German individuals, we identified 38 genetic loci to be associated with single bacteria and overall microbiome composition. Further analyses confirm the identified associations of ABO histo-blood groups and FUT2 secretor status with Bacteroides and Faecalibacterium spp. Mendelian randomization analysis suggests causative and protective effects of gut microbes, with clade-specific effects on inflammatory bowel disease. This holistic investigative approach of the host, its genetics and its associated microbial communities as a 'metaorganism' broaden our understanding of disease etiology, and emphasize the potential for implementing microbiota in disease treatment and management.

Postprandial factor VII activation does not increase plasma concentrations of prothrombin fragment 1 + 2 in patients with morbid obesity.

INTRODUCTION: Increased postprandial factor VII activation is observed after high-fat meals, but is not accompanied by thrombin formation in normal weight individuals. Obesity is associated with a higher circulating concentration of tissue factor (TF) and postprandial uptake of lipopolysaccharide (LPS), and this may increase thrombin formation after high-fat meals. We therefore compared postprandial effects of high-fat meals and low-fat meals on biomarkers of coagulation activation in patients with morbid obesity and investigated whether the response was associated with the gut bacteria composition. MATERIALS AND METHODS: A controlled cross-over study was conducted in obese patients (15 women, 5 men, mean BMI = 44.1 kg/m2), where high-fat meals (67 E% fat) and low-fat meals (16 E% fat) were served at 8:15 and 10:00 in a random order on two study days within one week. Blood samples were collected at 08:00 (fasting), 12:00, and 14:00 and analysed for triglycerides, activated FVII (FVIIa), TF, FVIIa-antithrombin (FVIIa-AT), prothrombin fragment 1 + 2 (F1+2), and TF pathway inhibitor (TFPI). The gut bacteria composition, measured as gram-negative bacteria and diversity, was analysed in faecal samples. RESULTS: Triglycerides, FVIIa, and FVIIa-AT increased significantly after high-fat meals, whereas F1 + 2 decreased equally and significantly over time after both meals. There were no significant changes in TF and TFPI over time. The postprandial changes in F1 + 2 and TFPI after high-fat meals were negatively correlated with diversity. CONCLUSIONS: Increased postprandial FVIIa is not accompanied by thrombin formation four hours after high-fat meals in patients with morbid obesity, possibly due to FVIIa-inhibition by AT.

Assessing the strength and sensitivity of the core microbiota approach on a highly diverse sponge reef.

Marine sponge reefs usually comprise a complex array of taxonomically different sponge species, many of these hosting highly diverse microbial communities. The number of microbial species known to occupy a given sponge ranges from tens to thousands, bringing numerous challenges to their analysis. One way to deal with such complexity is to use a core microbiota approach, in which only prevalent and abundant microbes are considered. Here we aimed to test the strength and sensitivity of the core microbiota approach by applying different core definitions to 20 host sponge species. Application of increasingly stringent relative abundance and/or percentage occurrence thresholds to qualify as part of the core microbiota decreased the number of 'core' OTUs and phyla and, consequently, changed both alpha- and beta-diversity patterns. Moreover, microbial co-occurrence patterns explored using correlation networks were also affected by the core microbiota definition. The application of stricter thresholds resulted in smaller and less compartmentalized networks, with different keystone species. These results highlight that the application of different core definitions to phylogenetically disparate host species can result in the drawing of markedly different conclusions. Consequently, we recommend to assess the effects of different core community definitions on the specific system of study before considering its application.

Antibiotics-induced monodominance of a novel gut bacterial order.

OBJECTIVE: The composition of the healthy human adult gut microbiome is relatively stable over prolonged periods, and representatives of the most highly abundant and prevalent species have been cultured and described. However, microbial abundances can change on perturbations, such as antibiotics intake, enabling the identification and characterisation of otherwise low abundant species. DESIGN: Analysing gut microbial time-series data, we used shotgun metagenomics to create strain level taxonomic and functional profiles. Community dynamics were modelled postintervention with a focus on conditionally rare taxa and previously unknown bacteria. RESULTS: In response to a commonly prescribed cephalosporin (ceftriaxone), we observe a strong compositional shift in one subject, in which a previously unknown species, UBorkfalki ceftriaxensis, was identified, blooming to 92% relative abundance. The genome assembly reveals that this species (1) belongs to a so far undescribed order of Firmicutes, (2) is ubiquitously present at low abundances in at least one third of adults, (3) is opportunistically growing, being ecologically similar to typical probiotic species and (4) is stably associated to healthy hosts as determined by single nucleotide variation analysis. It was the first coloniser after the antibiotic intervention that led to a long-lasting microbial community shift and likely permanent loss of nine commensals. CONCLUSION: The bloom of UB. ceftriaxensis and a subsequent one of Parabacteroides distasonis demonstrate the existence of monodominance community states in the gut. Our study points to an undiscovered wealth of low abundant but common taxa in the human gut and calls for more highly resolved longitudinal studies, in particular on ecosystem perturbations.

Marine Sponges as Chloroflexi Hot Spots: Genomic Insights and High-Resolution Visualization of an Abundant and Diverse Symbiotic Clade.

Members of the widespread bacterial phylum Chloroflexi can dominate high-microbial-abundance (HMA) sponge microbiomes. In the Sponge Microbiome Project, Chloroflexi sequences amounted to 20 to 30% of the total microbiome of certain HMA sponge genera with the classes/clades SAR202, Caldilineae, and Anaerolineae being the most prominent. We performed metagenomic and single-cell genomic analyses to elucidate the functional gene repertoire of Chloroflexi symbionts of Aplysina aerophoba. Eighteen draft genomes were reconstructed and placed into phylogenetic context of which six were investigated in detail. Common genomic features of Chloroflexi sponge symbionts were related to central energy and carbon converting pathways, amino acid and fatty acid metabolism, and respiration. Clade-specific metabolic features included a massively expanded genomic repertoire for carbohydrate degradation in Anaerolineae and Caldilineae genomes, but only amino acid utilization by SAR202. While Anaerolineae and Caldilineae import cofactors and vitamins, SAR202 genomes harbor genes encoding components involved in cofactor biosynthesis. A number of features relevant to symbiosis were further identified, including CRISPR-Cas systems, eukaryote-like repeat proteins, and secondary metabolite gene clusters. Chloroflexi symbionts were visualized in the sponge extracellular matrix at ultrastructural resolution by the fluorescence in situ hybridization-correlative light and electron microscopy (FISH-CLEM) method. Carbohydrate degradation potential was reported previously for "Candidatus Poribacteria" and SAUL, typical symbionts of HMA sponges, and we propose here that HMA sponge symbionts collectively engage in degradation of dissolved organic matter, both labile and recalcitrant. Thus, sponge microbes may not only provide nutrients to the sponge host, but they may also contribute to dissolved organic matter (DOM) recycling and primary productivity in reef ecosystems via a pathway termed the sponge loop. IMPORTANCE Chloroflexi represent a widespread, yet enigmatic bacterial phylum with few cultivated members. We used metagenomic and single-cell genomic approaches to characterize the functional gene repertoire of Chloroflexi symbionts in marine sponges. The results of this study suggest clade-specific metabolic specialization and that Chloroflexi symbionts have the genomic potential for dissolved organic matter (DOM) degradation from seawater. Considering the abundance and dominance of sponges in many benthic environments, we predict that the role of sponge symbionts in biogeochemical cycles is larger than previously thought.

The sponge microbiome project.

Marine sponges (phylum Porifera) are a diverse, phylogenetically deep-branching clade known for forming intimate partnerships with complex communities of microorganisms. To date, 16S rRNA gene sequencing studies have largely utilised different extraction and amplification methodologies to target the microbial communities of a limited number of sponge species, severely limiting comparative analyses of sponge microbial diversity and structure. Here, we provide an extensive and standardised dataset that will facilitate sponge microbiome comparisons across large spatial, temporal, and environmental scales. Samples from marine sponges (n = 3569 specimens), seawater (n = 370), marine sediments (n = 65) and other environments (n = 29) were collected from different locations across the globe. This dataset incorporates at least 268 different sponge species, including several yet unidentified taxa. The V4 region of the 16S rRNA gene was amplified and sequenced from extracted DNA using standardised procedures. Raw sequences (total of 1.1 billion sequences) were processed and clustered with (i) a standard protocol using QIIME closed-reference picking resulting in 39 543 operational taxonomic units (OTU) at 97% sequence identity, (ii) a de novo clustering using Mothur resulting in 518 246 OTUs, and (iii) a new high-resolution Deblur protocol resulting in 83 908 unique bacterial sequences. Abundance tables, representative sequences, taxonomic classifications, and metadata are provided. This dataset represents a comprehensive resource of sponge-associated microbial communities based on 16S rRNA gene sequences that can be used to address overarching hypotheses regarding host-associated prokaryotes, including host specificity, convergent evolution, environmental drivers of microbiome structure, and the sponge-associated rare biosphere.

Predicting the HMA-LMA Status in Marine Sponges by Machine Learning.

The dichotomy between high microbial abundance (HMA) and low microbial abundance (LMA) sponges has been observed in sponge-microbe symbiosis, although the extent of this pattern remains poorly unknown. We characterized the differences between the microbiomes of HMA (n = 19) and LMA (n = 17) sponges (575 specimens) present in the Sponge Microbiome Project. HMA sponges were associated with richer and more diverse microbiomes than LMA sponges, as indicated by the comparison of alpha diversity metrics. Microbial community structures differed between HMA and LMA sponges considering Operational Taxonomic Units (OTU) abundances and across microbial taxonomic levels, from phylum to species. The largest proportion of microbiome variation was explained by the host identity. Several phyla, classes, and OTUs were found differentially abundant in either group, which were considered "HMA indicators" and "LMA indicators." Machine learning algorithms (classifiers) were trained to predict the HMA-LMA status of sponges. Among nine different classifiers, higher performances were achieved by Random Forest trained with phylum and class abundances. Random Forest with optimized parameters predicted the HMA-LMA status of additional 135 sponge species (1,232 specimens) without a priori knowledge. These sponges were grouped in four clusters, from which the largest two were composed of species consistently predicted as HMA (n = 44) and LMA (n = 74). In summary, our analyses shown distinct features of the microbial communities associated with HMA and LMA sponges. The prediction of the HMA-LMA status based on the microbiome profiles of sponges demonstrates the application of machine learning to explore patterns of host-associated microbial communities.

Integrated metabolism in sponge-microbe symbiosis revealed by genome-centered metatranscriptomics.

Despite an increased understanding of functions in sponge microbiomes, the interactions among the symbionts and between symbionts and host are not well characterized. Here we reconstructed the metabolic interactions within the sponge Cymbastela concentrica microbiome in the context of functional features of symbiotic diatoms and the host. Three genome bins (CcPhy, CcNi and CcThau) were recovered from metagenomic data of C. concentrica, belonging to the proteobacterial family Phyllobacteriaceae, the Nitrospira genus and the thaumarchaeal order Nitrosopumilales. Gene expression was estimated by mapping C. concentrica metatranscriptomic reads. Our analyses indicated that CcPhy is heterotrophic, while CcNi and CcThau are chemolithoautotrophs. CcPhy expressed many transporters for the acquisition of dissolved organic compounds, likely available through the sponge's filtration activity and symbiotic carbon fixation. Coupled nitrification by CcThau and CcNi was reconstructed, supported by the observed close proximity of the cells in fluorescence in situ hybridization. CcPhy facultative anaerobic respiration and assimilation by diatoms may consume the resulting nitrate. Transcriptional analysis of diatom and sponge functions indicated that these organisms are likely sources of organic compounds, for example, creatine/creatinine and dissolved organic carbon, for other members of the symbiosis. Our results suggest that organic nitrogen compounds, for example, creatine, creatinine, urea and cyanate, fuel the nitrogen cycle within the sponge. This study provides an unprecedented view of the metabolic interactions within sponge-microbe symbiosis, bridging the gap between cell- and community-level knowledge.

An Enrichment of CRISPR and Other Defense-Related Features in Marine Sponge-Associated Microbial Metagenomes.

Many marine sponges are populated by dense and taxonomically diverse microbial consortia. We employed a metagenomics approach to unravel the differences in the functional gene repertoire among three Mediterranean sponge species, Petrosia ficiformis, Sarcotragus foetidus, Aplysina aerophoba and seawater. Different signatures were observed between sponge and seawater metagenomes with regard to microbial community composition, GC content, and estimated bacterial genome size. Our analysis showed further a pronounced repertoire for defense systems in sponge metagenomes. Specifically, clustered regularly interspaced short palindromic repeats, restriction modification, DNA phosphorothioation and phage growth limitation systems were enriched in sponge metagenomes. These data suggest that defense is an important functional trait for an existence within sponges that requires mechanisms to defend against foreign DNA from microorganisms and viruses. This study contributes to an understanding of the evolutionary arms race between viruses/phages and bacterial genomes and it sheds light on the bacterial defenses that have evolved in the context of the sponge holobiont.