Host Species and Environment Shape the Gut Microbiota of Cohabiting Marine Bivalves

Pacific oysters (Crassostrea gigas) and Mediterranean mussels (Mytilus galloprovincialis) are commercially important marine bivalves that frequently coexist and have overlapping feeding ecologies. Like other invertebrates, their gut microbiota is thought to play an important role in supporting their health and nutrition. Yet, little is known regarding the role of the host and environment in driving these communities. Here, bacterial assemblages were surveyed from seawater and gut aspirates of farmed C. gigas and co-occurring wild M. galloprovincialis in summer and winter using Illumina 16S rRNA gene sequencing. Unlike seawater, which was dominated by Pseudomonadata, bivalve samples largely consisted of Mycoplasmatota (Mollicutes) and accounted for >50% of the total OTU abundance. Despite large numbers of common (core) bacterial taxa, bivalve-specific species (OTUs) were also evident and predominantly associated with Mycoplasmataceae (notably Mycoplasma). An increase in diversity (though with varied taxonomic evenness) was observed in winter for both bivalves and was associated with changes in the abundance of core and bivalve-specific taxa, including several representing host-associated and environmental (free-living or particle-diet associated) organisms. Our findings highlight the contribution of the environment and the host in defining the composition of the gut microbiota in cohabiting, intergeneric bivalve populations.

Variations in Microbial Diversity and Metabolite Profiles of the Tropical Marine Sponge Xestospongia muta with Season and Depth

Xestospongia muta is among the most emblematic sponge species inhabiting coral reefs of the Caribbean Sea. Besides being the largest sponge species growing in the Caribbean, it is also known to produce secondary metabolites. This study aimed to assess the effect of depth and season on the symbiotic bacterial dynamics and major metabolite profiles of specimens of X. muta thriving in a tropical marine biome (Portobelo Bay, Panamá), which allow us to determine whether variability patterns are similar to those reported for subtropical latitudes. The bacterial assemblages were characterized using Illumina deep-sequencing and metabolomic profiles using UHPLC-DAD-ELSD from five depths (ranging 9-28 m) across two seasons (spring and autumn). Diverse symbiotic communities, representing 24 phyla with a predominance of Proteobacteria and Chloroflexi, were found. Although several thousands of OTUs were determined, most of them belong to the rare biosphere and only 23 to a core community. There was a significant difference between the structure of the microbial communities in respect to season (autumn to spring), with a further significant difference between depths only in autumn. This was partially mirrored in the metabolome profile, where the overall metabolite composition did not differ between seasons, but a significant depth gradient was observed in autumn. At the phyla level, Cyanobacteria, Firmicutes, Actinobacteria, and Spirochaete showed a mild-moderate correlation with the metabolome profile. The metabolomic profiles were mainly characterized by known brominated polyunsaturated fatty acids. This work presents findings about the composition and dynamics of the microbial assemblages of X. muta expanding and confirming current knowledge about its remarkable diversity and geographic variability as observed in this tropical marine biome.

The active bacterial assemblages of the upper GI tract in individuals with and without Helicobacter infection

OBJECTIVE

Patients infected with Helicobacter pylori develop chronic gastritis with a subgroup progressing to further complications. The role of microbiota from the oral cavity swallowed with saliva and either transiting the stomach or persisting in the gastric mucosa is uncertain. It is also not known whether the bacterial community differs in luminal and mucosal niches. A key question is whether H. pylori influences the bacterial communities of gastroduodenal niches.

DESIGN

Saliva, gastric and duodenal aspirates as well as gastric and duodenal biopsies were collected during oesophagogastroduodenoscopy from 24 patients (m:9, f:15, mean age 52.2±SD 14.5 years). RNA was extracted and the V1-V2 region of the retrotranscribed bacterial 16S rRNA amplified and sequenced on the Illumina MiSeq platform.

RESULTS

Overall, 687 bacterial phylotypes that belonged to 95 genera and 11 phyla were observed. Each individual comprised a unique microbiota composition that was consistent across the different niches. However, the stomach fluid enriched for specific microbiota components. Helicobacter spp were shown to dominate the mucosa-associated community in the stomach, and to significantly influence duodenal and oral communities.

CONCLUSIONS

The detailed analysis of the active global bacterial communities from the five distinct sites of the upper GI tract allowed for the first time the differentiation between host effects and the influence of sampling region on the bacterial community. The influence of Helicobacter spp on the global community structures is striking.

The Inner Workings of the Outer Surface: Skin and Gill Microbiota as Indicators of Changing Gut Health in Yellowtail Kingfish

The mucosal surfaces and associated microbiota of fish are an important primary barrier and provide the first line of defense against potential pathogens. An understanding of the skin and gill microbial assemblages and the factors which drive their composition may provide useful insights into the broad dynamics of fish host–microbial relationships, and may reveal underlying changes in health status. This is particularly pertinent to cultivated systems whereby various stressors may led to conditions (like enteritis) which impinge on productivity. As an economically important species, we assessed whether the outer-surface bacterial communities reflect a change in gut health status of cultivated Yellowtail Kingfish (Seriola lalandi). Active bacterial assemblages were surveyed from RNA extracts from swabs of the skin and gills by constructing Illumina 16S rRNA gene amplicon libraries. Proteobacteria and Bacteroidetes were predominant in both the skin and gills, with enrichment of key β-proteobacteria in the gills (Nitrosomonadales and Ferrovales). Fish exhibiting early stage chronic lymphocytic enteritis comprised markedly different global bacterial assemblages compared to those deemed healthy and exhibiting late stages of the disease. This corresponded to an overall loss of diversity and enrichment of Proteobacteria and Actinobacteria, particularly in the gills. In contrast, bacterial assemblages of fish with late stage enteritis were generally similar to those of healthy individuals, though with some distinct taxa. In conclusion, gut health status is an important factor which defines the skin and gill bacterial assemblages of fish and likely reflects changes in immune states and barrier systems during the early onset of conditions like enteritis. This study represents the first to investigate the microbiota of the outer mucosal surfaces of fish in response to underlying chronic gut enteritis, revealing potential biomarkers for assessing fish health in commercial aquaculture systems.

Using in situ hybridization to expand the daily egg production method to new fish species

The capacity to reliably identify fish eggs is critical in the application of the daily egg production method (DEPM) to estimate biomass of commercially important species. This application has largely been confined to species that have easily identifiable eggs. Various molecular strategies have been used to extend the DEPM to a broader range of species, with recent approaches like in situ hybridization (ISH) that preserves the integrity of whole eggs, embryos or larvae recommended as a suitable alternative over destructive procedures like PCR. Here, we designed and validated an ISH approach for the identification of whole eggs and larvae from Snapper (Chrysophrys auratus) from environmental samples using the mitochondrial 16S rRNA gene as a target for specific horseradish peroxidase (HRP)-conjugated oligonucleotide probes. This colorimetric assay allowed the highly specific detection of positive hybridization signals from intact C. auratus larvae and eggs from mixed-species samples comprising closely related taxa. Furthermore, evaluation of whole eggs across a range of developmental stages revealed the sensitivity of the approach for discerning early stages, thereby guiding staging and the identification of otherwise indistinguishable eggs from environmental samples. This approach represents a major advance from current molecular-based strategies as it is nondestructive and allows for the simultaneous identification and staging of fish eggs (and larvae). The resultant 100% egg identification certainty we have achieved allows the DEPM to be applied to a wider array of fish species and is particularly applicable to species in areas where morphologically similar eggs are being spawned at the same time.

Exploring the bacterial assemblages along the human nasal passage

The human nasal passage, from the anterior nares through the nasal vestibule to the nasal cavities, is an important habitat for opportunistic pathogens and commensals alike. This work sampled four different anatomical regions within the human nasal passage across a large cohort of individuals (n = 79) comprising individuals suffering from chronic nasal inflammation clinically known as chronic rhinosinusitis (CRS) and individuals not suffering from inflammation (CRS-free). While individuals had their own unique bacterial fingerprint that was consistent across the anatomical regions, these bacterial fingerprints formed into distinct delineated groups comprising core bacterial members, which were consistent across all four swabbed anatomical regions irrespective of health status. The most significant observed pattern was the difference between the global bacterial profiles of swabbed and tissue biopsy samples from the same individuals, being also consistent across different anatomical regions. Importantly, no statistically significant differences could be observed concerning the global bacterial communities, any of the bacterial species or the range of diversity indices used to compare between CRS and CRS-free individuals, and between two CRS phenotypes (without nasal polyps and with nasal polyps). Thus, the role of bacteria in the pathogenesis of sinusitis remains uncertain.

Functional soil metagenomics: elucidation of polycyclic aromatic hydrocarbon degradation potential following 12 years of in situ bioremediation

A culture-independent function-based screening approach was used to assess the microbial aerobic catabolome for polycyclic aromatic hydrocarbons degradation of a soil subjected to 12 years of in situ bioremediation. A total of 422 750 fosmid clones were screened for key aromatic ring-cleavage activities using 2,3-dihydroxybiphenyl as substrate. Most of the genes encoding ring-cleavage enzymes on the 768 retrieved positive fosmids could not be identified using primer-based approaches and, thus, 205 fosmid inserts were sequenced. Nearly two hundred extradiol dioxygenase encoding genes of three different superfamilies could be identified. Additional key genes of aromatic metabolic pathways were identified, including a high abundance of Rieske non-heme iron oxygenases that provided detailed information on enzymes activating aromatic compounds and enzymes involved in activation of the side chain of methylsubstituted aromatics. The gained insights indicated a complex microbial network acting at the site under study, which comprises organisms similar to recently identified Immundisolibacter cernigliae TR3.2 and Rugosibacter aromaticivorans Ca6 and underlined the great potential of an approach that combines an activity-screening, a cost-effective high-throughput sequencing of fosmid clones and a phylogenomic-routed and manually curated database to carefully identify key proteins dedicated to aerobic degradation of aromatic compounds.

Exploring the transcriptome of Staphylococcus aureus in its natural niche

Staphylococcus aureus is an important human pathogen and commensal, where the human nose is the predominant reservoir. To better understand its behavior in this environmental niche, RNA was extracted from the anterior nares of three documented S. aureus carriers and the metatranscriptome analyzed by RNAseq. In addition, the in vivo transcriptomes were compared to previously published transcriptomes of two in vitro grown S. aureus strains. None of the in vitro conditions, even growth in medium resembling the anterior nares environment, mimicked in vivo conditions. Survival in the nose was strongly controlled by the limitation of iron and evident by the expression of iron acquisition systems. S. aureus populations in different individuals clearly experience different environmental stresses, which they attempt to overcome by the expression of compatible solute biosynthetic pathways, changes in their cell wall composition and synthesis of general stress proteins. Moreover, the expression of adhesins was also important for colonization of the anterior nares. However, different S. aureus strains also showed different in vivo behavior. The assessment of general in vivo expression patterns and commonalities between different S. aureus strains will in the future result in new knowledge based strategies for controlling colonization.

Co-occurrence Analysis of Microbial Taxa in the Atlantic Ocean Reveals High Connectivity in the Free-Living Bacterioplankton

We determined the taxonomic composition of the bacterioplankton of the epipelagic zone of the Atlantic Ocean along a latitudinal transect (51°S–47°N) using Illumina sequencing of the V5-V6 region of the 16S rRNA gene and inferred co-occurrence networks. Bacterioplankon community composition was distinct for Longhurstian provinces and water depth. Free-living microbial communities (between 0.22 and 3 μm) were dominated by highly abundant and ubiquitous taxa with streamlined genomes (e.g., SAR11, SAR86, OM1, Prochlorococcus) and could clearly be separated from particle-associated communities which were dominated by Bacteroidetes, Planktomycetes, Verrucomicrobia, and Roseobacters. From a total of 369 different communities we then inferred co-occurrence networks for each size fraction and depth layer of the plankton between bacteria and between bacteria and phototrophic micro-eukaryotes. The inferred networks showed a reduction of edges in the deepest layer of the photic zone. Networks comprised of free-living bacteria had a larger amount of connections per OTU when compared to the particle associated communities throughout the water column. Negative correlations accounted for roughly one third of the total edges in the free-living communities at all depths, while they decreased with depth in the particle associated communities where they amounted for roughly 10% of the total in the last part of the epipelagic zone. Co-occurrence networks of bacteria with phototrophic micro-eukaryotes were not taxon-specific, and dominated by mutual exclusion (~60%). The data show a high degree of specialization to micro-environments in the water column and highlight the importance of interdependencies particularly between free-living bacteria in the upper layers of the epipelagic zone.

Bacterioplankton Biogeography of the Atlantic Ocean: A Case Study of the Distance-Decay Relationship

In order to determine the influence of geographical distance, depth, and Longhurstian province on bacterial community composition and compare it with the composition of photosynthetic micro-eukaryote communities, 382 samples from a depth-resolved latitudinal transect (51°S–47°N) from the epipelagic zone of the Atlantic ocean were analyzed by Illumina amplicon sequencing. In the upper 100 m of the ocean, community similarity decreased toward the equator for 6000 km, but subsequently increased again, reaching similarity values of 40–60% for samples that were separated by ~12,000 km, resulting in a U-shaped distance-decay curve. We conclude that adaptation to local conditions can override the linear distance-decay relationship in the upper epipelagial of the Atlantic Ocean which is apparently not restrained by barriers to dispersal, since the same taxa were shared between the most distant communities. The six Longhurstian provinces covered by the transect were comprised of distinct microbial communities; ~30% of variation in community composition could be explained by province. Bacterial communities belonging to the deeper layer of the epipelagic zone (140–200 m) lacked a distance-decay relationship altogether and showed little provincialism. Interestingly, those biogeographical patterns were consistently found for bacteria from three different size fractions of the plankton with different taxonomic composition, indicating conserved underlying mechanisms. Analysis of the chloroplast 16S rRNA gene sequences revealed that phytoplankton composition was strongly correlated with both free-living and particle associated bacterial community composition (R between 0.51 and 0.62, p < 0.002). The data show that biogeographical patterns commonly found in macroecology do not hold for marine bacterioplankton, most likely because dispersal and evolution occur at drastically different rates in bacteria.

Application of a Novel "Pan-Genome"-Based Strategy for Assigning RNAseq Transcript Reads to Staphylococcus aureus Strains

Understanding the behaviour of opportunistic pathogens such as Staphylococcus aureus in their natural human niche holds great medical interest. With the development of sensitive molecular methods and deep-sequencing technology, it is now possible to robustly assess the global transcriptome of bacterial species in their human habitat. However, as the genomes of the colonizing strains are often not available compiling the pan-genome for the species of interest may provide an effective method to reliably and rapidly compile the transcriptome of a bacterial species. The pan-genome of S. aureus and its associated core and accessory components were compiled based on 25 genomes and comprises a total of 65,557 proteins clustering into 4,198 Orthologous Groups (OGs). The generated gene catalogue was used to assign RNAseq-derived sequence reads to S. aureus in a variety of in vitro and in vivo samples. In all cases, the number of reads that could be assigned to S. aureus was greater using the OG database than using a reference genome. Growth of two S. aureus strains in synthetic nasal medium confirmed that both strains experienced strong iron starvation. Traits such as purine metabolism appeared to be more affected in a typical nasal colonizer than in a strain representative of the S. aureus USA300 lineage. Mapping sequencing reads from a metatranscriptome generated from the human anterior nares allowed the identification of genes highly expressed by S. aureus in vivo. The OG database generated in this study represents a useful tool to obtain a snapshot of the functional attributes of S. aureus under different in vitro and in vivo conditions. The approach proved to be advantageous to assign sequencing reads to bacterial strains when RNAseq data is derived from samples where strain information and/or the corresponding genome/s are unavailable.

High-resolution taxonomic profiling of the subgingival microbiome for biomarker discovery and periodontitis diagnosis

The oral microbiome plays a key role for caries, periodontitis, and systemic diseases. A method for rapid, high-resolution, robust taxonomic profiling of subgingival bacterial communities for early detection of periodontitis biomarkers would therefore be a useful tool for individualized medicine. Here, we used Illumina sequencing of the V1-V2 and V5-V6 hypervariable regions of the 16S rRNA gene. A sample stratification pipeline was developed in a pilot study of 19 individuals, 9 of whom had been diagnosed with chronic periodontitis. Five hundred twenty-three operational taxonomic units (OTUs) were obtained from the V1-V2 region and 432 from the V5-V6 region. Key periodontal pathogens like Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia could be identified at the species level with both primer sets. Principal coordinate analysis identified two outliers that were consistently independent of the hypervariable region and method of DNA extraction used. The linear discriminant analysis (LDA) effect size algorithm (LEfSe) identified 80 OTU-level biomarkers of periodontitis and 17 of health. Health- and periodontitis-related clusters of OTUs were identified using a connectivity analysis, and the results confirmed previous studies with several thousands of samples. A machine learning algorithm was developed which was trained on all but one sample and then predicted the diagnosis of the left-out sample (jackknife method). Using a combination of the 10 best biomarkers, 15 of 17 samples were correctly diagnosed. Training the algorithm on time-resolved community profiles might provide a highly sensitive tool to detect the onset of periodontitis.

Comparing the anterior nare bacterial community of two discrete human populations using Illumina amplicon sequencing

The anterior nares are an important reservoir for opportunistic pathogens and commensal microorganisms. A barcoded Illumina paired-end sequencing method targeting the 16S ribosomal RNA V1-2 hypervariable region was developed to compare the bacterial diversity of the anterior nares across distinct human populations (volunteers from Germany vs a Babongo Pygmy tribe, Africa). Of the 251 phylotypes detected, 231 could be classified to the genus level and 109 to the species level, including the unambiguous identification of the ubiquitous Staphylococcus aureus and Moraxella catarrhalis. The global bacterial community of both adult populations revealed that they shared 85% of the phylotypes, suggesting that our global bacterial communities have likely been with us for thousands of years. Of the 34 phylotypes unique to the non-westernized population, most were related to members within the suborder Micrococcineae. There was an even more overwelming distinction between children and adults of the same population, suggesting a progression of a childhood community of high-diversity comprising species of Moraxellaceae and Streptococcaceae to an adult community of lower diversity comprising species of Propionibacteriaceae, Clostridiales Incertae Sedis XI, Corynebacteriaceae and Staphylococcaceae. Thus, age was a stronger factor for accounting for differing bacterial assemblages than the origin of the human population sampled.

Comparative evaluation of establishing a human gut microbial community within rodent models

The structure of the human gut microbial community is determined by host genetics and environmental factors, where alterations in its structure have been associated with the onset of different diseases. Establishing a defined human gut microbial community within inbred rodent models provides a means to study microbial-related pathologies, however, an in-depth comparison of the established human gut microbiota in the different models is lacking. We compared the efficiency of establishing the bacterial component of a defined human microbial community within germ-free (GF) rats, GF mice, and antibiotic-treated specific pathogen-free mice. Remarkable differences were observed between the different rodent models. While the majority of abundant human-donor bacterial phylotypes were established in the GF rats, only a subset was present in the GF mice. Despite the fact that members of the phylum Bacteriodetes were well established in all rodent models, mice enriched for phylotypes related to species of Bacteroides. In contrary to the efficiency of Clostridiales to populate the GF rat in relative proportions to that of the human-donor, members of Clostridia cluster IV only poorly colonize the mouse gut. Thus, the genetic background of the different recipient rodent systems (that is, rats and mice) strongly influences the nature of the populating human gut microbiota, determining each model's biological suitability.

Strain-specific colitis susceptibility in IL10-deficient mice depends on complex gut microbiota-host interactions

BACKGROUND

Colitis susceptibility in Il10(-/-) mice depends on genetic background and microbiota composition. A major genetic locus mediating colitis susceptibility, Cdcs1, was transferred from susceptible C3Bir-Il10(-/-) to resistant B6-Il10(-/-) mice, resulting in susceptible congenic BC-R3-Il10(-/-) mice. The aim of this study was to determine the impact of microbiota on this differential colitis susceptibility using a Helicobacter hepaticus infection model.

METHODS

Parental C3Bir-Il10(-/-) , B6-Il10(-/-) , and congenic BC-R3-Il10(-/-) mice were inoculated with H. hepaticus and analyzed for inflammation. In parental Il10(-/-) mice, microbiota composition was determined by terminal restriction fragment length polymorphism (T-RFLP) and quantitative polymerase chain reaction (qPCR).

RESULTS

Most severe inflammation was observed in C3Bir-Il10(-/-) in the cecum, in BC-R3-Il10(-/-) in cecum and colon, and, unexpectedly, in B6-Il10(-/-) in the colon. C3Bir-Il10(-/-) and BC-R3-Il10(-/-) secreted significantly more interferon-gamma (IFNγ) and interleukin (IL)17 than B6-Il10(-/-) . T-RFLP analyses in C3Bir-Il10(-/-) and B6-Il10(-/-) mice revealed 1) a significant impact of H. hepaticus infection on species richness and diversity, and 2) strain differences in microbiota composition only after H. hepaticus infection. qPCR revealed higher numbers of Clostridia leptum and Bacteroides spp. in the cecum of infected C3Bir-Il10(-/-) mice, and Lactobacillus spp. in B6-Il10(-/-) mice.

CONCLUSIONS

Cdcs1 modifies the response to H. hepaticus infection. However, this infection alone does not reflect the original response to a complex colitogenic biota. H. hepaticus-induced inflammation altered intestinal microbiota in a mouse strain-specific manner. Bacteroides spp. became more abundant in susceptible C3Bir-Il10(-/-) , lactobacilli in B6-Il10(-/-) mice. Therefore, both host immune response and differential compositional changes of microbiota play a role in strain-specific colitis susceptibility in Il10(-/-) mice.

Validating T-RFLP as a sensitive and high-throughput approach to assess bacterial diversity patterns in human anterior nares

While recent works aimed to thoroughly characterize the bacterial community of the human anterior nares of a few candidates, this work sought to analyse a greater cross-section by sampling 100 volunteers. After optimizing and validating the method of terminal restriction fragment length polymorphism against six previously pyrosequenced samples, abundant species could be discriminated and their relative abundances measured in a high-throughput manner. The 100 volunteers could be statistically clustered into 12 groups, where two-thirds of volunteers shared more than 40% similarity in respect to their bacterial community structure, while the remaining third clustered into smaller groups being dominated by Dolosigranulum pigrum, Moraxella spp. or Staphylococcus aureus. Moraxella spp. was present predominantly in women rather than in men. The use of network analysis charting bacterial ecological co-occurrences revealed new evidence of likely positive associations between some core human nasal species. So, in the age of post 'omics' and 'deep sequencing', there is still a place for these well-tried and well-tested methods that can offer a rapid, reproducible and economical alternative, whereby also yielding valuable new information.

The temporal dynamics of bacterial communities across human anterior nares

The anterior nares (nostrils), are an important niche for bacterial colonization by both commensals and opportunistic pathogens. Here the temporal dynamics and variation of the global nasal bacterial community across 25 healthy volunteers was evaluated over 15 months. Overall, there was a global seasonal shift in bacterial community structure. Such a temporal shift was also strongly evident in the abundances of species such as Propionibacterium acnes and Staphylococcus epidermidis. However, such species dynamics over time was also inter-individual-dependent, and both individuals with highly stable communities and those with highly flexible communities could be defined. Even though the bacterial community of individual volunteers was thus generally variable over time and permanent carriage of a given species was seldomly observed, various species previously defined as constituting the core bacterial community could be identified as persistent in a subset of the volunteers suggesting that these same species also constitute to a 'temporal' core community.

In-vivo expression profiling of Pseudomonas aeruginosa infections reveals niche-specific and strain-independent transcriptional programs

Pseudomonas aeruginosa is a threatening, opportunistic pathogen causing disease in immunocompromised individuals. The hallmark of P. aeruginosa virulence is its multi-factorial and combinatorial nature. It renders such bacteria infectious for many organisms and it is often resistant to antibiotics. To gain insights into the physiology of P. aeruginosa during infection, we assessed the transcriptional programs of three different P. aeruginosa strains directly after isolation from burn wounds of humans. We compared the programs to those of the same strains using two infection models: a plant model, which consisted of the infection of the midrib of lettuce leaves, and a murine tumor model, which was obtained by infection of mice with an induced tumor in the abdomen. All control conditions of P. aeruginosa cells growing in suspension and as a biofilm were added to the analysis. We found that these different P. aeruginosa strains express a pool of distinct genetic traits that are activated under particular infection conditions regardless of their genetic variability. The knowledge herein generated will advance our understanding of P. aeruginosa virulence and provide valuable cues for the definition of prospective targets to develop novel intervention strategies.