Defining marine bacterioplankton community assembly rules by contrasting the importance of environmental determinants and biotic interactions.

Bacterioplankton communities govern marine productivity and biogeochemical cycling, yet drivers of bacterioplankton assembly remain unclear. Here, we contrast the relative contribution of deterministic processes (environmental factors and biotic interactions) in driving temporal dynamics of bacterioplankton diversity at three different oceanographic time series locations, spanning 15° of latitude, which are each characterized by different environmental conditions and varying degrees of seasonality. Monthly surface samples (5.5 years) were analysed using 16S rRNA amplicon sequencing. The high- and mid-latitude sites of Maria Island and Port Hacking were characterized by high and intermediate levels of environmental heterogeneity, respectively, with both alpha diversity (72%; 24% of total variation) and beta diversity (32%; 30%) patterns within bacterioplankton assemblages explained by day length, ammonium, and mixed layer depth. In contrast, North Stradbroke Island, a sub-tropical location where environmental conditions are less variable, interspecific interactions were of increased importance in structuring bacterioplankton diversity (alpha: 33%; beta: 26%) with environment only contributing 11% and 13% to predicting diversity, respectively. Our results demonstrate that bacterioplankton diversity is the result of both deterministic environmental and biotic processes and that the importance of these different deterministic processes varies, potential in response to environmental heterogeneity.

The Epidermal Microbiome Within an Aggregation of Leopard Sharks (Triakis semifasciata) Has Taxonomic Flexibility with Gene Functional Stability Across Three Time-points.

The epidermis of Chondrichthyan fishes consists of dermal denticles with production of minimal but protein-rich mucus that collectively, influence the attachment and biofilm development of microbes, facilitating a unique epidermal microbiome. Here, we use metagenomics to provide the taxonomic and functional characterization of the epidermal microbiome of the Triakis semifasciata (leopard shark) at three time-points collected across 4 years to identify links between microbial groups and host metabolism. Our aims include (1) describing the variation of microbiome taxa over time and identifying recurrent microbiome members (present across all time-points); (2) investigating the relationship between the recurrent and flexible taxa (those which are not found consistently across time-points); (3) describing the functional compositions of the microbiome which may suggest links with the host metabolism; and (4) identifying whether metabolic processes are shared across microbial genera or are unique to specific taxa. Microbial members of the microbiome showed high similarity between all individuals (Bray-Curtis similarity index = 82.7, where 0 = no overlap, 100 = total overlap) with the relative abundance of those members varying across sampling time-points, suggesting flexibility of taxa in the microbiome. One hundred and eighty-eight genera were identified as recurrent, including Pseudomonas, Erythrobacter, Alcanivorax, Marinobacter, and Sphingopxis being consistently abundant across time-points, while Limnobacter and Xyella exhibited switching patterns with high relative abundance in 2013, Sphingobium and Sphingomona in 2015, and Altermonas, Leeuwenhoekiella, Gramella, and Maribacter in 2017. Of the 188 genera identified as recurrent, the top 19 relatively abundant genera formed three recurrent groups. The microbiome also displayed high functional similarity between individuals (Bray-Curtis similarity index = 97.6) with gene function composition remaining consistent across all time-points. These results show that while the presence of microbial genera exhibits consistency across time-points, their abundances do fluctuate. Microbial functions however remain stable across time-points; thus, we suggest the leopard shark microbiomes exhibit functional redundancy. We show coexistence of microbes hosted in elasmobranch microbiomes that encode genes involved in utilizing nitrogen, but not fixing nitrogen, degrading urea, and resistant to heavy metal.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes.

BACKGROUND: Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown. RESULTS: Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes. CONCLUSIONS: This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

Rhodoliths holobionts in a changing ocean: host-microbes interactions mediate coralline algae resilience under ocean acidification.

BACKGROUND: Life in the ocean will increasingly have to contend with a complex matrix of concurrent shifts in environmental properties that impact their physiology and control their life histories. Rhodoliths are coralline red algae (Corallinales, Rhodophyta) that are photosynthesizers, calcifiers, and ecosystem engineers and therefore represent important targets for ocean acidification (OA) research. Here, we exposed live rhodoliths to near-future OA conditions to investigate responses in their photosynthetic capacity, calcium carbonate production, and associated microbiome using carbon uptake, decalcification assays, and whole genome shotgun sequencing metagenomic analysis, respectively. The results from our live rhodolith assays were compared to similar manipulations on dead rhodolith (calcareous skeleton) biofilms and water column microbial communities, thereby enabling the assessment of host-microbiome interaction under climate-driven environmental perturbations. RESULTS: Under high pCO2 conditions, live rhodoliths exhibited positive physiological responses, i.e. increased photosynthetic activity, and no calcium carbonate biomass loss over time. Further, whereas the microbiome associated with live rhodoliths remained stable and resembled a healthy holobiont, the microbial community associated with the water column changed after exposure to elevated pCO2. CONCLUSIONS: Our results suggest that a tightly regulated microbial-host interaction, as evidenced by the stability of the rhodolith microbiome recorded here under OA-like conditions, is important for host resilience to environmental stress. This study extends the scarce comprehension of microbes associated with rhodolith beds and their reaction to increased pCO2, providing a more comprehensive approach to OA studies by assessing the host holobiont.

Who bit the boat? New DNA collection and genomic methods enable species identification in suspected shark-related incidents.

Species identification following shark-related incidents is critical for effective incident management and for collecting data to inform shark-bite mitigation strategies. Witness statements are not always reliable, and species identification is often ambiguous or missing. Alternative methods for species identification include morphological assessments of bite marks, analysis of collected teeth at the scene of the incident, and genetic approaches. However, access to appropriate collection media and robust genetic assays have limited the use of genetic technologies. Here, we present a case study that facilitated a unique opportunity to compare the effectiveness of medical gauze readily available in first-aid kits, and forensic-grade swabs in collecting genetic material for shark-species identification. Sterile medical gauze and forensic-grade swabs were used to collect transfer DNA from the bite margins on a bitten surf ski which were compared to a piece of shark tissue embedded along the bite margin. Witness accounts and the characteristics of the bite mark impressions inferred the involvement of a Carcharodon carcharias (white shark). The morphology of a tooth found on the boat that picked up the surf ski, however, suggested it belonged to an Orectolobus spp. (wobbegong). Genetic analysis of DNA transferred from the shark to the surf ski included the application of a broad-target nested PCR assay followed by Sanger sequencing, with white shark contribution to the 'total sample DNA' determined with a species-specific qPCR assay. The results of the genetic analyses were congruent between sampling methods with respect to species identification and the level of activity inferred by the donor-specific DNA contribution. These data also supported the inferences drawn from the bite mark morphology. DNA from the recovered tooth was PCR amplified with a wobbegong-specific primer pair designed for this study to corroborate the tooth's morphological identification. Following the confirmation of gauze used for sampling in the case study event, two additional isolated incidents occurred and were sampled in situ using gauze, as typically found in a first-aid kit, by external personnel. DNA extracted from these gauze samples resulted in the identification of a white shark as the donor of the DNA collected from the bite marks in both instances. This study, involving three incidents separated by time and location, represents the seminal application of gauze as a sampling media after critical human-shark interactions and strongly supports the practical implementation of these methods in the field.

Bioenergetic mapping of 'healthy microbiomes' via compound processing potential imprinted in gut and soil metagenomes.

Despite mounting evidence of their importance in human health and ecosystem functioning, the definition and measurement of 'healthy microbiomes' remain unclear. More advanced knowledge exists on health associations for compounds used or produced by microbes. Environmental microbiome exposures (especially via soils) also help shape, and may supplement, the functional capacity of human microbiomes. Given the synchronous interaction between microbes, their feedstocks, and micro-environments, with functional genes facilitating chemical transformations, our objective was to examine microbiomes in terms of their capacity to process compounds relevant to human health. Here we integrate functional genomics and biochemistry frameworks to derive new quantitative measures of in silico potential for human gut and environmental soil metagenomes to process a panel of major compound classes (e.g., lipids, carbohydrates) and selected biomolecules (e.g., vitamins, short-chain fatty acids) linked to human health. Metagenome functional potential profile data were translated into a universal compound mapping 'landscape' based on bioenergetic van Krevelen mapping of function-level meta-compounds and corresponding functional relative abundances, reflecting imprinted genetic capacity of microbiomes to metabolize an array of different compounds. We show that measures of 'compound processing potential' associated with human health and disease (examining atherosclerotic cardiovascular disease, colorectal cancer, type 2 diabetes and anxious-depressive behavior case studies), and displayed seemingly predictable shifts along gradients of ecological disturbance in plant-soil ecosystems (three case studies). Ecosystem quality explained 60-92 % of variation in soil metagenome compound processing potential measures in a post-mining restoration case study dataset. With growing knowledge of the varying proficiency of environmental microbiota to process human health associated compounds, we might design environmental interventions or nature prescriptions to modulate our exposures, thereby advancing microbiota-oriented approaches to human health. Compound processing potential offers a simplified, integrative approach for applying metagenomics in ongoing efforts to understand and quantify the role of microbiota in environmental- and human-health.

Emergent community architecture despite distinct diversity in the global whale shark (Rhincodon typus) epidermal microbiome.

Microbiomes confer beneficial physiological traits to their host, but microbial diversity is inherently variable, challenging the relationship between microbes and their contribution to host health. Here, we compare the diversity and architectural complexity of the epidermal microbiome from 74 individual whale sharks (Rhincodon typus) across five aggregations globally to determine if network properties may be more indicative of the microbiome-host relationship. On the premise that microbes are expected to exhibit biogeographic patterns globally and that distantly related microbial groups can perform similar functions, we hypothesized that microbiome co-occurrence patterns would occur independently of diversity trends and that keystone microbes would vary across locations. We found that whale shark aggregation was the most important factor in discriminating taxonomic diversity patterns. Further, microbiome network architecture was similar across all aggregations, with degree distributions matching Erdos-Renyi-type networks. The microbiome-derived networks, however, display modularity indicating a definitive microbiome structure on the epidermis of whale sharks. In addition, whale sharks hosted 35 high-quality metagenome assembled genomes (MAGs) of which 25 were present from all sample locations, termed the abundant 'core'. Two main MAG groups formed, defined here as Ecogroup 1 and 2, based on the number of genes present in metabolic pathways, suggesting there are at least two important metabolic niches within the whale shark microbiome. Therefore, while variability in microbiome diversity is high, network structure and core taxa are inherent characteristics of the epidermal microbiome in whale sharks. We suggest the host-microbiome and microbe-microbe interactions that drive the self-assembly of the microbiome help support a functionally redundant abundant core and that network characteristics should be considered when linking microbiomes with host health.

Epidermal Microbiomes of Leopard Sharks (Triakis semifasciata) Are Consistent across Captive and Wild Environments.

Characterizations of shark-microbe systems in wild environments have outlined patterns of species-specific microbiomes; however, whether captivity affects these trends has yet to be determined. We used high-throughput shotgun sequencing to assess the epidermal microbiome belonging to leopard sharks (Triakis semifasciata) in captive (Birch Aquarium, La Jolla California born and held permanently in captivity), semi-captive (held in captivity for <1 year in duration and scheduled for release; Scripps Institute of Oceanography, San Diego, CA, USA) and wild environments (Moss Landing and La Jolla, CA, USA). Here, we report captive environments do not drive epidermal microbiome compositions of T. semifasciata to significantly diverge from wild counterparts as life-long captive sharks maintain a species-specific epidermal microbiome resembling those associated with semi-captive and wild populations. Major taxonomic composition shifts observed were inverse changes of top taxonomic contributors across captive duration, specifically an increase of Pseudoalteromonadaceae and consequent decrease of Pseudomonadaceae relative abundance as T. semifasciata increased duration in captive conditions. Moreover, we show captivity did not lead to significant losses in microbial α-diversity of shark epidermal communities. Finally, we present a novel association between T. semifasciata and the Muricauda genus as Metagenomes associated genomes revealed a consistent relationship across captive, semi-captive, and wild populations. Since changes in microbial communities is often associated with poor health outcomes, our report illustrates that epidermally associated microbes belonging to T. semifasciata are not suffering detrimental impacts from long or short-term captivity. Therefore, conservation programs which house sharks in aquariums are providing a healthy environment for the organisms on display. Our findings also expand on current understanding of shark epidermal microbiomes, explore the effects of ecologically different scenarios on benthic shark microbe associations, and highlight novel associations that are consistent across captive gradients.

Exploring the taxonomical and functional profile of As Burgas hot spring focusing on thermostable ß-galactosidases.

In the present study we investigate the microbial community inhabiting As Burgas geothermal spring, located in Ourense (Galicia, Spain). The approximately 23 Gbp of Illumina sequences generated for each replicate revealed a complex microbial community dominated by Bacteria in which Proteobacteria and Aquificae were the two prevalent phyla. An association between the two most prevalent genera, Thermus and Hydrogenobacter, was suggested by the relationship of their metabolism. The high relative abundance of sequences involved in the Calvin-Benson cycle and the reductive TCA cycle unveils the dominance of an autotrophic population. Important pathways from the nitrogen and sulfur cycle are potentially taking place in As Burgas hot spring. In the assembled reads, two complete ORFs matching GH2 beta-galactosidases were found. To assess their functional characterization, the two ORFs were cloned and overexpressed in E. coli. The pTsbg enzyme had activity towards o-Nitrophenyl-ß-D-galactopyranoside (ONPG) and p-Nitrophenyl-ß-D-fucopyranoside, with high thermal stability and showing maximal activity at 85 °C and pH 6, nevertheless the enzyme failed to hydrolyze lactose. The other enzyme, Tsbg, was unable to hydrolyze even ONPG or lactose. This finding highlights the challenge of finding novel active enzymes based only on their sequence.

Elasmobranch microbiomes: emerging patterns and implications for host health and ecology.

Elasmobranchs (sharks, skates and rays) are of broad ecological, economic, and societal value. These globally important fishes are experiencing sharp population declines as a result of human activity in the oceans. Research to understand elasmobranch ecology and conservation is critical and has now begun to explore the role of body-associated microbiomes in shaping elasmobranch health. Here, we review the burgeoning efforts to understand elasmobranch microbiomes, highlighting microbiome variation among gastrointestinal, oral, skin, and blood-associated niches. We identify major bacterial lineages in the microbiome, challenges to the field, key unanswered questions, and avenues for future work. We argue for prioritizing research to determine how microbiomes interact mechanistically with the unique physiology of elasmobranchs, potentially identifying roles in host immunity, disease, nutrition, and waste processing. Understanding elasmobranch-microbiome interactions is critical for predicting how sharks and rays respond to a changing ocean and for managing healthy populations in managed care.

The skin microbiome of elasmobranchs follows phylosymbiosis, but in teleost fishes, the microbiomes converge.

BACKGROUND: The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize (1) the skin microbiome will be host and clade-specific, (2) evolutionary difference in elasmobranch and teleost will correspond with a concomitant increase in host-microbiome dissimilarity, and (3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes. RESULTS: We show that the taxonomic and functional composition of the microbiomes is host-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch's average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes. CONCLUSION: Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence of extant vertebrate hosts and reveals that microbial relationships are not consistent across evolutionary timescales. Video abstract.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage.

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

The skin microbiome of the common thresher shark (Alopias vulpinus) has low taxonomic and gene function ß-diversity.

The health of sharks, like all organisms, is linked to their microbiome. At the skin interface, sharks have dermal denticles that protrude above the mucus, which may affect the types of microbes that occur here. We characterized the microbiome from the skin of the common thresher shark (Alopias vulpinus) to investigate the structure and composition of the skin microbiome. On average 618 812 (80.9% ± S.D. 0.44%) reads per metagenomic library contained open reading frames; of those, between 7.6% and 12.8% matched known protein sequences. Genera distinguishing the A. vulpinus microbiome from the water column included, Pseudoalteromonas (12.8% ± 4.7 of sequences), Erythrobacter (5. 3% ± 0.5) and Idiomarina (4.2% ± 1.2) and distinguishing gene pathways included, cobalt, zinc and cadmium resistance (2.2% ± 0.1); iron acquisition (1.2% ± 0.1) and ton/tol transport (1.3% ± 0.08). Taxonomic community overlap (100 - dissimilarity index) was greater in the skin microbiome (77.6), relative to the water column microbiome (70.6) and a reference host-associated microbiome (algae: 71.5). We conclude the A. vulpinus skin microbiome is influenced by filtering processes, including biochemical and biophysical components of the shark skin and result in a structured microbiome.

Aura-biomes are present in the water layer above coral reef benthic macro-organisms.

As coral reef habitats decline worldwide, some reefs are transitioning from coral- to algal-dominated benthos with the exact cause for this shift remaining elusive. Increases in the abundance of microbes in the water column has been correlated with an increase in coral disease and reduction in coral cover. Here we investigated how multiple reef organisms influence microbial communities in the surrounding water column. Our study consisted of a field assessment of microbial communities above replicate patches dominated by a single macro-organism. Metagenomes were constructed from 20 L of water above distinct macro-organisms, including (1) the coral Mussismilia braziliensis, (2) fleshy macroalgae (Stypopodium, Dictota and Canistrocarpus), (3) turf algae, and (4) the zoanthid Palythoa caribaeorum and were compared to the water microbes collected 3 m above the reef. Microbial genera and functional potential were annotated using MG-RAST and showed that the dominant benthic macro-organisms influence the taxa and functions of microbes in the water column surrounding them, developing a specific "aura-biome". The coral aura-biome reflected the open water column, and was associated with Synechococcus and functions suggesting oligotrophic growth, while the fleshy macroalgae aura-biome was associated with Ruegeria, Pseudomonas, and microbial functions suggesting low oxygen conditions. The turf algae aura-biome was associated with Vibrio, Flavobacterium, and functions suggesting pathogenic activity, while zoanthids were associated with Alteromonas and functions suggesting a stressful environment. Because each benthic organism has a distinct aura-biome, a change in benthic cover will change the microbial community of the water, which may lead to either the stimulation or suppression of the recruitment of benthic organisms.